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Controlling electrons at the level of elementary charge e has been demonstrated experimentally 
already in the 1980's. Ever since, producing an electrical current ef, or its integer multiple, at a 
drive frequency / has been in a focus of research for metrological purposes. In this review we first 
discuss the generic physical phenomena and technical constraints that influence charge transport. 
We then present the broad variety of proposed realizations. Some of them have already proven 
experimentally to nearly fulfill the demanding needs, in terms of transfer errors and transfer 
rate, of quantum metrology of electrical quantities, whereas some others are currently "just" wild 
ideas, still often potentially competitive if technical constraints can be lifted. We also discuss the 
important issues of read-out of single-electron events and potential error correction schemes based 
on them. Finally, we give an account of the status of single-electron current sources in the bigger 
framework of electric quantum standards and of the future international SI system of units, and 
briefly discuss the applications and uses of single-electron devices outside the metrological context. 
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I. INTRODUCTION 

Future definition of the ampere is foreseen to be based 
on fixing the elementary charge e. Its most direct re- 
alization would be the transport of a known number of 
electrons. Over the past quarter of a century, we have 
witnessed progress towards ever better control of individ- 
ual electrons. Since single-electron tunneling is by now 
a well established subject, several reviews of its different 



aspects exist in literature (Averin and Likharev 1991 



Averin and Nazarov 1992a[ |Durrani, 2009; Sohn et al. 



1997 van der Wiel et all 



|2002j). 



Several milestones have been achieved in the progress 
towards a single-electron current source since the ini- 
tial proposals of the metrological triangle in the mid 
80's (Likharev and Zorin 1985). The single-electron am- 



pere would be based on transporting an electron with 
charge e, or rather a known number N of electrons, Ne, 
in each operation of a control parameter that is cycli- 
cally repeated at frequency /, so that the output dc 
current is ideally equal to Nef. The needs of preci- 
sion metrology generally state that this operation has 
to be performed at a relative error level not larger than 
10~ 8 and at the same time the current level needs to 



be several hundreds of pico- amperes (Feltin and Pique 



mal 2009). Just a few years after the initial theoreti- 



cal proposals, the first metallic (Geerligs et al. 1990b 



Pothier 1991 Pothier et al.\ 1992 ) and semiconduct 



ing ( |Kouwenhoven et al. 1991b) single-electron turnstiles 
and pumps demonstrated currents / = Nef with an er- 
ror of a few percent, still orders of magnitude away from 
what is needed. Like often in precision metrology, the 
pursuit of higher accuracy has been a pacemaker for un- 
derstanding new physics, since the errors that need to 
be suppressed are often a result of interesting physical 
phenomena. For instance, quantum multi-electron pro- 
cesses and non-equilibrium phenomena have been inten- 
sively studied in order to improve the performance of 
single-electron sources. In five years, the accuracy of 
single-electron pumps was remarkably improved by an- 



other five to six orders of magnitude (Keller et al. 1996) 



by effectively suppressing the so-called co-tunneling cur- 
rent, however, at the expense of significantly increased 
complexity of the device and reduced overall magnitude 
of the output current (a few pico-amperes) of the pump. 



Alternative ideas were to be found. At the same time, 
single-electron conveyors in semiconducting channels us- 
ing surface-acoustic wave driving yielded promising re- 
sults, in particular in terms of significantly increased 
current level ( |Shilton et al. 1996b). Yet likely due to 



overheating effects in the channel, it may turn out to be 
difficult to suppress thermal errors to the desired level 
using this technique. 

Interestingly there was a decade of reduced progress in 
the field, until in the second half of the previous decade 
several new proposals and implementations were put for- 
ward. The most promising of these devices are unde- 
niably the sources based on a quantum dot (Blumen- 



(Kaestner et al. 



thai et aZ.[|2007[) , down to a single-parameter ac control 

and a SINIS turnstile (Pekola 



2008b) 



et al.\ 2008), which is a basic single-electron transis- 



tor with superconducting leads and normal metal island. 
These simple devices promise high accuracy, and a pos- 



sibility to run many of them in parallel (Maisi et al. 



2009). At around the same time, other promising ideas 



came out, for example a quantum phase slip based super- 
conducting current standard (Mooij and Nazarov, 2006)). 
Quantum phase slips provide the mechanism for the ex- 
istence of the Coulomb blockade effects in superconduct- 



ing wires without tunnel barriers (Astafiev et al. 2012), 



, and could potentially lead to current standards produc- 
ing larger currents. At the time of writing this review, 
we are definitely witnessing a period of intense activity 
in the field in a well-founded atmosphere of optimism. 



II. PRINCIPLES OF MANIPULATING INDIVIDUAL 
ELECTRONS 

A. Charge quantization on mesoscopic conductors 

In this section, we briefly summarize the essential con- 
cepts of single-electron device physics, with the empha- 
sis on the topics needed for the subsequent discussion 
of the quantized current sources. We focus mostly on 
metallic devices since those have an elaborated theory 
based in first principles. The main concepts can be how- 
ever adapted when considering semiconducting systems 



as well (Zimmerman et al. 2004) 



As is well known from the elementary treatments of the 
Bohr's model in quantum mechanics, electrostatic energy 
of an electron in a hydrogen atom is roughly equal to the 
kinetic energy of its confinement in the atomic orbital of 
the radius. The fact that the characteristic energy sep- 
aration of levels in the confinement energy spectrum de- 
creases much more rapidly than the electrostatic energy 
with the size of the confining region ensures then that 
in mesoscopic conductors which are large on the atomic 
scale, the electrostatic energy of individual electrons can 
be large even in the regime where the separation of the 
individual energy levels associated with quantum con- 
finement of electrons is negligible. As a characteristic 
estimate, electrostatic energy of charge e of one elec- 
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tron on a micrometer-size conductor is on the order of 
a milli-electron volt, or 10 K in temperature units, and is 
many orders of magnitude larger than the energy sepa- 
ration SE of electron confinement levels in the same con- 
ductor, which should be about nano-electron volt, well 
below all practical temperatures. As a result, at low 
but easily reachable temperatures in the kelvin and sub- 
kelvin range, the properties of mesoscopic conducting is- 
lands are dominated by electrostatic energy of individual 
electrons, while small SE provides one of the conditions 
that makes it possible to use macroscopic capacitances to 
quantitatively describe electrostatics of these conductors 
even in this "single-electron" regime. The charging en- 
ergy U of a system of such conductors can be expressed 
then as usual in terms of the numbers rij of excess elec- 
trons charging each conductor and the capacitance ma- 
trix C, see, e.g., (Landau and Lifshitz 1980a) 




FIG. 1 Schematic diagram of the basic circuit for manipulat- 
ing individual electrons, single-electron box (SEB): conduct- 
ing island carrying electric charge en, and an electrostatically 
coupled external electrode with the charge e(N — n) producing 
the gate voltage V g . 



U{{n j }) = e ~YjP~ l ^ n < 



(1) 



where the sum runs over all conductors in the structure. 

The electrostatic energy creates energy gaps sep- 
arating different charge configurations {nj} which pro- 
vide the possibility to distinguish and manipulate these 
charge configurations. Historically, one of the first obser- 
vations of distinct individual electron charges were the 
Millikan's experiments on motion of charged micrometer- 
scale droplets of oil, which produced the evidence that 
"all electrical charges, however produced, are exact mul- 
tiples of one definite, elementary, electrical charge" (Mil- 
likan 1911). In those experiments, the oil droplets were, 



however, charged randomly by an uncontrollable process 
of absorption of ions which exist normally in air. By con- 
trast, in mesoscopic conductors, the charge states rij can 
be changed in a controllable way. Besides the charging 
energy ([!), such a process of controlled manipulation of 
individual charges in mesoscopic conductors requires two 
additional elements. First, the tunnel junctions formed 
between the nearest-neighbor electrodes of the structure 
enable the electron transfer between these electrodes, and 
second, the possibility to control the electrostatic en- 
ergy g aps by continuous variation of charges on the junc- 
tions ( Averin and Likharev 1986 ) . The simplest way of 



varying the charges on the tunnel junctions continuously 
is by placing the electrodes in external electrical fields 
( Biittiker 1987 ) that create continuously- varying poten- 



tial differences between the electrodes of the structure. 
Externally-controlled continuous transport and gate volt- 
ages produced in this way can be used then to transfer 
individual electrons in the system of mesoscopic conduc- 
tors. 

A simple model of the sources of continuously-varying 
external voltages is obtained by taking some of the elec- 
trodes of the structure described by the energy Q to 
have very large self-capacitance and carry large charge, 
so that the tunneling of a few electrons does not affect 
the potentials created by them. For instance, the most 



basic single-electron structure, single-electron box (SEB) 



(Lafarge et al. 1991), can be simplified to the two elec- 



trodes, one main island carrying the charge en, and the 
electrode with the charge e(N — n) creating the gate volt- 
age V g (Fig.JTJa)). Quantitatively, the structure in Fig.[l] 
is characterized by the capacitance matrix 



C 



rig /-< 

-Cg C S 



(2) 



where Cf., Cs > C g are the total capacitances of the 
gate electrode and the island, respectively. In the limit 
N, C$ ->■ oo, with eN/Cft = V g , the energy {l) of the 
charges shown in Fig. [T] reduces for the capacitance ma- 
trix |2]) to the following expression: 



U = U Q + E c n l 



(3) 



Above, Uq is an n-independent energy of creating the 
source of the gate voltage, Uq — e 2 N 2 /2Cq in this case, 
Eq = e 2 /2Cs is the charging energy of one electron on 
the main electrode of the box, and en g = C g V g is the 
charge induced on this electrode by the gate voltage V g 
through the gate capacitance C g . As one can see from 
Eq. pi), the gate voltage V g indeed controls the energy 
gaps separating the different charge states n of the main 
island and therefore makes it possible to manipulate in- 
dividual electron transitions changing the island charge 
en. 

Figure[2ja) shows a scanning electron micrograph of an 
example of a realistic box structure, in which, in contrast 
to the schematic diagram of Fig. [T] one pays attention to 
satisfying several quantitative requirements on the box 
parameters. First of all, the capacitance Cs needs to be 
sufficiently small to have significant charging energy Ec, 
while one keeps the gate capacitance C g not very small in 
comparison to Cs, to be able to manipulate the charge en 
more easily and also to measure it. The box in Fig. [2](a) 
consists of the similar-size islands, and its equivalent elec- 
tric circuit is shown in Fig. [2jb) . The charging energy 



of the box is described by the same expression (J3|) , with 
en being the charge transferred from the left to the right 
island, and - the total mutual capacitance between 
the two islands, C s = C + C g , where Cg 1 = C^ 1 +C R 1 . 
Connecting the box islands to the source of gate voltage 
V g through capacitances Cl,r on both sides reduces cou- 
pling to parasitic voltage fluctuations in the electrodes of 
the structure, responsible for environment-induced tun- 
neling discussed below. Generally, a practical geometric 
structure of the box islands is also determined by the 
fact that the main contribution to the capacitance 
comes from the tunnel junction formed in the area where 
the "arms" of the islands [Fig. [2ja)] overlap. The size 
of this area should be minimized to increase Eq- At 
the same time, the islands themselves can be made much 
larger than the junctions, to increase the gate capaci- 
tance C g without affecting strongly the total capacitance 
Ce- Besides increasing the coupling to the gate voltage 
created by the two outside horizontal electrodes in the 
box structure shown in Fig. [2ja) , larger size of the box 
islands also increases the coupling to the single-electron 
transistor (discussed in more details below) which mea- 
sures the charge of the box and can be seen in the upper 
right corner of Fig.[2ja). 

As was mentioned above, the main qualitative prop- 
erty of the SEB is that it allows one to manipulate in- 
dividual electrons through variation of the gate voltage 
V g . Indeed, at low temperatures, T <C E c jk Bl the box 
occupies the ground state of the charging energy ^ . For 
a given gate- voltage-induced charge n g , the minimum is 
achieved when the number n of extra electrons on the 
island equals n g rounded to the nearest integer. This de- 
pendence of n on n g means that one electron is added or 
removed from the box island, changing n by ±1, when- 



ever n g passes through a degeneracy point, i. e. 



1/2 



modulo an integer, at which point the charging energies 
([3]) of the two charge states that differ by one electron 
transition, 8n = 1, are equal. If the gate voltage increases 
monotonically, the de pendence n(n Q ) has the shape of the 
"Coulomb staircase" (Lafarge et al. 1991 ), with each step 



of the staircase corresponding to the addition of one elec- 
tron with gate voltage increase by 8n g = 1. If the gate 
voltage oscillates in time around the degeneracy point 
n g = 1/2, as in Fig.[2[c), with an appropriate amplitude 
(Srig ~ 1), it inducesback-and-forth electron transitions 
between the two charge states separated by one electron 
charge, 6n = 1, which can be seen in Fig. [2fc) as the 
two-level telegraph signal of the detector measuring the 
box charge. Thus, Fig. [2^c) gives a practical example of 
manipulation of an individual electron transition in the 
SEB. 

One of the most interesting dynamic manifestations of 
the manipulation of individual electrons in a system of 
mesoscopic conductors is the possibility to arrange the 
system dynamics in such a way that electrons are trans- 
ferred through it one by one, in a correlated fashion. This 
can be achieved, for instance, if the gate voltage V g of the 
SEB grows in time at a constant rate such that effectively 
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FIG. 2 Practical SEB: (a) scanning electron micrograph of 
a realistic box structure, (b) its equivalent electric circuit, 
and (c) single-electron transitions in the box illustrating the 
"charge quantization" : a time-dependent gate- voltage V g (t) 
(sinusoidal curve) of an appropriate amplitude drives individ- 
ual electron transitions changing the box state between the 
two discrete charge configurations, electron on the left or on 
the right island. These two charge states are detected via the 
shown two-level detector output current synchronous with the 
osci llating V q {t) . (Figur e adapted from Refs. Saira et al. 2010 
and Saira et al. 2012b I. 



a constant dc "displacement" current I = en g is injected 
in the box junction. The same dynamic would be ob- 
tained if the real dc current / flows into a mesoscopic 
tunnel junction. In this case, correlated successive trans- 
fer of electrons one by one through the junctio n gives 
rise to the "single-electron tunneling" oscillations ( Averin 



and Likharev 1986| Bylander et al. 2005 ) of voltage on 
the junction, dU/d(en) = e(n — n g )/Cs, with frequency 
/ related to the current by the fundamental equation 

I = ef. (4) 

More complex structures than SEB or an individual tun- 



nel jun ction, like s ingle-electron turnstile ( Geerligs et al. 
|1990bb and pump (Keller et aUI 19961 IPothier et a/.l|1992 



discussed below, make it possible to "invert" this relation 
and transfer one electron per period of the applied gate 
voltage oscillation with frequency /. The above discus- 
sion of the manipulation of individual electrons in the 
SEB shows that the charge states n, while controlled by 
the gate voltage V g , remain the same in a range of vari- 
ation of V g . Physically, such "quantization of charge" 
results from the fact that an isolated conductor can con- 
tain only an integer total number of electrons, with the 
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charging energy producing energy gaps separating differ- 
ent electron number states. Charge quantization enables 
one to make the accuracy of manipulation of individual 
electron charges in structures like SEB very high, in prin- 
ciple approaching the metrological level. Such potentially 
metrological accuracy also extends to the transport in 
turnstiles and pumps, making the current sources based 
on single-electron tunneling promising candidates for cre- 
ation of the quantum standard of electrical current. 




B. Sequential single-electron tunneling 

As discussed above, one of the key elements in manipu- 
lating individual electrons in systems of mesoscopic con- 
ductors is a tunnel junction, which provides the means to 
transfer electrons along the system thus creating the dc 
current / through it. A tunnel junction ( Giaever 1960 ) is 



a system of two conductors separated by a layer of insu- 
lator that is sufficiently thin to allow electrons to tunnel 
between the conductors (see Fig. [3]). For normal conduc- 
tors, the current through the junction depends linearly on 
the applied transport voltage, and is characterized by the 
tunnel conductance Gt = 1/-Rt- I n the single-electron 
devices, Gt should satisfy two contrasting requirements. 
To increase the current / driven through the structure, 
e.g., to increase the allowed range of frequencies / for 
which Eq. Q is satisfied accurately, one should maxi- 
mize Gt- On the other hand, charge quantization on the 
electrodes of the structure requires that they are well- 
isolated from each other, i.e., Gt should be small. The 
latter condition can be formulated more quantitatively 
requiring that the characteristic charging energy Ec of 
the localized charge states is well-defined despite the fi- 
nite lifetime of these states, ~ Gt / C, where C is the typi- 
cal junction capacitance in the structure: Ec 3> %Gt/C. 
This condition can be expressed as Gt <C 1/Rq, where 
Rq = irh/e 2 ~ 13.9 kf2 is the characteristic "quantum" 
resistance. When this condition is satisfied, the local- 
ized charge states provide an appropriate starting point 
for the description of a single-electron structure, while 
electron tunneling can be treated as perturbation. In 
what follows, we mostly concentrate on such a regime of 
"strong Coulomb blockade" which is necessary for imple- 
mentation of precise transport of individual electrons, as 
required for quantized current sources. 

Majority of practical metallic structures employ tun- 
nel junctions based on barriers formed by either ther- 
mal or plasma oxidation of aluminum. The main reason 
for this are the superior properties of the aluminum ox- 
ide layer, in terms of its uniformity and electrical and 
noise properties. A typical barrier structure is shown 
in Fig. [3] that includes a high resolution transmission- 
electron-microscopy (TEM) image of a cross-section of an 
aluminum-based junction with amorphous ALjOj, tunnel 
barrier. From the point of view of the Landauer-Biittiker 
formula for electric conductance of a mesoscopic conduc- 
tor, the junction tunnel conductance can be expressed 




FIG. 3 High resolution TEM image of a cross-section of an 



aluminium oxide tunnel junction (from Ref. Prunnila et al. 
20101. 



as Gt = (1/ ' RQ)^2jTj 1 where the sum is taken over 
spin-degenerate electron transport channels propagating 
across the junction, and Tj is the quantum mechanical 
transmission coefficient of the insulator barrier for elec- 
trons in the jth channel. Condition of strong Coulomb 
blockade, Gt <C 1/Rq, implies that all individual trans- 
mission coefficients are small, Tj <C 1. Although the 
transmission coefficients Tj are sensitive to the atomic- 
scale structure of the junction, the fact that the alu- 
minum oxide layer is relatively uniform on an intermedi- 



ate space scale larger than the individual atoms ( Greibe 
et al.\ |2011[ ) allows transport properties to be estimated 
semi-quantitatively from the "bulk" properties of the 
barrier. 

Since the tunnel current depends exponentially on 
the barrier parameters, the measured electron tunneling 
rates in high-resistance junctions and in the large volt- 
age range allow one to estimate parameters of the alu- 
minum oxide barrier (see, e.g., Tan et al. 2008]): they 
yield a barrier height U ~ 2 eV, and effective electron 
mass m e ff ~ 0.5 m e in terms of the free electron mass 
m e . While the dimensions of the typical tunnel junc- 
tions need to be small - of the order of ~ 100 nm, also 
cf. Fig. M&) - in order to make the junction capacitance 
sufficiently low, they are still quite large on the atomic 
scale. In this regime, discreteness of the spectrum of the 
transverse modes j is negligible, and the tunnel conduc- 
tance Gt is proportional to the junction area A. For the 
value of specific junction resistance A/Gt ~ 10 kil /im 2 
typical for the tunnel junctions, estimates using the bar- 
rier parameters and the simplest assumption of ballistic 
transport in the junction give for the barrier transparency 
T ~ 10~ 6 corresponding to barrier thickness close to 2 
nm (cf. Fig. [3]). Barrier with this thickness transmits ef- 
fectively only the electrons impinging on it orthogonally. 
This "focusing" effect means that the tunnel conductance 
can be expressed in terms of one maximum value of the 
transmission coefficient, Gt = NT/Rq, where the effec- 



G 



tive number N of the transport channels in the junction 
is not determined directly by the density of states in the 
electrodes, but depends also on the characteristic "traver- 
sal energy" eo of the barrier which gives the energy scale 
on which the barrier transparency changes with energy, 
N ~ Ameo/2Trfi 2 . For the parameters of the aluminum 
oxide barrier mentioned above, this gives for the area per 
transport channel, A/N ~ 1 nm 2 . As will be discussed in 
the next subsection, some of the higher-order transitions 
in the single-electron structures, e.g. Andreev reflection, 
depend explicitly on the barrier transmission coefficients 
Tj. In contrast to this, the lowest-order electron tun- 
neling depends only on the total junction conductance 
Gt- 

The most straightforward approach to the description 
of tunneling in the single-electron structures in the strong 
Coulomb blockade regime is based on the tunnel Hamil- 



tonian method (M. and C. 1962), in which the junction 



is modeled with the following Hamiltonian: 



H — Hi -\- H2 ~\~ Ht , Hj 



+ /i.e.] . (5) 



Here i?i ( 2 are the Hamiltonians of the junction elec- 
trodes, Ht is responsible for tunneling, with Ck and c p 
denoting electron destruction operator in the two elec- 
trodes, respectively, and T^ v are the tunneling ampli- 
tudes. As was discussed in the previous section, dis- 
creteness 5E of the single-particle electron states is com- 
pletely negligible in a typical metallic mesoscopic con- 
ductor, so these states form a continuum with some den- 
sity of states which in a normal metal, is constant on a 
small energy scale of interest in the structures. In this 
case, one can treat Ht using Fermi's golden rule to ob- 
tain the rate T(E) of a tunneling process that changes 
the charge configuration {rij} on the system of meso- 
scopic conductors by transferring one electron through a 
tunnel junction between the two conductors. As a re- 
sult, changing the electrostatic energy ([I]) by an amount 
E = ^({ftj,in}) — ^({ n j,fin}), where {nj t i n } is the initial 
and {rijfi n } the final charge configuration, we obtain 



T(E) = 



Gn 



def(e)[l - f(e + E^ejv^e + E) . (6) 



In this expression, /(e) is the equilibrium Fermi distri- 
bution function, and v^(e) is the density of the single- 
particle states in the jth electrode of the junction, 
j = 1,2, in units of the normal density of states pj, 
which together with the average of the squares of the 
tunneling amplitudes determine the tunnel conductance, 
Gt = 47re 2 (|Tfc p | 2 )pi j 02- Equation (16J) assumes that the 
energy E ~ Ec is much smaller than all internal en- 
ergies of the junction in the normal state, in particular 
the traversal energy eo, the condition very well satisfied 
for practical metallic structures in which Ec ~ 1 meV, 
while eo ~ 1 eV. Using the standard properties of the 
Fermi distribution functions, one can see directly that the 
rate (|6| of tunneling between the two equilibrium elec- 
trodes satisfies the necessary detailed balance condition 



Tx->2(-E) = e- E / kBT T 2 ^i{E). If, in addition, the den- 
sities of states are symmetric with respect to the chemi- 
cal potentials of the electrodes, the tunneling rate is also 
symmetric, Fi_j.2(^) = ^2^i{E), and the detailed bal- 



ance condition simplifies to T(—E) = e 



_ p -E/k B T 



T(E). The 



detailed balance condition makes it possible to express 
the tunneling rate (|6| in terms of the current-voltage 
characteristic I(V) of the junction at fixed bias voltage 
V: 



r(E) = I(E/e)/[e(l-e 



-E/k B T 



)] 



(7) 



For the NIN junctions, when both electrodes are in the 
normal (N) states, Vj{t) = 1, Eq. (|6| gives, in agreement 
with Eq. ([7]), for the tunneling rate 



T{E) 



Gn 



E 



1 



-E/k B T 



(8) 



Tunneling of individual electrons with the rate ([8]) is an 
irreversible dissipative process which converts the electro- 
static energy change E into internal energy of the elec- 
tron gas inside the junction electrodes. In accordance 
with this understanding, at small temperatures T, the 
rate ^ vanishes as e E / kBT for energetically unfavorable 
transitions, E < 0, when the energy for the transition is 
taken from the thermal fluctuations of the electron reser- 
voirs. In the regime of allowed transitions, E ~ Ec > 0, 
the magnitude of the typical transition rate T ~ Gt/C 
for the realistic values of the parameters, Gt ~ 1 MSI, 
C ~ 10~ 16 . . . 1CT 15 F, is quite high, in the GHz range. 

In the SIN junctions, when one of the junction elec- 
trodes is a superconductor (S), the BCS density of states: 
Vi (e) = |e|/(e 2 - A 2 ) 1 / 2 for |e| > A, and vanishing oth- 
erwise, implies that at temperatures well below the su- 
perconducting energy gap A, the tunneling rate ^ is 
strongly suppressed and can be reduced into the kHz and 
even Hz range. Indeed, evaluating the integral in Eq. ^ 
for the SIN junction assuming fcgT, E <C A, one gets 



G 



T(E) = -^-y/2nAk B Te- 



A/k B T sinh{E/k B T) 

I _ e -E/k B T 



(9) 



Figure [4] shows the tunneling rate ([9]) measured in an 
SIN junction in the configuration of a "hybrid" SEB [see 
Fig. Eta)], in which one of the islands of the box is a 
superconductor (aluminum), the other one being normal 
metal (copper). The electrostatic energy change E in 
the case of the box follows from Eq. ^ as E = U(n = 
0) — U(n = 1) = 2Ec{n g — 1/2), is proportional to the 
deviation of the gate voltage of the box from the degener- 
acy point n g = 1/2. The measurements can be described 
well by Eq. ^ with reasonable values of parameters in- 
cluding the superconducting energy gap A of aluminum. 

Since the tunneling transitions described quantita- 
tively by the rates ([6|-([9| are inherently random stochas- 
tic processes, dynamics of the structures in the strong 
Coulomb blockade regime and electron transport prop- 
erties including the dc current /, current noise, or even 
full statistical distribution of the transferred charge, can 
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FIG. 4 (a) Measured thermally activated rates of forward 
T(E) and backward T(-E) tunneling in a "hybrid" SIN 
single-electron box at different temperatures, as a function 
of the gate voltage offset from the degeneracy point related 
to the energy change E in tunneling as E = 2Ec(n g — 1/2). 
Solid lines are the theory prediction according to Eq. (JoJ) with 
fitted parameters: E c = 157 ^eV, A = 218 fj,eV, 1/G T = 100 
Mfi. (b) The tunneling rate at degeneracy, E = 0, as a func- 
tion of temperature (square markers), and best fit (solid line) 
to Eq. (Figure adapted from Ref. Saira et al. 2012b ) 



be obtained from the time evolution of the probabilities 
P({ n j}) °f various charge configurations {rij} governed 
by the standard rate equation for the balance of the prob- 
ability fluxes. The most basic single-electron system that 
allows for the flow of dc current through it, and gives an 
example of such an equation is the single-electron transis- 



tor (SET) ( Averin and Likharev 


198G 


Fulton and Dolan 


1987 |Likharev 


1987 


). The transistor can be viewed as a 



generalization of the SEB, and consists of a mesoscopic 
conducting island connected by two tunnel junctions to 
the bulk electrodes that provide the transport voltage V 
across it. The island is also coupled capacitively to the 
source of the gate voltage V g which controls the flow of 
current / through the transistor between the two elec- 
trodes. Equivalent circuit of the transistor is shown in 
Fig. [5j and an example of its geometric structure can be 
seen in the upper right corner of Fig. [2ja), where it is 
used to measure the charge state of the SEB. The charge 
configuration of the transistor is characterized simply by 
the number n of extra electrons on its central island, and 
accordingly, the rate equation describing its dynamics is: 



p(n) 



£[p(n±l)r 



( + ) 
3 



(nil) 



■p(n)Tf\n)\, (10) 



where p(n) is the probability distribution of the charge 
en on the central island of the transistor, and the rates 
Tj(n) describe the tunneling processes in junction j 
with the tunnel conductance Gj out of the state n in 
the direction that increases (+) or decreases (— ) n by 
one. The rates are given by Eqs. pi or d9j), or their gen- 



eralizations, depending on the nature of the transistor 
electrodes. They depend on the indices of T's in Eq. ( 10 ) 



through the change of E of the charging energy U of the 
transistor, which is a function of all these indices. The 
transistor energy U consists of two parts, one that coin- 
cides with the charging energy ^ of the SEB in which 
Cs = C\ + C 2 + C g , and the other, Uy, is created by the 
transport voltage V: 



U v = -eNV - enV(C 2 + C g /2)/C T 



(11) 



Here iV is the number of electrons that have been trans- 
ferred through the transistor. Both the dc current I 
through the transistor (Averin and Likharev 1991[ ) and 
the current noise (Korotkov 1994) can be calculated 



starting from Eq. (10). The main physical property 
of the transistor transport characteristics is that they 
depend periodically on the gate voltage, in particular 
I(n g + 1) = I (rig). This dependence of the transistor 
current on the charge en g induced on its central island 
makes the SET a charge detector, with sub-el ectron sen- 
sitivity approaching (10 -5 . . . 10 -6 ) e/Hz 1 / 2 (Krupenin 



1998 Roschier et al. 



2001 Zimmerli et al. 1992). As 



a result, the SET is the most standard charge detector 
for measurements of, e.g., individual electron dynamics 
in other single-electron structures [cf. Fig. 2jja)]. 

The hybrid SINIS or NISIN transistors have an addi- 
tional important feature that distinguishes them from the 
SETs with normal electrodes. They provide the possibil- 
ity to realize the regime of the quantized current I Q, 
when driven by the ac gate voltage V g (t) of frequency / 
(IPekola et a/. 1 120081). This property of the hybrid SETs 



is one of the main topics of this review discussed below. 



C 1' G 1 




FIG. 5 Equivalent electric circuit of an SET. 

The basic expression ^ for the tunneling rates as- 
sumes that the electrodes of the tunnel junction are in 
equilibrium at temperature T, with the implied assump- 
tion that this temperature coincides with fixed tempera- 
ture of the whole sample. Since each electron tunneling 
event deposits an amount of heat AU into the electron 
system of the electrodes, this condition requires that the 



relaxation processes in the electrodes are sufficiently ef- 
fective to maintain the equilibrium. The relaxation rates 
decrease rapidly with decreasing temperature, e.g., oc T 5 
for electron-phonon relaxation in an ordinary metal - see, 
e.g. ( |Giazotto et al. 2006 ). This makes the relaxation in- 
sufficient and causes the overheating effects to appear at 
some low temperature, in practice around 0.1 K. There- 
fore, the overheating sets a lower limit to the effective 
temperature of the transitions, in this way limiting the 
accuracy of control over the individual electron transport. 

One more assumption underlying the expression ^ for 
the tunneling rate in the SIN junctions is that the elec- 
tron distribution function is given by the Fermi distri- 
bution function /(e). As known from the statistical me- 
chanics, even in equilibrium, this requires that the total 
effective number of particles that participate in forming 
this distribution is large. In normal-metal islands, this 
requirement is satisfied at temperatures much larger than 
the single-particle level spacing, T ^> SE/ks, as is the 
case for practically all metallic tunnel junctions. In con- 
trast to this, in superconducting islands, this condition 
can be violated at temperatures below the superconduct- 
ing energy gap, T <C A/ks, when the total number of 
the quasiparticle excitations in the electrode is no longer 
large. The temperature scale of the onset of this "indi- 
vidual quasiparticle" regime can be estimated from the 
relation J~ def(e)v(e)/8E ~ 1. The main qualitative 
feature of this regime is the sensitivity of the electron 
transport properties of a superconducting island to the 



parity of the total number of electrons on it ( Averin and 
Nazarov 1992b |Tuominen et aL[|1992[ ). In particular, 
the charge tunneling rate (|9| in the SIN junction should 
be modified in this case into the rates of tunneling of in- 
dividual quasiparticles. For T 3> SE/ks, these rates are 
still determined by the same average {\Tj~ p \ 2 ) over many 
single-particle states of the squares of the tunneling am- 
plitudes |5| which give the tunnel conductance Gt, and 
therefore can be expressed through Gt- In the regime of 
"strong" parity effects, when T < T* = A/k B \nN eff , 
where N eff = {2irk B T A) 1 ' 2 / 5E is the effective^ num- 



ber of states for the quasiparticle excitations (Tuomi 



nen et al. 19921, an ideal BCS superconductor should 
reach the state with no quasiparticles, if the total number 
Nq of electrons in the superconductor is even, and pre- 
cisely one unpaired quasiparticle if N is odd. Although 
many non-equilibrium processes in realistic superconduc- 
tors lead to creation of a finite density of quasiparticle ex- 
citations which do not "freeze out" at low temperatures 
(see, e.g., (de Visser et al. 2011 )), one can realize the sit- 



uation wit 



i the number of quasiparticles controlled as in 



an ideal BCS superconductor, as e.g. in i 


Lafarge et al. 


1993 


Saira et al. 


2012a 


Tuominen et al. 


1992 


). In this 



regime, the rates of sequential charge tunneling between 
the normal metal electrode and a superconducting island 
depend on th e parity of the total numb er iVo of electrons 
in the island flSchon and Zaikin[|l994|). For T < T* , the 



tunneling rates both to and from the island are domi- 
nated by the one quasiparticle that exists on the island 



for odd Nq. When this quasiparticle is equilibrated to 
the edge of the quasiparticle spectrum at energy A, the 
rates of tunneling to and from the island (i.e., increasing 
and decreasing the charge en of the island) coincide, and 
for \E\ -C A are independent of the electrostatic energy 
change E: 



r 



odd 



G T SE 
4e 2 



(12) 



For even No, when there are no quasiparticles on the 
island, tunneling necessarily involves the process of cre- 
ation of a quasiparticle making the tunneling rates energy 
dependent on the energy change E: 



T odd N -(A-E)/k B T 



(13) 



In the hybrid superconductor/normal metal structures, 
these tunneling rates determine the electron transport 
properties through the rate equation (10). 



C. Co-tunneling, Andreev reflection, and other higher-order 
processes 

Sequential tunneling discussed above represents only 
the first non-vanishing order of the perturbation the- 
ory in the tunnel Hamiltonian In the strong 
Coulomb blockade regime Gt *C 1/Rq, this approxima- 
tion provides an excellent starting point for the descrip- 
tion of electron transport, accounting quantitatively for 
main observed properties of these structures. However, 
more detailed picture of the transport should include also 
the tunneling processes of higher order in Ht, which in- 
volve transfers of more than one electron in one or several 
tunnel junctions. Although for Gt <C 1/ Rq, the rates of 
these more complex multi-step electron " co- tunneling" 
processes are small in comparison with the rates of the 
single-step sequential electron tunneling, they are fre- 
quently important either because they provide the only 
energetically allowed transport mechanism, or because 
they limit the accuracy of the control of the basic se- 
quential single-electron transitions. The simplest exam- 
ple of the co-tunneling is the current leakage in the SET 
in the CB regime (Averin and Odintsov 1989 Geerligs 
et al. 1990a[ ), when the bias voltage V is smaller than the 
CB threshold and any single-step electron transfer that 
changes the charge en on the transistor island by ±e (see 
Fig. [5| would increase the charging energy ^ and is sup- 
pressed. In this regime, only the two-step co-tunneling 
process that consists of electron transfers in both junc- 
tions of the transistor in the same direction gains the bias 
energy (11). It achieves this by changing the number N 



of electrons transported through the transistor by 1 with- 
out changing the charge en on the island. Qualitatively, 
this process represents a quantum tunneling through the 
energy barrier created by the charging energy. Because of 
the discrete nature of charge transfer in each step of the 
co-tunneling, its rate is not suppressed exponentially as 
for the usual quantum tunneling, and is smaller only by 



9 



a factor GtRq <C 1 than the rate of sequential tunneling 
processes. 

In a hybrid SIN junction, in addition to the charging 
energy, superconducting energy gap A provides an extra 
energy barrier to tunneling of individual electrons, sup- 
pressing the sequential tunneling rate Q at low temper- 
atures T <C A/fes- The gap A exists only for individual 
electrons, while pairs of electrons with zero total energy 
and momentum can enter a superconductor as a Cooper 
pair, in the process called Andreev reflection (Andreev 



superconductor obtaining the total AR tunneling rate as 
T AR = 1 rJ2\ A ( e ''^l)\ 2 f(^)f^l)S(e k +e l +E), (15) 



k.l 



where E is the electrostatic energy change due to the 
complete AR tunneling process. The sum over all states 



p in the superconductor in Eq. ( 14 ) implies that the con 



1964). In tunnel junctions, Andreev reflection (AR) can 



be described similarly to the co-tunneling, as a pertur- 
bative two-step tunneling process, in which the transfer 
of the first electron is virtual and only the second elec- 
tron transfer makes the process energetically favorable 
and real. Quantitatively, the rates of such multi-step 
transitions can be determined through their higher-order 
transition amplitudes constructed according to the stan- 
dard rules of the perturbation theory (see, e.g., (Landau 



tribution of the individual quasiparticles (which is impor- 
tant in the parity-dependent transition rates ( 12 ) and 
(13)) is negligible in the amplitude A, and individual 
quasiparticles affect Tab. only through the change of the 
charging conditions of tunneling. 

Result of summation over different single-particle 
states in Eqs. (14) and (15) depends on the geometric 
structure of the SIN junction. For instance, quadratic 
dependence of the AR amplitude A ( 14 ) on the tunncl- 



and Lifshitz 1980b)). For instance, in the simplest exam- 
ple of a two-step AR process in a hybrid single-electron 
box, the elementary amplitude A(e k ,€i) of the process 
that takes two electrons in the normal electrode with en- 
ergies efc, ei and transfers them into the superconductor 
as a Cooper pair can be written as: 



ing amplitudes makes the magnitude of Andreev reflec- 
tion sensitive not only to the total tunnel conductance 
Gt but also to the distribution of the barrier transmis- 
sion probabilities. Two main qualitative features of the 
aluminium oxide tunnel junctions (Fig. [3]), which are the 
focus of the main part of this review, are the relatively 
thick insulator barrier characterized by the focusing ef- 
fect on the tunneling electrons, and low resistance of the 
junction electrodes. The simplest junction model that 
takes into account both features assumes ballistic elec- 
tron motion that can be separated into different trans- 
port channels throughout the junction. In this case, the 



A{e kl ei) 



p 



Up Vp T k p T[p( 



+ Ei — e/g fl p + Ei 



states k, I in Eqs. ( 15 ) belong to the same transport chan- 
nel, and summation over different channels can be done 
directly and gives the effective number TV of the chan- 
nels which, as discussed in the section |II.B| is limited 
by the angular dependence of the barrier transmission 
probabilities (see, e.g., 



ballistic approximation, Eqs 



Averin and Bardas 
Wand (|15 



The two-step process goes through an intermediate state 
obtained as a result of the first step of the process. The 
intermediate states differ by the order of transfer of the 
two electrons and by the single-particle state of energy 

e p in the superconductor in which the virtual quasiparti- AR tunneling rate ( |Averin and Pekola] |2008[ ) 
cle with excitation energy a p = (A 2 + e 2 ) 1 / 2 is created. 
In addition to f2 p , the energy of the intermediate state 
includes the charging energy barrier Ei to the transfer of 
one electron from the normal electrode to the supercon- 
ductor. The standard BCS factors v p = \(l-e p /Q, p )/2] 1 / 2 
and u n = [(1 + e p /O p )/2] 1 / 2 enter Eq. (14) because v p is 



1995| . In the 
I give for the 



t.f2 A 2 r 

TaR= ^nJ def(e-E/2)f(-e-E/2) 



Y / a(±e-E l -E/2)\", (16) 



the amplitude of state p being empty in the BCS ground 
state, thus allowing the first electron transfer, while u p 
is the overlap of the doubly occupied orbital state p with 
the BCS ground state, which gives the amplitude of re- 
turn to the ground state after the second electron trans- 
fer. Since no trace of the intermediate states is left in 
the final state obtained after the whole AR process is 
complete, they should be summed over coherently, at the 
level of the amplitude A(ek,£i)- By contrast, the initial 
states efc, e; of the electrons in the transition are left 
empty in the final state, and can be used to distinguish 
between different transition processes. This means that 
they should be summed over incoherently, in the expres- 
sion for the tunneling rate At small temperatures 
ksT <C A, one can neglect thermal excitations in the 



where 

a(e) = (e 2 



A 2 )- 1 / 2 In 



A 2)l/2 



A-e-(e 2 -A 2 ) 1 /2 



Equation ( 16 1 is well defined if the relevant energies in the 
amplitude a(e) do not approach the edge of the supercon- 
ducting energy gap, e ~ A, which gives logarithmically 
divergent contribution to Tar- This singularity can be 
smeared by many mechanisms, e.g., the non- uniformity 
of the gap A or finite transmission probability of the bar- 
rier. In the single-electron tunneling regime, the main 
broadening mechanism should be the lifetime of the in- 
termediate charge state in the AR process, and can be 
accounted for by replacing in Eq. ( 14 ) the energy Ei with 
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FIG. 6 Real time detection of Andreev tunneling in an iso- 
lated SEB shown as red in the scanning electron micrograph 
of (a) and its schematic in panel (b). The electrometer, shown 
in blue is used for counting the single-electron and Andreev 
tunneling rates. Panel (c) shows the tunneling rate for AR 
as blue and red dots for forward and backward directions. 
The green lines are theoretical calculations where the non- 
uniformity of the t unnel barrier is tak en into account. Figure 
adapted from Ref. (Maisi et al. 2011 1. 



Ej, — «7/2, where 7 is the rate of sequential lowest-order 
tunneling out of the intermediate charge configuration. 

Experimentally, individual AR processes can be ob- 
served directly in time domain in the hybrid SEB ( Maisi 



et al. 2011 1. Such observation allows one to extract the 
rates of AR tunneling shown in Fig. [6] as a function of 
the normalized gate voltage n g , which determines the en- 
ergies Ei and E of the transition, together with the the- 
oretical fit based on Eq. (16). One can see that Eq. (16) 



describes very well the shape of the curves. It requires, 
however, considerably, roughly by a factor 15, smaller ef- 
fective number N of the transport channels. In practice, 
the magnitude of AR tunneling rates that is larger by 
roughly a factor of 10 than the theoretical expectation 
for a given tunnel conductance Gt, is the usual feature 



of the tunnel junctions (see also, e.g., (Greibe et al. 2011 



Pothier et al. 1994)), and in principle could be accounted 



by the variation of the barrier thickness over the junction 
area. Unfortunately, there is so far no experimental or 
theoretical evidence that the barrier non-uniformity is 
indeed the reason for the discrepancy between the mag- 
nitude of the lowest-order and AR tunneling. 

In the structures without superconducting electrodes, 
multi-step electron transitions, in contrast to the AR pro- 
cesses, involve electron transfers in different directions 
and/or across different tunnel junctions, since a transfer 
of the two electrons in the same junction and the same 
direction can not make the process energetically favor- 
able in the absence of the pairing gap A. In the simplest 
example of the NININ SET, the two-step co-tunneling 



process in the CB regime discussed qualitatively at the 
beginning of this Section, consists of two electron trans- 
fers across the two junctions of the transistor. Quantita- 
tively, the rate of this process is dominated by inelastic 
contribution Ti n , in which the single-particle states e p , e q 
of electrons in the central island of the transistor involved 
in the transfers are different ( |Averin and Odintsov 1989 1. 
As a result, the occupation factors of these states are 
changed, i.e. electron-hole excitations are created after 
the process is completed, the fact that implies that contri- 
butions to the tunneling rates from different e p , e q should 
be summed over incoherently. The elementary amplitude 
A of this process consists then only of a sum over the two 
possibilities, one in which an electron is first transferred 
onto the island increasing the charging energy of the in- 
termediate state by E^, and the other, in which an 
electron tunnels first from the island still increasing the 
charging energy but by a different amount E^: 



1 



A - Tk M €n + E (+) 



1 



£fc €i+E(-) 



-)• (17) 



The total rate Ti n is given then by the sum of \A\ 2 over all 
single-particle states involved with the appropriate equi- 
librium occupation factors, and can be expressed directly 
through the junction conductances as: 

r m (£) = ^ZZT f de k de p de q deif(e k )f(e q ) 



2-rre 4 

x [1 - f(e p )} [1 - f(e t )]6(E - e p + e k - e, + e q ) 



' e p + E(- 



e k ei + E(-) 



-I 2 , (18) 



where E — eV (see Fig. [5]) is the energy gain due to the 
transfer of electron charge e through both junctions of 



the transistor. Equation ( 18 ) shown of the second-order 



electron co-tunncling that involves one virtual intermedi- 
ate stage is indeed smaller than the rate ^ of sequential 
single-electron tunneling roughly by a factor RqGt <C 1. 
The derivation above also makes it clear that the rate 
of the multi-step electron transitions that go through n 
virtual intermediate stages with larger n would be sup- 
pressed much more strongly, by a factor (RqGt)™. 

If the energy gain E and thermal energy fc^T are 
smaller than the charging energy barriers £ (±) , Eq. (fl8|) 
for the inelastic co-tunneling rate can be simplified to 



hCiG-2 



1 



1 



127re 4 y E(- 



E(~) 



2 E[E* 



{2-Kk B Tf 



1 _ e -E/k B T 



This equation shows that, as a result of creation of exci- 
tations in the process of inelastic co-tunneling, its rate 
decreases rapidly with decreasing E and T. At very 
low energies, the process of co-tunneling in the NININ 
transistor will be dominated by the elastic contribution, 
in which an electron is added to and removed from the 
same single-particle state of the transistor island, with- 
out creating excitations on the island. Because of the 
restriction on the involved single-particle states, the rate 



11 



of such elastic contribution contains an additional factor 
on the order of SE/Ec ( |Averin and Nazarov 1990) and 
can win over T in only at very low temperatures, practi- 
cally negligible for the structures based on the /j,m-scale 
metallic islands considered in this review. 

The approach to multi-step electron transitions in 
the single-electron structures illustrated in this Section 
with the examples of Andreev reflection and electron co- 
tunneling can be directly extended to other higher-order 
tunneling process es, e.g. co-tunneling of a Cooper pair 
and an electron (Averin and Pekola 2008), which to- 



gether with Andreev reflection and electron co-tunneling 
limit in general the accuracy of control over sequential 
single-electron transitions. 



D. Coulomb blockade of Cooper pair tunneling 

In contrast to the tunneling processes considered 
above, which involve electrons in the normal-metal elec- 
trodes, tunneling of Cooper pairs in a junction between 
two s uperconductors is intrinsically a dissipationlcss pro- 
cess (Josephson 19621) . As such, it should not be char- 



acterized by a tunneling rate but a tunneling amplitude. 
Quantitative form of the corresponding term in the junc- 
tion Hamiltonian can be written most directly at low 
energies, fcgT , Eq <C A, when the quasiparticles can not 
be excited in the superconducting electrodes of the junc- 
tion, and the tunneling of Cooper pair is well-separated 
from the tunneling of individual electrons. In this regime, 
a superconductor can be thought of as a Bose-Einstein 
condensate of a "mesoscopically" large number of Cooper 
pairs which all occupy one quantum state. Transfer of 
one pair between two such condensates in the electrodes 
of a tunnel junction does not have any non-negligible 
effects on the condensates besides changing the charge 
Q = 2en on the junction capacitance by ±2e. There- 
fore, the part of the Hamiltonian describing the tunnel- 
ing of Cooper pairs should contain the terms accounting 
for the changes of the charge Q. Using the standard no- 
tation —Ej/2 for the amplitude of Cooper-pair tunnel- 
ing and including the charging energy p), one obtains 
the Hamiltonian of an SIS tunnel junction or, equiv 
alently, Cooper-pair box in the following form (Averin 



at al. 1985 Buttiker 1987) 



H = AE c {n-n g f 



Ej \ - 
2 4^ 



n){n± 1| 



(20) 



Here n is the number of Cooper pairs charging the to- 
tal junction capacitance, and n g is the continuous (e.g., 
gate-voltage-induced) charge on this capacitance nor- 
malized now to the Cooper pair charge 2e. Similarly 
to the sequential tunneling rates, the Cooper-pair tun- 
neling amplitude in the Hamiltonian ( 20 1 is a macro- 
scopic parameter which receives contributions from all 
Cooper pairs in the condensate, and can be expressed 
directly through the t unnel conductance Gt of the junc - 
tion, Ej = 7rGTA/2e (Ambegaokar and Baratoff 1963), 



in agreement with the simple fact that the amplitude of 
the two-electron tunneling should have the same depen- 
dence on the barrier transparency as the rate of tunneling 
of one electron. In the situation of the junction ( 20 1 real- 



ized with the actual Bose-Einstein condensates of atoms, 
such "Bose Josephson junction" can contain relatively 
small total number of particles, and the tunnel ampli- 
tude varies then with the difference n of the number of 
particles in the two condensates - see, e.g., (Averin et al. 



2008 Cheinet et al. 2008 Foiling et al. 



2007). 



Dependence on the ground state of the Hamiltonian 
(20) on the induced charge n g allows for qualitatively 



similar control of the individual Cooper pairs as indi- 
vidual electrons in the normal-state SEB discussed in 
Sec. |II.B| If Ej <C E c , precisely one Cooper pair is trans- 
ferred through the junction, changing n by ±1, when- 
ever n g passes adiabatically through a degeneracy point 
n g = 1/2 modulo an integer. This leads to the same 
staircase-like dependence n(n g ) as in the normal case, 
but with each step corresponding to the transfer of one 
more Cooper pair with the increase of n g by 1. The 
main new element of the superconducting situation is 
that the SIS junction is intrinsically a coherent quantum 
system without dissipation, and if extrinsic sources of de- 
coherence can be made sufficiently weak, should exhibit 
reversible dynamics of a simple quantum system. For in- 



stance, close to the degeneracy point n g 



1/2 the 



two charge states with the same electrostatic energy, n 
and n + 1, are coupled by coherent quantum mechani- 
cal tunneling of a Cooper pair, and the junction behaves 
as the very basic quantum two-state system ( Bouchiat 



et al. 1998 Nakamura et al. 19991. Such two-state dy- 



namics and general coherent quantum dynamics of the 
Hamiltonian ( 20 ) serve as the basis for the development 
of superconductin g quantum information devices - for 
reviews, see, e.g., (Averin 2000 



Makhlin et al. 2001) 



Superconducting junctions also exhibit the dynamics 
similar to the single-electron tunneling oscillations. If the 
induced charge n g grows in time at a constant rate, so 
that effectively a dc "displacement" current / = 2en g is 
injected into the junction, Cooper pairs are transferred 
Cooper pairs are transferred through it in a correlated 
manner, one by one, giving rise to the "Bloch" oscilla- 
tions ( Averin et al. |1985 1 of voltage across the junction, 
with frequency / related to the current I: 



I = 2ef. 



(21) 



The Hamiltonian ( 20 1 and its extensions to the multi- 



junction systems can be used to design time-dependent 
periodic dynamics with frequency / which would to 
transfer precisely one Cooper pair per period and there- 
fore produce the dc current quantized according to 
Eq. (21). Although the system dynamics employed for 



such Coo per-pair pumping ca n be of different kinds - 
see, e.g., (Hoehne e t al. 2012), the most typical is the 

in which the 



1991) 



adiabatic dynamics ( Geerligs et al. 
pumped charge i s rela t ed (|Aunola and Toppari[ 2003 
Mottonen et al.\ 120081 IPekola et al.\ 119991) to Berry's 
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PHOTON 
ABSORPTION 
and TUNNELING 




FIG. 7 A simple schematic showing photon absorption by a 
generic tunnel junction, and the inelastic electron tunneling 
from the left side of the barrier to right. 



phase or (Faoro et al. 2003) its non-abelian extensions 



E. Influence of environment on tunneling 

Tunneling in small junctions is influenced by electro- 
magnetic environment. The tunneling rates are modi- 
fied by photon absorption or emission; Figure [7] depicts 
schematically a process where tunneling rate in a generic 
junction is enhanced by absorption of a photon from the 
environment. 

The general theoretical framework of how this hap- 
pens was put forward in seminal works by Devoret et 



al. (Devoret et al. 1990) and Girvin et al. (Girvin et al. 



1990P , an d later expanded in the rev iew by Ingold and 
Nazarov (Ingold and Nazarov 1992 1. The golden rule 



type tunneling rates discussed in the earlier sections get 
modified as 



r= -jp- / / dEdE' v x {E - 5E)v 2 {E') x 



h{E-8E)[l- f 2 (E')]P{E-E'), 



(22) 



where Vi(E), i = 1,2 are the normalized densities of 
states (DOS) in the two electrodes, fi(E) are the cor- 
responding energy distributions in the electrodes, and 
SE is the energy cost in the tunneling event. The func- 
tion P(E) can be interpreted as the probability density 
to emit energy E to the environment, which becomes a 
delta-function in the special case of a junction with per- 
fect voltage bias. The P(E) can be calculated as the 
following transformation using the phase-phase correla- 
tion function J(t): 



P{E) = 



1 

2nh 



exp (J(i) + T Et)dt. 



(23) 



By modeling the environment by a frequency u>/2ir de- 
pendent impedance Z(w) in thermal equilibrium at tem- 
perature Tenv, one obtains 



J(t) = 



duj Re[Z(w) 



[coth( 



2k B T a 



)(cos tot — — i sin ut I ,(24) 



where Rk = h/e 2 is the resistance quantum. 

Often one can assume that the unintentional environ- 
ment can be modeled as a wide-band dissipative source 
in form of and RC circuit. For a purely resistive and 
capacitive environment 



Re[Z(oj)} = R/[l + (ojRC) 2 



(25) 



where R is the resistance of the environment and C is the 
total capacitance including the junction capacitance and 
parallel shunt capacitors. This rather simple model has 
been successfully applied to explain several experimental 
observations, see, e.g., (|Hergenrother et al. 



tinis and Nahum 19931. For a system with intentionally 



enhanced capacitance, it could be used to account for ex- 
perimental improvement of the characteristics of an NIS 



junction and of a single-electron turnstile (Pekola et al. 



2010), and further improvements have been obtained in 



Refs. (Saira et al. 2012a 2010), see Sec. 3.2. We show 



in Fig. [8jan NIS-junction and its current-voltage charac 
teristics under different experimental conditions. 

Focusing to the single-electron sources, the environ- 
ment has at least two effects to be considered, (i) The 
coupling of the black-body radiation of the hot surround- 
ing environment can induce photon-assisted tunneling, 
(ii) The intentionally fabricated on-chip environment in 
the immediate vicinity of the single-electron circuit serves 
as a filter against the external noise. Moreover, it can 
influence the tunneling rates in a way improving the per- 
formance (Bubanjal |2011[ ILotkhov et al\ |2001[ IZorin 



et al. 2000p . Detailed discussion on the error processes 



m pumps, including those due to coupling to the envi- 
ronment, is given further in this review. 



F. Heating of single-electron devices 

Single-electron circuits operate optimally at low tem- 
peratures. The standard condition is that ksT <C Ec, 
where Ec is the characteristic charging energy scale. An- 
other condition in superconductor based devices is that 
kgT <C A, where A is the energy gap of the supercon- 
ductor. Since thermal errors in synchronized transfer of 
electrons are typically proportional to e~ E ° R / kBT , where 
-E-CR is the characteristic energy (in the previous exam- 
ples Ec or A), it is obvious that temperature needs to 
be more than an order of magnitude below Ecn/kB- At 
low temperatures the overheating becomes a critical is- 
sue (Giazotto et al. 2006). The energy relaxation be- 



tween the electron system and the bath, typically formed 
by the phonons, becomes increasingly slow towards low 
temperatures. Moreover, the various heating rates are 
typically not scaling down similarly with decreasing tem- 
perature. 

Heat is injected to the electron system, first and fore- 
most, as Joule heating due to current in a biased circuit. 
Other sources of heat include application of dissipative 
gate voltages or magnetic flux injection, thermal radia- 
tion discussed in Section flLEI and shot-noise induced dis- 
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FIG. 8 NIS junctions influenced by hot environment, (a) Ge- 
ometry of a NIS junction made of aluminum (low contrast) 
as the superconductor and copper (high contrast) as the nor- 
mal metal. The tapered ends lead to large pads, (b) Typical 
IV characteristics, measured at 50 mK for a junction with 
Rt = 30 kfi. Linear leakage, i.e. non-vanishing sub-gap cur- 
rent due to coupling to the environment, can be observed. The 
cyan line is the corresponding theoretical line from the P(E) 
theory and RC environment with dissipation R at T cnv —4.2 
K. (c) Measured IV curves of an NIS junction with Rt = 761 
kf2 on a ground plane providing a large protecting capacitance 
against thermal fluctuations (solid symbols) and of a similar 
junction with Rt = 627 kQ without the ground plane (open 
symbols). Solid lines present the theoretical results for capac- 
itance C = 10 pF (red line) and C = 0.3 pF (green line). The 
resistance and the temperature of the environment are set to 
R = 2 Q and T cnv = 4.2 K, respectively. The inset shows 
IV curves based on the full P{E) calculation as functions of 
the shunt capacitance C. The red and green lines are repro- 
duced on this graph from the main figure. Figure adapted 
from Ref. (iPekola et aZ.||2010b. 



sipation by back-action from a charge or current detector. 
The steady state temperature of the electron system is 
determined by the balance between the input powers and 
the heat currents via different relaxation channels. The 
injected energy relaxes to phonons via electron-phonon 
relaxation, to the leads by heat transport through the 
tunnel junctions, and by radiation to other dissipative 
elements in the cold circuit. We will discuss these pro- 



cesses in more detail below. 

Joufe heating and cooling: In a biased circuit, the total 
Joule power is P = IV where V is the overall voltage and 
/ the current. This power can, however, be distributed 
very un-evenly in the different parts of the circuit: in 
an extreme example, some parts may cool down whereas 
the others are heavily overheated. Let us focus on dis- 
sipation in biased tunnel junctions. The basic example 
is a tunneling process in a junction between two conduc- 
tors with essentially constant density of states, which is 
the case presented by normal metals. At the finite bias 
voltage V the tunneling electron leaves behind a hole-like 
excitation and it creates an excited electron in the other 
electrode, i.e., both electrodes tend to heat up. Quantita- 
tively we can write the expression of the power deposited 
in the, say, right electrode as 



Pt, 



1 



e 2 R T 



dE E[f L (E - eV) - f R (E)} 



(26) 



Here, Jl,r refer to the energy distributions on the left 
(L) and right (R) side of the junction, respectively. 
P R = when fi, = fn, i.e., when the temperatures of 
the two sides are equal. By symmetry, or by direct cal- 
culation we can verify that the same amount of power is 
deposited to the left electrode in this situation. Thus the 
total power dissipation equals P = Pl + Pr = V 2 /Rt = 
IV, as it should. 

If one of the conductors is superconducting, the 
current-voltage characteristics are non-linear and the 
power deposited to each electrode is given by 



P 



1 



N,S 



e 2 R T 



dEE N/s u s (E)[f N (E - eV) - f s (E)] 



(27) 

Here E N = E — eV, and Eg = E, where N, S now refer 
to the normal and superconducting leads, respectively. 
The overall heating is again given by IV, but in this case, 
under bias conditions eV ~ A, Pn on the normal side can 

" 2006| ) 



become negative (NIS cooling, (Giazotto et 



and Ps on the superconductor side is always positive, 
i.e. it is heated up. 

Finally, if both sides are superconducting, the current- 
voltage characteristics are highly non-linear, but due to 
symmetry P L = P R = IV/2 = P/2. 

Other heating sources: Overheating of a single-electron 
circuit can be caused by various other sources. Ac gate- 
voltages or ac magnetic fluxes can induce dissipative cur- 
rents and heating due to dielectric losses, and single- 
electron electrometry or electrometry by a quantum point 
contact detector can cause effective heating due to the 
shot-noise back-action coupling to the single-electron cir- 
cuit, just to mention a few possibilities. 

Energy relaxation by conduction to leads: If a differ- 
ence between the electronic temperatures TV, and T R of 
the left and right leads exists, AT = T^—Tr, heat Pl^r. 
can flow electronically through the tunnel barrier. In the 
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case of a normal-normal junction, we have 



Pi 



1 



e 2 R T 



dEE[f L {E)-f R {E)] 



6e 2 R T K ' 



Tl). 



(28) 



where in the last step we assume that the junction is not 
biased. For a small temperature difference AT about the 
mean T = (Tl +Tr)/2 of the two temperatures, we may 
then write the thermal conductance Gth = Pl^r/AT 
of a NIN tunnel junction as 



Gth = 



ir 2 k 2 B T 
3e 2 i? T ' 



(29) 



which is the Wiedemann-Franz law for a conductor with 
resistance Rt- For either an NIS or SIS junction, heat 
conductance is exponentially small at low temperatures 
due to A. Another mechanism for the heat flow is the 
diffusion in the leads. It is discussed, in particular in 
superconducting ones, in Section |ilI.B.3| 

Electron-phonon relaxation: Electron phonon relax- 
ation is one of the dominant and in many systems one of 
the best understood relaxation mechanisms. For a nor- 
mal metal conductor with a uniform temperature T that 
differs from the bath phonon temperature Tq, one may 



write quite generally (Wellstood et al. 19941 



P c -p = £V(T 5 - T 5 ), 



(30) 



where S is a material con stant of order 10 9 
WK" 5 m- 3 dGiazotto et aLl |2006b, and V is the volume 



of the conductor. This equation holds amazingly well at 
sub-kelvin temperatures for various metals, irrespective 
of their dimensions. In single-electron devices, we typi- 
cally consider dissipation in a small Coulomb-blockaded 
region, whose volume and thus, according to Eq. (301, 
the coupling to the phonon bath is weak. Due to the 
small dimensions, one can typically assume a spatially 
uniform energy distribution on the conductor, moreover, 
assumption of overheating with a well defined electron 
temperature is also justified quite generally. 

In some cases these assumptions are not necessarily 
valid. An important exception is given by superconduc- 
tors where energy relaxation via phonon emission be- 
comes extremely weak due to the energy gap. At low 
temperatures the relaxation is limited by the emission of 
2 A phonons corresponding to the recombination of quasi- 
particles into Cooper pairs ( |Rothwarf and Taylor 19671. 
In the past few years, several experiments have measured 
the relaxation rate in this context, see, e.g., (Barends 



et al. 2008), and the corresponding energy release rate 



was measured recently in ( |Timofeev et al. 2009a). Ac- 
cording to the latter measurement the recombination- 
related heat flux is strongly suppressed from that given 
in Eq. (30), being about two orders of magnitude weaker 

0.3T C , 
At 



than in the normal state at the temperature T ■ 
where Tc is the critical temperature of aluminium. 



(a) (b) 



6V, 



r& 



6V, 



FIG. 9 (a) Radiative heat flow is caused by the photons which 
carry energy between resistors Ri and R2 at temperatures Ti 
and T2 respectively. The heat transport can be modelled by 
having voltage fluctuations SVi as shown in (b). Here we have 
assumed total transmission. The assumption can be relaxed 
by adding a non-zero impedance to the loop. 



even lower temperatures the heat current is further sup- 
pressed, eventually exponentially as oc e~ A / kE,T . Be- 
sides recombination, also the diffusive heat conduction 
is strongly suppressed in a superconductor at T <C Tc- 
This means that a superconductor is a poor material as 
a lead of a single-electron source, where non-equilibrium 
quasiparticles are injected at the rate /. The situation 
can be improved by inserting so-called quasiparticle traps 
into the circuit, discussed in section [III.B.3| Yet a fully 
superconducting Cooper-pair pump can be dissipation- 
less ideally. 

Heating and cooling by radiation: Coupling of a junc- 
tion to the electromagnetic environment is associated 
with heat exchange. Hot environment can induce pho- 



ton assisted tunneling as discussed in subsection II. E 
above. The basic concept of radiative heat transport in 
an electric circuit has been known since the experiments 



of Johnson and Nyquist (Johnson 1928 Nyquist 1928) 



more than eighty years ago. Electromagnetic radiation 
on a chip has recently turned out to be an important 



channel of heat transport at low temperatures (|Meschke 


et al. 


2006 


Schmidt et al. 


2004 Timofeev et al.\ 


2009b). 



If two resistors Ri and R2 at temperatures T\ and T 2 are 
connected directly to each other in a loop the heat ex- 
change between them can be modeled by Langevin type 
circuit analysis, as indicated in Fig. [9] by the voltage 
sources producing thermal noise. Assuming an idealized 
quantum limit, where the circuit transmits all frequen- 
cies up to the thermal cut-off at wth = ksTi/h, the net 
heat current between the two resistors is given by 



P 7 = 



R1R2 



(Ri 



7Tfcg . 2 _ rp2\ 

R 2 yuh (l 2h 



(31) 



This is an interesting limit which applies for circuits 
on a chip where the stray capacitances and inductances 
are small enough such that the circuit low-pass cut-off 
frequency exceeds w t h- Equation (31) has an impor- 
tant limit for the maximum coupling with R\ = R2 
and for small temperature differences \Ti — T 2 \, namely 
Gth = P 1 /{T l -T 2 ) = (i:k 2 B T)/(6h) = G Q , the so-called 
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quantum of thermal conductance (Pendry 19831. Equa- 



tion (31 1 can be applied also in the case where one of 



the resistors is replaced by a NIN tunnel junction with 
the corresponding resistance. Another important case is 
that of a hot resistor R at temperature T env , discussed in 
the previous subsection. In this limit, with RC cut-off as 
discussed in |ILE| one finds that the heat absorption rate 
by a resistor or a normal tunnel junction (at T <C T env ) 



is given approximately by (Pekola and Hekking 2007[ ) 



RtC 



(32) 



III. REALIZATIONS 

A. Normal metal devices 



Single-electron tunneling effects provide means to 
transport electrons controllably one by one. In this re- 
spect the obvious choices are metallic single-electron cir- 
cuits and semiconducting quantum dots. The metallic 
ones can be either in their normal state or superconduct- 
ing, or as hybrids of the two. The quantum dots, metal- 
lic hybrids and superconducting circuits will be discussed 
in later sections. The first single-electron source was a 
metallic (non-superconducting) turnstile with four tun- 



nel junctions and one active gate (Geerligs et al. 1990b ). 
The word "turnstile" refers to a device that is voltage 
biased at Vb between the external leads, but where the 
transport of electrons is impeded under idle conditions 
because of an energy gap. Under the active gate op- 
eration, electrons are transported synchronously one at 
a time. The finite voltage determines the direction of 
charge transport, at the expense that the device is also 
dissipative. We will discuss a more recent version of a 
turnstile in the subsection IIII.BI 

The most impressive results at the early days of single- 
electron sources were obtained by metallic multi-junction 
pumps, operating in a non-superconducting state. A pro- 
totype of them, featuring the main principle is the three- 
junction pump, with two islands and a gate to each of 
them, see Fig. 10 (a). This kind of a pump was suc- 
cessfully operated in 1991 by Pothier at al 



1991 Pothier et al. 1991 1992). Figure 10 (b) demon- 



(Pothier 



strates the stability diagram of a three-junction pump, 
which is essentially the same as that of the more com- 
mon double-island quantum-dot circuit. The two axes 
here are the two gate voltages, n g i,ro g 2, normalized by 
the voltage corresponding to charge displacement of one 
electron, i.e., n g i = C g iV g ^/e, where C g t is the gate ca- 
pacitance of island i. The stability diagram consists of 
lines separating different stable charge states on the is- 
lands, indicated by indices (ni,n 2 ) in the figure. The 
important property of this stability diagram is the ex- 
istence of the nodes where three different charge states 
become degenerate, being the three states with the low- 
est energy. The pump is operated around such a node, 
setting the working point at this node by applying dc 




b) V = 




FIG. 10 A three-junction pump. Schematics shown in (a), 
where the pump is biased by voltage V and with gate voltages 
Ui and U2- (b) The stability diagram of the three-junction 
pump on the plane of the gate voltages at zero bias voltage. 
For operation of the pump, see text. Figure from Ref. ( |Poth-| 
ierl 119911). 



voltages to the two gates. Let us focus on one such node, 
that at rigi = n g 2 = 1/3. Now, if the temporally varying 
gate voltages with frequency / added to these dc gate 
biases are such that the cyclic trajectory encircles the 
node at n g \ = n g 2 = 1/3 counterclockwise, one electron 
is transported through the pump from left to right. The 
simplest implementation of such a cyclic trajectory is a 
circle around the node, which is represented by two equal 
amplitude (in n g i) sinusoidal voltages applied to the two 
gates, phase-shifted by 90 degrees. Let us take point A 
as the starting point of the cycle. There the system is 
in the charge state (0,0). Upon crossing the first degen- 
eracy line, the new stable charge state is (1,0), meaning 
that if an electron has to tunnel from the left lead to 
island 1, while moving in this part of the stability dia- 
gram. Reversible pumping is achieved when / is so slow 
that the transition occurs right at the degeneracy line. If, 
however, the pumping frequency is too fast, the tunnel- 
ing is not occuring before meeting the next degeneracy 
line, and the pumping fails. Roughly speaking, the tun- 
neling process is stochastic, where the decay time of the 
Poisson process is determined by the junction resistance 
(see section II. B I, and if the pumping frequency becomes 
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Ref. (Averin et al. 1993 Jensen and Martinis 1992). 



FIG. 11 Predicted relative co-tunneling induced error ver- 
sus inverse temperature for the multi-junction pumps, with 
N = 4 (circles) and N = 5 (squares). The computer simu- 
lations (points) and the predictions of analytic results (lines) 
are shown. Parameters are Rt = 20R K , f = 4 • l Q- i /{R T C), 
and CV/e = —0.15. Figure adapted from Ref. ( jJensen and| 



Martinis 1992) 



comparable to the inverse decay time for tunneling, the 
desired event can be missed. In the successful cycle, on 
the contrary, the system next crosses the degeneracy line 
between charge states (1,0) and (0,1), and under the 
same conditions, the system transits to the new stable 
charge state by an electron tunneling from the left is- 
land to the right one. In the remaining part of the cycle, 
on crossing the last degeneracy line, an electron tunnels 
from the right island to the right lead completing the cy- 
cle where charge e (one electron) has been transported 
from the left lead to the right one. By cyclically repeat- 
ing this path at frequency /, an average current / = ef 
runs from right to left, and this current can be read, for 
instance, by a regular transimpedance amplifier. 

One of the advantages of the three junction pump over 
the early turnstiles is that the device can be operated, 
in principle, reversibly, since no external bias voltage 
is needed. It can pump even against moderate bias. 
Another difference between the pump and the turnstile 
above is that in the pump there are no unattended is- 
lands, on which the charge would be poorly controlled. 
Yet early realizations using fully normal metal conduc- 
tors in both the turnstiles and pumps suffered from other 
error sources which made these devices relatively inaccu- 
rate, on the level of 1%, even at low operation frequencies. 
A fundamental error source in this case is co-tunneling, 
discussed in Section [II. C| To circumvent this problem, a 
pump with a longer array of junctions is desirable: the 
error rate due to co-tunneling is effectively suppressed by 
increasing the number of junctions in the array. 

Theoretical analysis of co-tunneling in multi- 
junction pumps in form of N junctions in series 
with non-superconducting electrodes was performed in 



Thermal co-tunneling errors were analyzed with focus 
on N = 4 and N — 5. The conclusion of the analysis 
was that under realistic experimental conditions, N = 4 
pump fails to produce accuracy better than about 10~ 5 , 
insufficient for metrology, whereas N = 5 should be 
sufficiently good at low operation frequencies, as far 
as c o-tu nneling is concerned. This is demonstrated in 
where relative error of 10~ 8 was predicted for an 
5 pump at the operation frequency of / — 1.3 MHz, 



Fig. 

N 



for a pump with junctions having Rrp = 500 kO and 
C = 0.6 fF, if we assume that the working temperature 
is T = 50 mK. All these parameters are quite realistic. 
However, one notes already here that the frequency at 
which such a multi-junction pump can be operated is 



very low. In the subsequent experiments ( Martinis et al. 



1994) the error rate of about 0.5 ppm could be achieved, 
but still orders of magnitude above the prediction 
based on co-tunneling for their circuit parameters and 
experimental conditions. The authors conjectured this 
discrepancy to arise from photon-assisted tunneling, see 
Section 
e.g.. 



|II.E[ considered in this context theoretically, 
in Refs. ( |Martinis and Nahum[ |1993| |White and| 



Wagner 19931. After the experiment in (Martinis et al. 
1994), focus was turned into N — 7 pump where further 
improved results, 15 pbb, could be obtained at pumping 
frequencies of about 10 MHz (|Keller et q7T] |1996|>. This 



impressive result, depicted in Fig. 12 however of limited 
applicability in metrology of redefining ampere due 
to the small magnitude of the current, was proposed 
to present a capacitance standard based on electron 
counting (Keller et al. 19991. Further analysis and 



experiments on erro rs in multi-junction pumps can be 
found 



e-g-, m 



Refs. 



2003 Kautz et al. 



. (Covington et al. 


1999 


2000 


)• 



2000 Jehl et al. 



Another succesful line of metallic single-electron 
pumps relies on a smaller number of junctions with the 
help of dissipative on-chip environment used to suppress 



harmful cotunneling and photon assisted tunneling (Ca- 
marota et al. 2012 Lotkhov et al. 20011). These devices 



will be discussed further in Sections llII.H.2l and IIV.D.2I 

B. Hybrid superconducting-normal metal devices 

1. Operating principles 

The hybrid turnstile, originally proposed and demon- 
strated in Ref. (Pekola et al. 2008), is based on a 



single-electron transistor where the tunnel junctions are 
formed between a superconductor and a normal metal, 
see Fig. [13] on top left. In principle, it can be realized 



either in a SINIS or NISIN configuration (Averin and 



Pekola 2008 Kemppinen et al. 2009a). However, it has 



turned out for several reasons that the former one is the 
only potential choice out of the two for accurate synchro- 
nized electron transport purposes (Averin and Pekolal 



2008). One reason is that in the NISIN structure, tun- 
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FIG. 12 The seven junction pump. The plot on the left 
shows the schematic of the pump, with six islands, each with 
a gate. The electrons are pumped to/from the external is- 
land on the top, and the charge on the island is detected by a 
single-electron electrometer. The plot in the middle shows the 
voltage V p on the external island vs time when pumping ±e 
with a wait time of 4.5 s in between. The graph on the right 
shows the pumping error vs temperature of the measurement, 
demonstrating the 15 ppb accur acy at temperatures below 100 
mK. Figure adapted from Ref. (Keller et al.\ 19961. 



neling strongly heats the island due to Joule power and 
weak energy relaxation in the small superconducting is- 
land, whereas in the SINIS case the island is of normal 
metal, better thermalized to the bath, and under proper 



operation, it can be cooled, too (Kafanov et al., 2009 
The NISIN turnstile may also suffer from unpredictab 
le/2e periodicity issues. Furthermore, a detailed analy- 
sis of the higher-order tunneling processes shows that co- 
tunneling limits the fundamental accuracy of the NISIN 
turnstile, whereas unc ertainties below 10~ 8 are p redicted 
for the SINIS version (lAverin and Pekolal 120081). Hence 



we focus on the SINIS turnstile here. 

The stability diagram of a conventional single-electron 
transistor is composed of Coulomb diamonds on the gate 
voltage V g - drain-source voltage 14 plane, see Fig. 14 
Gate voltages V s are again written in dimcnsionlcss form, 
normalized by the voltage corresponding to charge dis- 
placement of one electron, n g . In this case the adjacent 
diamonds touch each other at a single point at Vb = 0, 
implying that the charge state is not locked for all gate 
voltage values. The operation of the SINIS-turnstile, on 
the contrary, is based on the combined effect of the two 
gaps: the superconducting BCS gap expands the stabil- 
ity regions of the charge states and the neighboring re- 
gions overlap. The principle of operation of the turnstile 
is illustrated in Fig. 



14 



When the gate charge n g (t) is 
alternated between two neighboring charge states, elec- 
trons are transported through the turnstile one by one. 
A non-zero voltage, which yields a preferred direction 
of tunneling, can be applied since the idle current is 
ideally zero in the range |eVt>| < 2 A at any constant 
gate charge value. If the gate signal is extended to 
span k + 1 charge states, one obtains current plateaus 
with k electrons pumped per cycle. However, the first 



FIG. 13 The hybrid NIS turnstile. Top left: A scanning elec- 
tron micrograph of a SINIS turnstile, which is a hybrid single- 
electron transistor with superconducting leads and a normal 
metal island. Top right and bottom: current of a turnstile 
under rf drive on the gate at different operation points with 
respect to the dc gate position and rf amplitude of the gate. 
Figur e adapted from Refs. (fKemppincn 2009, Pckola ~f al.\ 
20081. 



plateau around the symmetric (degeneracy of two neigh- 
bouring charge states) dc position of the gate is optimal 
for metrology. Note that if a non-zero bias voltage is ap- 
plied across a normal-state SET, a gate span between dif- 
ferent charge states always passes a region where none of 
the states is stable and where the current can freely flow 
through the device (white region in Fig. [l4ja)). Hence 
the normal-state SET cannot act as a turnstile even in 
principle, except for an experimentally unfeasible gate- 
sequence where n g jumps abruptly between its extreme 
values in the Coulomb blockaded parts of the stability 
diagram. 



Figure 13 presents data obtained by a basic tu rnstile 
operated under various conditions (Pekola et al. 2008 1 . 



Several wide current plateaus with increasing gate ampli- 
tude A g can be seen. The gate drive is expressed here as 
n g (t) = n g o+A g w(t) where n g o and A g are the gate offset 
and drive amplitude, respectively. The gate waveform of 
unit amplitude, is denoted by w(t). The optimal gate 
drive is symmetric with respect to the two charge states: 
therefore in later sections we assume that n g o = 1/2. In 
these first experiments, accuracy of synchronized charge 
transport as I = Nef, with N the integer index of a 
plateau, could be verified within about 1%. 

A rough estimate for the optimal bias voltage Vb 
is obtained by considering the dominant thermal er- 
rors ( |Pekola et al.\ |2008[ ) . The probability of an electron 
tunneling against bias, i.e., "in the wrong direction" is 
given approximately by ~ exp(— eV^/k^T). This error 
would lead to no net charge transferred during a pump- 
ing cycle, but it can be suppressed by increasing Vb. On 
the other hand, increasing Vb increases the probability 
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FIG. 14 Schematic picture of pumping (a) with a normal 
SET, (b) with a hybrid SET with E c = A, and (c) with a 
hybrid SET with Ec = 2A. The shaded areas are the stability 
regions of the charge states n = and n = 1. The edges 
of the normal SET stability regions are drawn to all figures 
with dashed black lines. The long coloured lines represent 
the transition thresholds from states n = and n = 1 by 
tunneling through the left (L) or the right (R) junction in 
the wanted forward (F, solid line) or unwanted backward (B, 
dashed line) direction. The thick black line corresponds to 
pumping with constant bias vol tage eVb /A = 1 an d a varying 
gate voltage. Figure from Ref. (Kemppinen 20091. 
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of transporting an extra electron in the forward direc- 
tion. The magnitude of this kind of an error can be 
estimated as ~ exp(— (2A — eVbV&B? 1 ), since there is an 
energy cost given by the voltage distance from the con- 
duction threshold at 2A/e. Combining these conditions, 
we obtain a trade-off eVb ~ A as the optimum bias volt- 
age, where the thermal error probability is approximately 
~ exp(— A/k^T). The combined thermal error probabil- 
ity is ~ 10~ 9 at realistic temperatures of about 100 mK 
and with the BCS gap of aluminum, A/fee ~ 2.5 K. The 
exact optimum of the bias close to the value given here 
depends on many other processes to be discussed below. 
Experimentally, however, the choice eVh = A is a good 
starting point. 

The optimal gate drive amplitude A g lies somewhere 
between the threshold amplitudes for forward and back- 
ward tunneling which are A g j t = A/AE C and A gi bt = 
3A/4_E C for the optimum bias voltage, respectively. The 
sub-gap leakage is maximized at the degeneracy point 
rigo = 1/2. In this respect, a square-wave signal is op- 
timal. On the other hand, passing the threshold for 
forward tunneling too quickly tends to heat the island, 
whereas a sine signal can also cool it. Hence, the optimal 
waveform is of some intermediate form. 

The SINIS turnstile presents a choice of a single- 
electron source which is easy to manufacture and operate, 
and whose characteristics can be analyzed theoretically 
into great detail. It promises high accuracy, as will be dis- 
cussed in the next section. Its operation in a parallel con- 
figuration is straightforward thanks to the simple element 
of a single turnstile, and therefore it can yield higher cur- 
rents than the other fixed-barrier single-electron sources 
presented in the previous subsection. Thus it can be 
considered as a very promising candidate in providing a 
realization of ampere. 



FIG. 15 (a) Stability diamonds for single-electron tunneling 
(blue solid lines) and Andreev tunneling (red dotted lines) 
for sample with Ec > A. (b) Stability diamonds for Ec < 
A. (c) First pumping plateau of the high Ec device as a 
function of the gate voltage amplitude A a . The solid symbols 
show pumped current with / = 10 MHz and three different 
bias voltages. Dotted lines are the simulated traces with the 
corresponding biasing, (d) The same data as in (c) but now 
for the low Ec device showing excess current due to Andreev 
tunneling. 



2. Higher-oder processes 

As for the fully normal metal pumps, the idealized pic- 
ture of electron transport based on single electron tunnel- 
ing is disturbed by simultaneous tunneling of several elec- 
trons. Owing to the gap in the quasiparticle excitation 
spectrum of a BCS superconductor, elastic co-tunneling 
takes place only when the bias voltage over the device 
exceeds 2A/e. The turnstile operation is achieved with 
voltages well below this threshold and hence co-tunneling 
is suppressed, in contrast to purely normal metal devices. 
As a general rule, any process that leaves behind an un- 
paired electron on a superconducting electrode incurs an 
energy penalty equal to A. 

For hybrid structures, the lowest order tunneling pro- 
cess where the energy cost of breaking a Cooper pair can 
be avoided is Andreev tunneling \ Andreev 1964), i.e., a 
complete Cooper-pair tunneling through a junction. An- 
dreev tunneling has been studied thoroughly with single 



NIS junctions ((Blonder et al. 


1982 


Eiles et al. 


1993 


Greibe et al. 


2011 


Hckking and Nazarov[ 


1994: 


Lafarge 


et al. 


1993; Maisi et al. 


2011 


Pothier et al. 


1994 


Rajau- 


ria et al. 


2008) as well as in so-called Cooper pair split- 



ters where the electrons of a Cooper pair tunnel to differ- 



ent normal metal regions ( Herrmann et al. 2010 Hofstet 
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ter et al.\ [2009{ |Wei and Chandrasekhar| |2010[ ). In the 



(a) 



case of a SINIS turnstile, Andreev tunneling manifests 
itself as two electrons being added to or removed from 
the island. Consecutively, increasing the charging energy 
of a device makes Andreev tunneling energetically unfa- 



vorable, suppressing it (Averin and Pekola 2008 Maisi 



et al.\ [2011h . The impact of Andreev tunneling on the 
accuracy of a turnstile has been directly observed on the 
pumped current (Aref et al. 2011). In Fig. [l5|a,b) 



sta- 
bility diamonds for single electron and Andreev tunneling 
are shown for a high Eq and low Eq device, respectively. 
The pumping plateau of the high Eq device, shown in 
panel (c), is free of Andreev tunneling whereas the low 
Eq sample exhibits it as seen in panel (d). 

For high-charging energy devices where Andreev tun- 
neling is supressed, the process limiting the accuracy of 
the SINIS turnstile is co-tunneling of a Cooper pair and 
a single electron (Averin and Pekola 2008). In this pro- 



cess, the island will be charged or discharged by a single 
electron while another electron effectively passes through 
the device. The net energy change is that of the corre- 
sponding single-electron process, plus the energy gained 
in transporting the Cooper pair from one electrode to an- 
other, which equals 2eVb in the forward direction. Hence, 
the process cannot be made energetically unfavorable in 
a working turnstile. However, it can be supressed rela- 
tive to the first-order processes by making the junctions 
opaque enough. Ideally, to obtain an accuracy of 10~ 7 , 
one needs to limit the speed of an aluminum based turn- 
stile to a few tens of pA (Averin and Pekola 2008). This 



theoretically predicted maximum operation speed is ex- 
pected to slow down by an additional factor of 3 due to 



nonuniformity of the tunnel barriers (Aref et al. 2011 



Maisi et al. 2011 ). Thus 10 pA is expected to be the op- 
timum yield per aluminum based turnstile. In addition to 
the Cooper-pair electron co-tunneling, the co-tunneling 
of two Cooper pairs through the device increases the 
leakage current (Zaikin 1994). In optimized devices dis- 



cussed above, the Cooper-pair electron co-tunneling is 
nevertheless the dominant one limiting the accuracy since 
its threshold is exceeded in the turnstile operation and it 
is of lower order than the Cooper pair co-tunneling. 



3. Quasiparticle thermalization 

Single-electron tunneling to or from a superconductor 
will generate quasiparticle excitations. Once created, the 
excitations carry an energy of A, which enables them to 
cross the tunnel barrier to the normal metal if the electro- 
static energy cost is lower than A. Hence, they constitute 
a potential source of pumping errors for the hybrid turn- 
stile. Typically the excitations are injected close to the 
gap edge. Also, the quasiparticles relax quickly internally 
compared to the weak recombination rate, so that at low 
temperatures we can assume them to lie close to the gap 
edges and have a temperature T qp which is higher than 
the phonon bath temperature of the system. With this 
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FIG. 16 Two typical geometries for a superconducting bias 
lead: (a) A lead having a constant cross-section determined by 
the thickness d and width w. The length of the line is I. (b) A 
sector shaped lead characterized by an opening angle 9, initial 
radius ro, and final radius r t . For the picture 9 is set to 180°. 
The brown parts on top denote a quasiparticle trap connected 
via an oxide barrier. Panel (c) shows quasiparticle density n qp 
along a constant cross section line with various oxide trap 
transparencies k, and panel (d) along an opening line. In the 
plots, n qp is scaled by rii — DlPi n j /Ai , where the injenction 
area Ai equals wd for panel (c) and 9rod for panel (d). For the 
leads in (b) and (d), we also use notation x — (r — ro)/l with 
I — rt — ro, and have used values ro = 20 nm and rt — 5 /im. 



assumption, we can calculate the density of the quasipar- 
ticle excitations to be 



/>oo 

2D(E F ) / dEv s {E)e~ fiE 
V2ttD(E f )A^=, 



(33) 



where = l/(fc B T qp ), 



and we have assumed e 



/3A <<; 1 

19721. The 



and a negligible branch imbalance (Clarke 
tunneling rate caused by the excitations can be calculated 
from the orthodox theory expressions (see Sec. II. B I. It 
depends linearly on the density, and is independent of 



the biasing at low energies: T, 



•(iii 



l[2e 2 R T D{E F )} 



It should be compared to the rate at which we pump 
electrons. As discussed above, we can obtain roughly 
10 pA from a turnstile free of higher order tunneling 
errors at accuracy of 10~ 7 . The tunneling resistance 
of such a device is approximately Rt — 1 MSI. To 
ensure that the quasiparticle excitations do not cause 
errors on this level, we require the tunneling rate to 
satisfy T qp < 10 -7 • 10 pA/e. With parameter values 
D(E F ) = 1.45 • 10 47 J^m" 3 and A = 200 ^eV, we need 
n qp < 0.04 /um -3 . Such a level is demonstrated in an 
experiment without active driving of the system ( |Saira| 
et al. 2012a[ ), and is sensitive to the filtering and shield- 
ing of the sample. Also, the trapping of quasiparticles 
was shown to be important in this experiment. 

Next, we consider the relaxation of the quasiparti- 
cles. In turnstile operation, injection of hot quasipar- 
ticles through the tunnel junction drives the quasipar- 
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ticle system of the superconductor actively out of equi- 
librium. We model the quasiparticle relaxation in the 
superconductor in terms of heat flow and obtain a diffu- 
sion equation for n qp . Such an approach has been used to 



model several experiments (|Knowles et al. 


2012 


O'Neil 


et al\ 2011 


3 eltonen et al.\ |2011| |Rajauria et al. 


|2009l 


Ullom et al. 


1998 


|. The heat flow of quasiparticles J 



follows the equation V • J = —p, where p is the power 
per unit volume removed from the quasiparticles. We use 
Fourier's law of heat conduction: J = — KsVT qp , where 
(/3A) 2 e~^ A is the heat conductivity of a 



6 L T„- 



superconductor ( |Bardeen et al. 1959 ) . Here L is Lorenz 
number and p n the resistivity in normal state. By taking 
the derivatives only over strong exponential dependencies 
and using Eq. (33), we obtain a diffusion equation 



DV 1 



P- 



(34) 



where the coefficient D = 



constant. 



- , , , . , is assumed to be 
To write down the source term on the right side 
of Eq. ( 34 ) , we consider the available mechanisms of heat 
conduction. Electron-phonon coupling is an inherent re- 
laxation mechanism for quasiparticles inside a supercon- 
ductor. However, it is so weak that the resulting decay 
length of n qp is t ypically on the millimeter scale (Mar- 



tinis et aL| |2009| |Peltonen et aZ.[ |2011| ). Typically, to 
enhance the relaxation, one uses so-called quasiparticle 



traps ( 


O'Neil et al. |2011 


Pekola et al.\ 


2000 


Rajau- 


ria et al. 


2009 


1, which are normal metallic regions con- 



nected to the superconductor either directly or via an ox- 
ide layer. Once the hot quasiparticles enter the trap, the 
stronger electron-phonon relaxation in a normal metal 
removes their excess energy. A perfect quasiparticle trap 
forces the quasiparticle temperature at the interface to 
equal the electronic temperature of the normal metal. In 
the context of Eq. (34), this can be implemented as a 
boundary condition for n qp . The boundary condition at 
the junction is obtained by setting the heat flow equal to 
the power injected by quasiparticle current. 

In the case the trap is connected via an oxide bar- 
rier, the heat is carried by quasiparticle tunneling. The 
orthodox theory result for the source term in such a con- 
figuration is 



2a T 

17 J 1 

e 2 D(E F )d 



dEEn s {E){e- pE -e- p ° E ) 



(35) 



which is obtained by setting the chemical potential dif- 
ference of the trap and the superconductor to zero and 
assuming the diffusion to take place in two dimensions 
which is well justified for the thin films typically used in 
the samples. We also assumed fcsT qp <C A. Here, <jt 
is the electrical conductance per unit area of the trap, 
d thickness of the superconducting film, f3 — l/(fcBT qp ), 
and Pq = 1/ (fog To), where Tq is the temperature of the 
normal metal electrons. We denote by n qp o the quasipar- 



ticle density of a fully thermalized superconductor, i.e., 
one where T qp = To. 

Let us consider some typical geometries of supercon- 
ducting leads used in devices. First, take a lead with a 
constant cross-section, as shown in Fig. |Io|a). We as- 
sume that a heat flow Pi n j is injected at one end of the 
line, and that the other end is thermally anchored by a di- 
rect trap. For the lead itself, we assume a trap connected 
via an oxide barrier to be located on top. We can solve 
Eqs. (34H35) analytically in one dimension to obtain 



It 



qp{x) '- 

Here, k 



l 

DVk y 



-I Pi, 



y/2(k B T np A)V 2 d 

the wire starting at the injection side (x 



l)-^+«9P0- 

and x is the coordinate along 
0) and end- 



16] (c), we show 
ues of gt- The 



ing at the direct trap (x = I). In Fig. 
the quasiparticle density for various vab 
lowest gt corresponds to the case where the quasiparti- 
cles only diffuse through the wire and then relax at the 
direct contact. At higher transparencies, the oxide trap 
starts to help for the relaxation as well. If we use pa- 
rameter values, d = 50 nm, w = 100 nm, I = 1 pm and 
T qp = 130 mK which are typical for fabricated samples, 
we see that a typical injection power of P{ n j = 2 fW will 



yield n q 



''qpO 



10 /im without the oxide trap and 



even with the highest transparency <tt = (100 flpm 2 ) 1 
which is possible to fabricate without pinholes ( Brenning 



et al. 2004), we get only an order of magnitude improve- 
ment. 

To decrease the quasiparticle density to the acceptable 
level discussed above, one needs to optimize the lead 
geometry as well. Therefore, let us consider a lead 
that widens as shown in Fig. 16 ]b). In this case, 
we can solve a one-dimensional diffusion equation in 
polar coordinates. The junction is assumed to be 
located at radius r = tq, and the direct contact trap 
to begin at radius r = r t . Thickness of the lead 
and the overlaid trap are as in the previous example. 
The solution of Eqs. (34)-(35) can be expressed with 
modified Bessel functions I a and K a as n qp (r) = n qp o + 
l Pi, 

DVk 8r 

(h(Vkr ) 
Fig. pTd) 



* MM^^r'loiVkr)}. In 
we show n qp (r) for various transparen- 
cies of the oxide trap. The lowest transparencies, again, 
correspond to pure diffusion limit. Note that the quasi- 
particle density at the junction depends only weakly on 
the transparency of the trap: Due to the logarithmic 
dependence, changing the transparency by several orders 
of magnitude makes less than an order of magnitude 
difference to n qp (ro). In a widening lead, heat sinking is 
made efficient by spreading the heat to a larger volume, 
and the area of the trap contact is also increased. By 
using realistic parameter values, d = 50 nm, 8 = n/2, 
ro = 50 nm, r t = 5 pm, T qp = 130 mK, p n = 10 nf2m 
and Pi n j — 2 fW, we see that it is possible to reach 
n qp < 1 /im~ 3 at the junction even without an oxide 
trap. By increasing the thickness of the electrode by a 
factor of ten would then start to be sufficient for the 
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metrological accuracy requirements. 



Several experiments ( 


Knowles et al. 2012 O'Neil 


et al. 


2011 


Rajauria et al. 




2009; 


Ullom et al. 




1998) 



show that the above diffusion model is valid for quasipar- 
ticle densities of the order of n qp ~ 10 /zm~ 3 . A smaller 
quasiparticle density required for metrological applica- 
tions implies that the absolute number of quasiparticlcs 
in the conductors becomes very small. With a typical 
volume of a lead, 100 nm x 100 /J.m 2 , the quasiparticle 
number is N < 1 with n qp < 0.1 /im~ 3 . It is not cur- 
rently obvious if such a situation can be treated with the 
diffusion model or whether a more elaborated theory is 
required. Pumping experiments on metrological accuracy 
can provide a way to shed light on such a situation. 



C. Quantum-dot-based single-electron pumps and turnstiles 

In this section, we introduce semiconducting quantum 
dots and review their applications as single-electron cur- 
rent sources. 



1. Indroduction to quantum dots 

In contrast to conventional three-dimensional bulk 
conductors or two-dimensional more exotic conductors 
such as quantum Hall systems or graphene, semiconduct- 
ing quantum dots can be regarded as zero-dimensional 
conductors, for which the electrons are tightly confined 
in all three spatial dimensions. Thus quantum dots show 
truly discrete excitation spectra that are reminiscent to 
those of natural atoms. One of the early key experiments 



on these artificial atoms ( Kastner 1993 ) was the obser 



vation of discrete quantum levels (Johnson et al. 1992 



Reed et al. 1988 Su et al. 1992 1 and the shell structure 



in the filling of the electron states (Tarucha et al. 1996 1 . 

The conceptual difference between small metallic is- 
lands studied in the previous sections and quantum dots 
is that the Fermi level and hence the conduction elec- 
tron density in the metallic islands is high, making the 
energy spacing between the spatially excited electron 
states extremely small. The metallic system can be typ- 
ically described by a constant density of states as op- 
posed to the strongly peaked density of states in quan- 
tum dots. Furthermore, quantum dots can contain a 
low number of electrons in the conduction band ranging 
from zero (|Ashoori et al. 1993[ Elzerman et al. 



2003 



Lim et al. |2009 1 to more than hundreds similar to nat- 
ural atoms, whereas the corresponding number is orders 
of magnitude higher for metallic systems. In fact, the 
sharp potential created by a single donor atom in silicon 
can also be considered to be an ultra-small quantum dot. 
By connecting such natural atoms to electron reservoirs, 
for example, single-electron transistors (SETs) have been 



fabricated (Lansbergen et al. 2008 Tan et al. 2010) 



Figure |17| shows different types of quantum dot archi- 
tectures. The most conventional quantum dots are based 
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FIG. 17 (a) Lateral and (b) vertical quantum dot arrange- 
ments. All quantum dot pumps and turnstiles discussed fur- 
ther in this paper are in the lateral arrangement. The elec- 
trons tunnel between the dot and the source and drain reser- 
voirs. The tunnel barriers between the dot and the reservoirs 
are created either by (a) the electrostatic potentials of nearby 
gate electrodes or (b) by materials choise such as AlGaAs. 
The gate arrangement for (c) the accumulation and (d) de- 
pletion mode quantum dots in the lateral arrangement. 



on two-dimensional degenerate electron gas (2DEG) that 
forms either naturally to the interface between AlGaAs 
and GaAs (Chang et al. 1974) or that is induced to the 



interface between silicon and silicon oxide by an external 



be fabricated from epitaxially grown nanowires 


Fasth 


et al. 


2007 


Nadj-Perge et al. 


2010 


Ohlsson et al. 


2002). 



In the conventional dots, the confinement is very strong 
in the direction perpendicular to the interface. Etching 
techniques, local anodic oxidation (Held et al. 19971, or 
metallic electrodes [see Fig. 17 'c) and (d)J can be em- 
ployed to provide the electrostatic potential defining the 
well for the electrons in the plane of the interface. The 
in-plane diameter of this type of dot can vary from tens 
of nanometers to several micrometers. Thus there are 
plenty of atoms and electrons in the region of the dot 
but most of them lie in the valence band and require 
an energy of the order of one electronvolt to be excited. 
Since the relevant energy scales for the spatial excita- 
tions and the single-electron charging effects are orders 
of magnitude lower, the occupation of the valence states 
can be taken as fixed. 



In the effective mass approximation (Ando et al. 



1982), the details of the electrostatic potential and the 
effects of the valence electrons in the solid are coarse 
grained such that only the electrons in the conduction 
band are taken into account, and these electrons are 
treated as particles in the smooth potential defining the 
dot. This description has proved to reproduce several 
important experimental findings both qualitatively and 
quantitatively (Ando et al. 1982), and it provides insight 



into the single-electron phenomena in quantum dots. Es- 
pecially, the potential barriers arising from the gates 
defining the dot can be visualized just for the small num- 
ber of electrons in the conduction band. 
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FIG. 18 The first single-electron current source based on 
quantum dots by Kouwenhoven et al. ( 1991a|b| . (a) SEM 
image of the device from the top, (b),(c) operation principle, 
and (d) measured TV curves reported. The gate configuration 
corresponds to the case in Fig.flTjd). The different IV curves 
are measured while driving the turnstile with different center 
gate [gate C in panel (a)] voltages, rf amplitudes, and phase 
differences. The curves are not offset and the dashed lines 
show the current levels nef with n = — 5, . . . , 5. 



2. Pioneering experiments 

Although quantum dots hold a much smaller number 
of electrons than metallic islands, probably their great- 
est benefit is that the tunnel barriers can be formed 
by electrostatic potentials and controlled externally by 
gate voltages. Thus the height of the potential barrier, 
through which the electrons tunnel to the source and 
drain reservoirs, can be controlled in situ. This property 
provides fruitful grounds for electron pumping since the 
dependence of the tunneling rate on the barrier height 
and hence on the voltage of the gate electrode is typi- 
cally exponential. 

The first experiments employing quantum dots for 
fre quency-locked single-electron transport were reported 
by Kouwenhoven et al. ( 1991a|b l [see also Kouwenhoven 
( 1992p j. Here, they used surface-gated GaAs dots as 
shown in Fig. [l8^a). The negative voltages on Gates 
C, F, 1, and 2 deplete the 2DEG that is located 100 nm 
below the surface, thus defining the quantum dot in the 
center with a radius of about 300 nm and charging en- 
ergy 2E C = e 2 /2C s = 0.34 meV. (Gates 3 and 4 are 
grounded and do not deplete the 2DEG.) 

In addition to dc voltages defining the dot, 180-degree- 
phase-shifted sinusoidal rf drive is superimposed to Gates 
1 and 2, lowering one barrier at a time. This rf drive in- 
duces a turnstile operation as shown in Figs. |18[b) and 
18 'c) for negative bias voltage on the left side of the dot: 
when the voltage at Gate 1 is high (low tunnel barrier) 
and low at Gate 2 (high tunnel barrier), an excess elec- 
tron enters the dot through the left barrier [Fig. [l8^b)], 
and when the voltage at Gate 1 is low and high at Gate 2, 
the electron escapes through the right barrier [Fig.|18[c)]. 
Thus the average dc current through the device in the 
ideal case is given by I p — ef, where / is the operation 
frequency. For bias voltages greater than the charging 
energy, \eV\ > Eq, more than a single electron can be 
transported in a cycle yielding ideally I p — nef, where n 
is an integer. Signatures of this type of current quanti- 
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FIG. 19 (a) Schematic illustration of the device and (b) 
observed current plateaus during the turnstile operation by 
Nagamune et al. (1994). 



zation were observed in the experiments (Kouwenhoven 
|1992| Kouwenhoven et al] 1991a|b ) and are illustrated in 
Fig. |18[d). The current through the device as a function 
of the bias voltage tends clearly to form a staircase-like 
pattern with the step height ef. This was the first exper- 
imental demonstration of current quantization in quan- 
tum dot structures. Note that in addition to turnstile 
operation, Fig. 18 'd) also shows pumping of electrons 
against the bias voltage for certain phase differences of 
the driving signals. The error in the pumped current is a 
few percent, falling somewhat behind of the first exper- 
iments on metallic structures reported by |Geerligs et al. 
(|X990bb. 



The second set of experiments on single-electron turn- 
stiles based on quantum dots were published by |Naga- 



mune et al. (1994). Here, the quantum dot forms into a 
gallium arsenide 2DEG that is wet etched into a shape 
of 460- nm- wide wire as illustrated in Fig. [l9^a). Two 
230-nm-wide metallic gates are deposited perpendicular 
to the wire with a distance of 330 nm. This different bar- 
rier gate configuration and the higher charging energy of 
2Ec — 1.7 meV resulted in a clear improvement of the 
staircase structure as shown in Fig. 19 b) . However, the 
authors report that a parallel channel forms due to the 
rf operation and the effect of this channel is substracted 
from Fig. 19 'b). The authors estimate the accuracy of 
their device to be about 0.4% if the correction from the 
parallel channel is taken into account. 

In 1997-2001, a series of experiments was carried 
on so-called multiple-tunnel junction devices as electron 



et al. 



(Altebaeumer and Ahmed 




2001 




Altebaeumer 


2001 


Tsukagoshi et al. 


1998 


1997). Here, the 



phorus doped silicon that was etched such that a central 
region is connected to source and drain reservoirs by nar- 
row strips as shown in Fig(20ja) . The side gates near the 
strips are set to a constant potential and an rf drive on the 
central side gate induces a current that depends linearly 
on frequency as shown in Fig. 20 b). The explanation of 
this type of operation is that the dopants and disorder 
in the strips function as Coulomb blockade devices them 
selves rather than as single tunnel junctions, which gives 
rise to the term multiple-tunnel junction. Since these 
experiments were more motivated by applications in in- 
formation processing with only few electrons rather than 
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FIG. 20 (a) SEM image of the device and (b) pumped current 
through it in the experiments by Altcbacu mer and Ahmed] 
( |2001[ ). Different values of the current correspond to different 
dc voltages Vgi [see panel (a)]. 




20 -10 10 20 30 

v D (mV) 



FIG. 21 (a) SEM image of the device and (b) observed cur- 
rent plateaus up to I MHz pumping frequency by |Ono and] 
|Takah ashi (2003) on a silicon quantum dot. 



finding a metrological current source, the accuracy of the 
device was not studied in detail. 



3. Experiments on silicon quantum dots 

The first step toward single-electron pumping in silicon 
was taken by Fujiwara and Takahashi (2001 ) as they pre- 



sented an ultra small charge-coupled device (CCD) and 
demonstrated that it could be used to trap and move indi- 
vidual holes controllably at the temperature of 25 K. This 
device was fabricated with silicon-on-insulator techniques 
( |Takahashi et al.\ [1995) and had two adjacent polysili- 
con gates acting as metal-oxide-semiconductor field ef- 
fect transistors (MOSFETs). Subsequently, a rather sim- 
ilar device with charging energy 2Ec — 30 meV shown 
in Fig. [2l]( a) was uti lized for electron pumping by |Ono 
and Takahashi (2003) at the temperature of 25 K. They 
obtained an accuracy of the order of 10 -2 up to 1 MHz 
pumping that was also the limitation set by the calibra- 
tion of their measurement equipment. Here, the electron 
pumping was based on two sinusoidal driving signals that 
are offset less than 180 degrees, which renders the chem- 
ical potential of the dot to move during the cycle. 



In addition to pumping, Ono et al. ( 2003 1 utilized the 



device shown in Fig. 21 &) as a single-electron turnstile. 



The operational principle is the same as in the pioneer- 
ing experiments wit h GaAs quantum dots described in 
Fig. [I8jb) and (c). |Ono et al.\ p003j) observed current 



steps of ef up to / = 1 MHz operation frequencies [see 
Fig. J22^a)]. The flatness of the plateaus was of the order 
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FIG. 22 (a) Measured current plateaus for different fre- 
quencies of the turnstile operation with the device shown in 
Fig. [2TJ a) by |Ono et al] |2003"] ) . (b) SEM image of the silicon 
quantum dot device and a schematic measurement setup em- 
ployed in the experiments by Chan et al.| ( |201l] ). (c) Measured 
current plateaus (solid line) and the corresponding theoretical 
curve (dashed line) by Chan et al. (20111. The insets show 
zooms at the n = (bottom) and n = 1 (top) plateaus. The 
red dashed lines show ±10~ 3 relative deviation from the ideal 
ef level. 



lie aluminum gates to define a silicon quantum dot in 
the electron accumulation layer of the device as shown in 
Fig. 22 ^b) . Although the relative variation of the current 
at the plateau they measured was below 10~ 3 for a broad 
range of source-drain voltages [see Fig.[22|c)], they could 
not strictly claim better than 2% relative accuracy in the 
current due to the inaccurate calibration of the gain of 
the transimpedance amplifier employed. These experi- 
ments were carried out at 300 mK phonon temperature 
bu t the sequential t unneling model used to fit the data 
by Chan et al. (2011) suggested that the electron tem- 



perature of the dot rose up to 1.5 K. It is to be studied 
whether 1.5 K was due to power dissipated at the sur- 
face mount resistors in the vicinity of the sample or due 
to the direct heating of the 2DEG from the electrostatic 
coupling to the driven gate potentials. 

Fujiwara et al. (2008) introduced a single-electron 



10 2 measured at 25 K. Chan et al. (2011 ) used metal- Ref. (Nishiguchi et al. 2006 ) I . In contrast to the pioneer- 



ratchet based on a silicon nanowire quantum dot with two 
polysilicon gates working as MOSFETs [see also Ref. ( Fu- 
jiwara et al. 2004)]. In general, ratchets generate direc- 
tional flow from a non-directional drive due to the asym- 
metry of the device. Here, an oscillating voltage is ap- 
plied to one of the gates such that an electron captured 
through it near the maximum voltage, i.e., minimum bar- 
rier height, and ejected through the other barrier near 
the minimum voltage. In fact, the number of electrons 
pumped per cycle depends on the applied dc voltages and 
current plateaus up to 5e/ were reported. Furthermore, 
nanoampere pumped current was observed at the 3e/ 
plateau with the pumping frequency / = 2.3 GHz. The 
error in the current was estimated to be of the order of 
10~ 2 for the experiment carried out at 20 K temperature. 

The first error counting experime nt s in silicon were 
carried out by Yamahata et al. (2011) [see also 
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ing error counting experiments by Keller et al. (19961, 



only a single silicon nanowire quantum dot was used as 
the current source and the electrons were steered into 
and out of a quantum dot coupled to a charge sensor. 
By opening the MOSFET separating the node from the 
drain reservoir, it was possible to use the same device as 
a dc current source. The observed pumping error was 
of the order of 10 -2 and was reported to be dominated 
by thermal errors at the 17 K temperature of the exper- 
iments. 



Very recently, Jehl et al. (2012) reported on frequency- 



locked single-electron pumping with a small quantum 
dot formed in metallic NiSi nanowire interrupted by two 
MOSFETs controlled by barrier gates. With rf drives on 
the barrier gates, they were able to pump currents be- 
yond 1 nA but the accuracy of the pump was not stud- 
ied in detail. The MOSFET channels in this device had 
very sharp turn-on characteristics requiring only about 
4.2 mV of gate voltage to change the conductivity of the 
channel by a decade, which can be important in reduc- 
ing unwanted effects from the gate voltage drive such as 
heating. 
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FIG. 23 Current plateau in the electron pumping experi- 
ments by Blumenthal et al. ( 2007 1 as a function of the middle 
gate voltage at 547 MHz operation frequency. The dashed 
lines show a = ±10 f A uncertainty in the electrometer cal- 
ibration. The top left inset shows the device used as the 
electron pump. The top right inset shows current plateaus 
at 1 GHz pumping frequency and the bottom inset shows the 
pumped current as a function of the operation frequency. 



4. Experiments on gallium arsenide quantum dots 



After the pioneering experiments discussed in 



Sec. III.C.2 the focus in single-electron sources based on 
gallium arsenide moved toward the idea of using surface 
acoustic waves (SAWs) to drive the single electrons in 
a one-dimensional channel — a topic to be discussed in 



Sec. III.D In this section, we focus on gate-controlled 
GaAs pumps. 



The seminal work by Blumenthal et al. (2007) took 
gate-controlled GaAs quantum dots a leap closer to 
a metrological current source, namely, they reported 
547 MHz (87.64 pA) single-electron pumping with one- 
standard-deviation (la) relative uncertainty of 10 -4 (see 
Fig. [23| . However, the authors did not report the full 
dependence of the pumping errors as functions of all con- 
trol parameters. As the device, they employed a chemi- 
cally etched AlGaAs-GaAs wire with overlapping metal- 
lic gates as shown in Fig. [23] Only the three leftmost 
gates L, M, and R were used such that 180-degree-phase- 
shifted sinusoidal driving signals were applied to the gates 
L and R in addition to dc voltages applied to all of the 
three gates. The amplitudes of the rf signals were cho- 
sen asymmetric such that the device can work as a pump 
rather than a turnstile. The charging energy of the device 
was estimated to be 2Ec = 1 meV and the experiments 
were carried out at the bath temperature of 300 mK. 

With a similar device architecture as shown in Fig. [23 
but using o nly two gates instead of three, |Kaestner 
et al. ( 2008a ) demonstrated that frequency- locked single- 



electron pumping can be carried out with a single sinu- 
soidal driving voltage, thus decreasing the complexity of 
the scheme. This type of single-parameter pumping with 
two gates is employed in the rest of the works discussed in 



this section. Maire et al. (20081 studied the current noise 
of a similar single-parameter pump at / = 400 MHz and 
estimated based on the noise level that the relative pump- 
ing error was below 4 %. Kaestner et al. (2008b) studied 



the robustness of the current plateaus as functions of all 
control parameters of the pump except the source-drain 
bias. They showed that single-parameter pumping is ro- 
bust in the sense that wide current plateaus appear in 
the parameter space but their measurement uncertainty 
was limited to about 10 -2 , and hence a detailed study 
of the behavior of the accuracy as a function of these 
parameters was not available. 



Wright et al. ( 2008 1 made an important empirical ob- 



servation that the accuracy of the single-parameter pump 
can be improved by an application of perpendicular-to- 
plane magnetic field [see also Ref. ( Wright et al. |2009 ) 
and Fig. 24jb)]. They applied fields up to 2.5 T and 
demonstrated that the n = 1 plateau as a function of 
the dc voltage on the non-driven gate widens noticeably 
with increasing magnetic field. In the further studies 



by Kaestner et al. (2009) and Leicht et al. (2011) up to 



magnetic fields of 30 T, a great widening on the plateau 
was observed, but it essentially stopped at 5 T. On the 
contrary, high-resolution measurements on the pumped 
current up to 14 T by Fletcher et al. (20111 showed a 



continuous improvement on the pumping accuracy with 
increasing field [see also Fig. |24^b)]. This discrepancy is 
possibly explained by the different samples used in the 
different sets of experiments. 



Giblin et al. (2010) employed a magnetic field of 5 T 
and reported 54 pA of pumped current with la = 15 ppm 
relative uncertainty with a single parameter sinusoidal 
drive. They were able to measure at such a low uncer- 
tainty with a room-temperature current amplifier since 
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they subtracted a reference current from the pumped cur- 
rent and passed less than 100 fA through the amplifier. 
Thus the uncertainty in the gain of the amplifier did not 
play a role. The reference current was created by charg- 
ing a low-loss capacitor and was traceable to primary 
standards of capacitance. 

To date, the most impressive results on single-electron 



lin et al.\ ( 


2012 


). c 


Ref. ( 


Giblin et al. 



2010) 



they made several changes to 
improve the results. They used higher magnetic field of 
14 T and an advanced generation of samples with a litho- 
graphically defined place for the quantum dot in both 
directions in the plane [see Fig. [24|a)]. Instead of using 
sinusoidal waveforms, they also tailored the drive volt- 
age so that the cycle time was distributed more evenly 
for the different parts of the cycle. To make traceable 
measurements, a reference current was created by an ac- 
curate temperature-controlled 1 Gf2 resistor, a voltage 
source and a high-precision voltmeter. The voltmeter and 
the resistor were calibrated through intermediate steps 
against the Josephson voltage standard and the quan- 
tum Hall resistance standard, respectively. In this work. 



Giblin et al. ( 2012 ) reported 150 pA pumped current with 
relative la uncertainty of 1.2 ppm [see Fig. |24^c)]. Most 
of the uncertainty, 0.8 ppm, arose from the calibration 
of the 1 Gfi resistor. Thus it is possible that the the 
electron pumping was actually even more accurate, as 
suggested b y fitting the results to a so-called decay cas- 
cade model (Kashcheyevs and Kaestner 2010[). However, 



there can be processes that are neglected by the model 
and since there is no experimental evidence on lower than 
1.2 ppm uncertainty, it remains the lowest demonstrated 
upper bound for relative pumping errors for quantum dot 
single-electron pumps. Error counting, as demonstrated 



in silicon by Yamahata et al. (2011) and in aluminum 
by |Keller et al.\ ( |1996[ ),~ is a way to measure the pump- 
ing errors to a very high precision independent of the 
other electrical standards and remains to be carried out 
in future for the GaAs quantum dot pumps. 



D. Surface-acoustic-wave-based charge pumping 

After the pioneering experiments on single-electron 
sources based on GaAs discussed in Sec. |III.C.2| the fo- 
cus in this field moved toward the idea of using surface 
acoustic waves (SAWs) to drive the single elect rons in 
a one-dimensional channel (IShilton et al.\ 1996a I. Here, 



the sinusoidal potential created for the electrons in the 
piezoelectric GaAs by SAW forms a moving well that can 
trap an integer number of electrons and transport them 
in a one-dimensional channel. 

The first experiments on this kind of SAW electron 



pumps were carried out by (Shilton et al. 1996b). They 



employed a SAW frequency of 2.7 GHz and observed a 
corresponding n = 1 current current plateau at 433 pA 
with the uncertainty of the order of 10~ 2 at 1 K tern- 
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FIG. 24 (a) SEM image of the device with a schematic mea- 
surement setup employed by Giblin et al. (2012 1, (b) Current 



plateaus obtained by usings a sine wave drive at different fre- 
quencies and magnetic fields, (c) Relative difference of the 
pumped current from ef using a sine waveform and a tailored 
arbitrary waveform at different frequencies. The rightmost 
data point denoted by an asterisk shows the result with the 
potential of the entrance gate shifted by 10 meV from the 
optimal operation point. 



perature. Talyanskii et al. (19971 carried out more de- 



tailed experiments on similar samples at two different 
SAW frequencies and the results were in agreement with 
ef scaling law. Furthermore, several current plateaus 
were observed as a function of the gate voltage corre- 
sponding to different integer values of pumped electrons 
per cycle. However, the experimental uncertainty at the 
plateau was again of the order of 10 -2 and sharp current 
peaks were observed at various gate voltage values. 

After these first experiments, |Janssen and Hartland| 
(2000a|b"t studied the accuracy of the SAW pump and 
reported 431 pA current at the center of the plateau 
with 200 ppm relative deviation from the ideal value. 



Ebbecke et al. ( 2000 ) demonstrated SAW pumping up to 
4.7 GHz frequencies and with two parallel channels to 
increase the current, but the measurement accuracy was 
rather limited here. To improve the quality of the plateau 
Janssen and Hartland ( 2001[ ) decreased the width of the 
one-dimensional channel, which helps in general. How- 
ever, they observed that the required rf power to drive the 
elcetrons increases with decreasing channel width caus- 
ing sever rf heating of the sample. This heating caused 
the quality of the plateau to drop and the conclusion was 
that materials with lower losses due to rf are needed [see 

In fact 



also (Utko et al. 



2006)] 



Ebbecke et al. (20031; 



Flensberg et al. ( 1999 ) reported that the accuracy of the 
SAW current is fundamentally limited in one-dimensional 
channels because of tunneling of electrons out from a 
moving dot. 

To overcome the limitation pointed out in 
Ref. (Ebbecke et al. 2003), the charging energy 



was increased in the system by defining a quantum 



2G 



dot with surface gates rather than utilizing an open 
one-dimensional channel (Ebbecke et al. 2004). Thus 



the applied SAWs modulate both the tunnel barriers 
between the dot and the reservoirs and the electro- 
chemical potential at the dot. With this technique, 
current plateaus were observed at a SAW frequency 
of 3 GHz and the reported relative deviation from the 
ideal value was of the order of 10 -3 . Although the 
results by Janssen and Hartland ( 2000a|b ) remain the 
most accurate ones reported to date with SAW electron 
pumps, and hence are not valuable for a metrological 
current source, single-electron transfer with SAWs can 
be useful in other applications. For example, |McNeil 
et al. (2011 ) have shown that an electron taken by SAWs 



coill ^ (b) 1 cp 



from a quantum dot can be captured by another dot at 
distance. This kind of electron transport can potentially 
be used to transport single spins working as quantum 



bits in a spin-based quantum computer (Hanson et al. 
20071 IMorello et al] 120101). 



E. Superconducting charge pumps 

The envisioned advantage in pumping Cooper pairs in- 
stead of electrons is that the supercurrent produced by 
the Cooper pair pumps is inherently dissipationless and 
the BCS gap protects the system from microscopic exci- 
tations. Thus the operation frequency of the pump can 
possibly be high with the system still remaining at very 
low temperature. Another advantage of the supercurrent 
is that it can sustain its coherence, and hence be virtu- 
ally noiseless, in contrast to single-electron current that 
is based on probabilistic tunneling. Furthermore, since 
the charge of a single Cooper pair is 2e, single-Cooper- 
pair pumps yield twice the current compared with single- 
electron pumps operated at the same frequency. Despite 
these advantages, the lowest uncertainties i n the achieved 
Cooper pair current is at the perce nt level ( Gasparinetti 
et al. |2012 Vartiainen et al. 2007 ) . One reason for this 
is the low impedance of the device rendering it suscepti- 
ble to current noise. 

Two types of Cooper pair pumps exist in the litera- 



ture: arrays of superconducting islands (Geerligs et al. 



1991 ) with source and drain leads, all separated by sin- 



gle Josephson jun ctions with fixed tunnel cou plings, and 
a so-called sluice ( Niskanen et al. 2005 . |2003 ) that com- 
poses of a single island connected to the leads by two 
SQUIDs that function as tunable Josephson junctions, 
see Fig. [25] As in the case of single-electron pumps, the 
device operation is based on Coulomb blockade effects al- 
lowing the controlled transfer of individual Cooper pairs, 
which means in the case of array pumps that the fixed 
Josephson energies of the junctions must be much lower 
than the Cooper pair charging energy of the correspond- 
ing islands. In the sluice, it is sufficient that the min- 
imum obtainable Josephson energy is much lower than 
the charging energy. For the arrays, the thermal energy, 
fcfiT, must be much lower than the Josephson energy that 
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FIG. 25 (a) Scanning electron micrograph of the sluice used 
in the experiments by Vartiainen et al\ ( 2007 1 with a simpli- 
fied measurement setup, (b) Magnified view of the island of 
the device shown in panel (a) with four Josephson junctions, 
(c) Measured pumped current with the sluice (solid lines) as 
a function of the magnitude of the gate voltage ramp such 
that n corresponds to the ideal number of elementary charges 
e pumped per cycle. The inset shows the step-like behavior 
observed in the pumped current. 



defines the energy gap between the ground state and the 
excited state of the quantum system at charge degener- 
acy. For the sluice, the maximum Josephson energy of 
the SQUIDs yields the minimum energy gap of the sys- 
tem thus relaxing the constraint on temperature. 

The first experiment demonstrating Cooper pair 
pumping was performed by Geerligs et al. (1991). The 



device is a linear array of three Josephson tunnel junc- 
tions. The two superconducting islands separated by 
the junctions are capacitively coupled to individual gate 
electrodes. Except in the vicinity of the charge degen- 
eracy points in the gate voltage space, the number of 
Cooper pairs on these islands is rather well defined by 
the gate voltages because the Coulomb blockade regime 
is employed. By biasing the device and applying sinu- 
soidal ac voltages with appropriate amplitudes to the 
gates, one obtains a continuously repeated cycle, during 
which a Cooper pair is transferred through the device, 
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i.e., Cooper pairs are pumped one by one. Ideally, this 
yields a dc current I = 2ef that is proportional to the 
pumping frequency /. The driving voltage at each gate 
should have the same frequency and a phase difference 
of 7r/2. The pumping direction can be reversed if the 
difference is changed by ir. Thus the pumping princi- 
ple is the same as for a normal pump discussed in Sec- 
tion |III.A| The height of the measured current plateau 
follows rather well the predicted relation / = 2ef at low 
pumping frequencies, but deviates strongly at higher fre- 
quencies. This is explained by several mechanisms. The 
uncertainty of the device was not assessed in detail but it 
seems to lie at least on the percent level with picoampere 
currents. One of the error mechanisms is the Landau- 
Zener tunneling when the system is excited to the higher 
energy state without transferring a Cooper pair. This 
was the dominant mechanism at the high-end of the stud- 
ied pumping frequencies in the experiment by |Geerligs| 



et al. 



( 1991 ) thus imposing an upper limit on the opera- 



tion frequency of the device. Another error source in the 
device is the tunneling of non-equilibrium quasiparticles, 
photon-excited tunneling and relaxation of the excited 
states produced by Landau-Zener tunneling. In addi- 
tion, co-tunneling of Cooper pairs through the two junc- 
tions produces a step-like feature to the current plateaus, 
thus reducing the pumping accuracy. Later, a similar 
three-junction C ooper pair pump was studied by |Top-| 
pari et al. { 2004 ) and essentially the same conclusions on 



the pumping accuracy were made. In both experiments, 
no 2e-periodicity was observed in the dc measurements, 
which suggests a substantial presence of nonequilibrium 
quasiparticles in the system. 

The effect of quasiparticles on Cooper pair pump- 
ing was also observed in the seven-junction Cooper pair 
pump (Aumentado et al. 20031. The device is basically 



the same as the one used for pumping single electrons in 



the earlier experiments in the normal state (Keller et al. 



1996). The pump consists of six micrometer-scale alu- 



minum islands linked by aluminum-oxide tunnel barri- 
ers. Investigation of this circuit in the hold and pumping 
modes revealed that besides 2e tunneling events, there 
is a significant number of le events associated with the 
quasiparticle tunneling. All these experiments show that 
in order to obtain accurate Cooper pair pumping, one 



must suppress u n wanted quasiparticle tunneling. Leone 
and Levy (2008); Leone et al. (2008) propose topologi- 
cal protection in pumping Cooper pairs. The charge is 
expected to be strictly quantized determined by a Chern 
index. The idea has not been tested experimentally as 
far as we know. 

In order to increase the output dc current and accu- 
rac y of a single pump, th e sluic e pump was introduced 
by Niskanen et al.\ ( |2005 2003). In the pumping cy- 
cle, the two SQUIDs separating the single island work 
in analogy with valves of a classical pump and the gate 
voltage controlling the island potential is analogous to a 
piston. At each moment of time, at least one SQUID is 
closed (minimum critical current). The gate voltage is 



used to move Cooper pairs through open SQUIDs (max- 
imum critical current). If the pairs are taken into the 
island through the left SQUID and out of the island 
through the right SQUID, the resulting dc current is ide- 
ally / = N2ef, where the number of pairs transported 
per cycle N is determined by the span of the gate volt- 
age ramp. In practise, the critical current of the SQUIDs 
is controlled by flux pulses generated by superconduct- 
ing on-chip coils. Since each operation cycle can transfer 
up to several hundreds of Cooper pairs, |Vartiainen et al.\ 
( 2007 ) managed to pump roughly 1 nA current with un- 



certainty less than 2% and pumping frequency of 10 MHz, 
see Fig. 25 The investigated high-current Cooper pair 



pump demonstrated step-like behavior of the pumped 
current on the gate voltage, however, its accuracy was 
affected by the residual leakage in the tunnel junctions 
and the fact that the SQUIDs did not close completely 
due to unequal Josephson junctions in the structure. 

The leakage current in the sluice can be suppressed by 
working with a phase bias instead of voltage bias which 
was applied in (INiskanen et al. 2005 Vartiainen et al. 



2007 ) . The only experiment reported for a phase biased 
pump was carried out by Mottonen et al. (20081. They 



connected a sluice in a superconducting loop with an- 
other Josephson junction. By measuring the switching 
behavior of this junction from the superconducting state 
to the normal branch with forward and backward pump- 
ing, they were able to extract the pumped current of 
the sluice. However, this type of current detection did 
not turn out to be as sensitive as the direct measure- 
ment with a transimpedance amplifier used in the case 
of voltage bias. A potential way to improve the sensitiv- 
ity is, instead of the switching junction, to use a cryogenic 
current comparator coupled inductively to the supercon- 
ducting loop. This type of an experiment has not been 
carried out to date. Instead, IGasparinetti et al. (2012) 



have measured a sluice in the vicinity of vanishing voltage 
bias, where they demonstrated single-Cooper-pair pump- 
ing plateaus in both the bias voltage and dc level of the 
gate voltage, see Fig. [26j The quasiparticle poisoning 
was reported to be suppressed compared with the pre- 
vious experiments, and hence they observed clear 2ef 
spacing of the current plateaus. 

In addition to the above-mentioned pumping schemes, 



Nguyen et al. ( 2007 1 have studied how a superconducting 



quantum bit referred to as quantronium can be used to 
detect the gate charge ramp arising from a current bias on 
the gate electrode of the island of the device. The accu- 
racy of this technique on converting the bias current into 
frequency remains to be studied in detail. |Hoehne et al.\ 
(2012) studied another type of a quantum bit, a charge 



qubit, for pumping Cooper pairs non-adiabatically. The 
aim here was to increase the pumping speed compared 
to adiabatic schemes but due to the accumulation of the 
errors from a pumping cycle to another, a waiting period 
between the cycles needed to be added. Furthermore, 



Giazotto et al. (20111 have showed experimentally how 



phase oscillations arising from the ac Josephson effect 
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FIG. 26 Pumped average charge by Gasparinetti et al. ( 2012 1 



for a single pumping cycle of a sluice pump near vanishing 
voltage bias as a function of the gate charge offset, n g , and 
span during the pumping cycle, An 9 . 



can drive Cooper pairs in a system with no tunnel junc- 
tions. However, this kind of pumping was found to be 
very inaccurate in this proof-of-the-concept experiment. 



F. Quantum phase slip pump 

There is a proposal to build a source of quantized cur- 
rent based on the effect of quantum phase slip (QPS) 
in nanowires made of disordered superconductors ( Mooij 
2006| ). Phase slip events occur in thin su- 



and Nazarov 



perconducting wires where thermodynamic f luctuations 
of the order parameter become significant (Arutyunov 



et al. 2008). During the phase slip, the superconducting 



order parameter vanishes at a certain instance and po- 
sition in the wire, and the phase difference between the 
wire ends changes by 27r. This gives rise to a voltage 
pulse, in accordance with the Josephson relation. If the 
phase slips happen frequently, they produce a finite dc 
voltage, or a finite resistance. 

Phase slips caused by thermal activation broaden the 
temperature range of superconducting phase transition 
and produce a resistive tail below the critical temper- 
ature of a superconductor (Tinkham 1996). At suffi- 



ciently low temperatures quantum fluctuations take over, 
and the residual temperature independent resistivity of 
a nanowire can be attributed to the quantum phase slips 



(see (Arutyunov et al. 2008) and references therein) 



dissipation associated with every switching event is sup- 
pressed. This can be achieved in superconductors with 
strong disorder, in which Cooper pairs localize before the 



superconducting transition takes place ( Feigel'man et al. 
2007). Such a localization behavior has been observed by 



scanning tunneling microscopy in amorphous TiN and 
InO x films ( Sacepe et al. 2010 2011 ), which are believed 



to be the most promising materials for the observation 
of QPS. 

The key parameter describing a nanowire in the quan- 
tum phase slip regime is the QPS energy -Eqps = ^TqpSj 
where Tqps is the QPS rate. Consider a supercon- 
ducting nanowire of length L, sheet resistance Rq and 
made of a superconductor with the superconducting tran- 
sition temperature T c , coherence length £ = {S.otj 1 / 2 , 
where £o is the BCS coherence length and I the elec- 
tron mean free path (£ <C £o)- Although there is no 
a commonly accepted e xpression for Pqps, it is agreed 
( Arutyunov et al] |2008| |Mooij and Harmans 2005) that: 



TQPS oc exp( 
der unity and Rq 



Thermally activated phase slips are inherently incoher- 
ent. Quantum phase slips may be coherent provided 



), where A is a constant of or- 

irh/e 2 . Clearly, the exponential 
dependence of Tqpg on the wire resistance on the scale 
of £, i?£ = Rut;, requires extremely good control of the 
film resistivity as well as the wire cross-sectional dimen- 
sions. For the nanowires to be in the quantum phase slip 
regime rather than in the thermally activated regime, 
-Eqps should exceed the energy of thermal fluctuations 
k^T. For a typical measurement temperature of 50 mK, 
rqps/27r should be higher than 1 GHz. Although the ex- 
act estimation of R^ is rather difficult, especially in the 
case of strongly disordered films, the experimental data 
presented by (Astafiev2012) for InO x films agrees with 
the following values: Fqps/27t ps 5 GHz, R^ = 1 kfi and 
£ = 10 nm. 

The first experiment reporting the indirect observa- 
tion of coherent QPS in nanowires was performed by 
(Hongisto2012). They studied a transistorlike circuit 
consisting of two superconducting nanowires connected 
in series and separated by a wider gated segment. The 
circuit was made of amorphous NbSi and embedded in 
a network of on-chip 30 nm thick Cr microresistors en- 
suring a high external electromagnetic impedance. The 
NbSi film had a superconducting transition temperature 
of ps 1 K and normal-state sheet resistance about 550 
O per square. Provided the nanowires are in the regime 
of QPS, the circuit is dual to the dc SQUID. The sam- 
ples demonstrated appreciable Coulomb blockade voltage 
(analog of a critical current of the dc SQUID) and pe- 
riodic modulation of the blockade by the gate voltage. 
Such a behavior was attributed to the quantum interfer- 
ence of voltages in two nanowires that were in the QPS 
regime. This is completely analogous to the quantum 
interference of currents in a dc SQUID. 

An unambiguous experimental evidence of coherent 
QPS was provided in the work by (Astafiev2012). Co- 
herent properties of quantum phase slips were proven by 
a spectroscopy measurement of a QPS qubit, which was 
proposed earlier by (Mooij2005). The qubit was a loop 
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that had a 40-nm-wide and about l-/im-long constric- 
tion. The loop was made of a 35-nm-thick supercon- 
ducting disordered InO x film with T c = 2.7 K and a 
sheet resistance of 1.7 kSl per square slightly above T c . 
The qubit was coupled inductively to a step-impedance 
coplanar waveguide resonator, which was formed due to 
the impedance mismatch between an indium oxide strip 
and Au leads to which it was galvanically connected. The 
ground planes on both sides of the strip were made of Au. 
At the qubit degeneracy point at a flux bias (m+l/2)5>0j 
where m is an integer, there is an anticrossing in the 
qubit energy spectrum with a gap Sqps = HTqps- At 
this flux bias, the two quantum states \m) and |m + 1) 
corresponding to the loop persistent currents circulating 
in the opposite directions are coupled coherently, which 
gives a gap -Eqps between the lowest energy bands of the 
qubit. With the flux offset 5$ from degeneracy, the gap 
evolves as AE = ((2/ p <S$) 2 + £q PS ) 1/2 , where I p is the 
persistent current in the loop. This gap was revealed in 
the spectroscopy measurements by monitoring the res- 
onator transmission as a function of external magnetic 
field and microwave frequency. When the microwave fre- 
quency matched the qubit energy gap, a dip in the trans- 
mission was observed. The width of the dip ~ 260 MHz 
close to the degeneracy point indicated rather strong de- 
coherence whose origin is still to be understood. 

Based on the exact duality of QPS and the Josephson 
effects, it is argued that it should be possible to build a 
QPS electric current standard, which is dual to the ex- 



isting Josephson voltage standard (Mooij and Nazarov 
2006). When biased resistively and irradiated by a high 
frequency signal, QPS junctions should exhibit current 
plateaus, which could provide the basis for the funda- 
mental standard of the electric current. When an ac sig- 
nal of frequency / is applied to a Josephson junction, 
Shapiro voltage steps V n = n^, where n is an integer, 
are observed. Similarly, when an ac signal is applied to a 
QPS junction, an equivalent of Shapiro steps will occur 
in form of plateaus at constant current levels I n = n2ef. 
One should note, however, that error mechanisms have 
not yet been analyzed for this type of quantized current 
source: thus it is not clear at the moment how accurate 
this source will be. 

From the practical point of view, realization of a QPS 
current source looks rather challenging because it re- 
quires fabrication of nanowires with an effective diameter 
~ 10 nm as well as precise control of the sheet resistance 
Ra of the nanowire, which is in the exponent of the ex- 
pression for -Eqps- Various approaches to the nanowire 
fabrication including step decoration technique, sputter- 
ing of a superconductor on a suspended carbon nanotube, 
trimming of a nanowire by argon milling, etc. are de- 
scribed in Ref. (Arutyunov et al. 20081. Another issue 



nanowire with the cross-sectional dimensions 20 nm x 
20 nm, sheet resistance 1 kO per square and carrying a 
dc current of 100 pA will have the effective electron tem- 
perature of the order of 250 mK, which is high enough 
to smear the current plateaus. 

One of the first attempts to observe current plateaus 
on the current- voltage characteristics of superconducting 
nanowires under the rf radiation was reported by ( Lehti- 
nen et aT\ 2012). The nanowires were made of Ti and 
had length up to 20 /Ltm and effective diameter from 40 
nm down to about 15 nm. The nanowire sheet resistance 
varied from about 20 ft up to 1.9 kfl per square. They 
were biased through high-Ohmic Ti or Bi leads having 
total resistance of 15 kfl and 20 Mf2, respectively. The 
low-Ohmic samples biased through 15 kil exhibited weak 
Coulomb blockade. The estimated -Eqps was ~ 0.1 /zeV 
only. More resistive nanowires (R\j = 180 ft, effective 
diameter ps 24 nm) biased through 20 Mfl leads had pro- 
nounced Coulomb blockade with a critical voltage of up 
to 0.4 mV. The thinnest nanowires (Rq — 1.9 kfi, effec- 
tive diameter < 18 nm) exhibited a Coulomb gap of a 
few hundred mV with the largest gap exceeding 600 mV. 
These gaps did not vanish above T c of Ti, from which 
the authors concluded that some weak links were unin- 
tentionally formed in the thinnest nanowires. Despite 
the fact that nanowires had large variations of param- 
eters, all their current-voltage or dV/dl characteristics 
exhibited some quasiregular features under the external 
rf radiation. Those features were interpreted as being 
current steps formed due to the phase-locking of intrin- 
sic oscillations by the external signal. 

It is interesting to note that physics of QPS in su- 
perconducting nanowires resembles physics of QPS in 
Josephson junction arrays ( Fisher| |1986[ ). A nanowire 
can be modeled as a ID array of small superconduct- 
ing islands connected by Josephson junctions. Forma- 
tion of isolated superconducting regions within a nomi- 
nally uniform disordered film was confirmed experimen- 



tally ( Sacepe et al. , 2010 2011 ). Such a weakly connected 



is the overheating of the nanowire electron system. As- 
suming the phase slip region becomes normal (which is 
true, for example, for Ti nanowires) , for the estimation of 
the electron temperature, one can use the expression ( |30[ ) 
for the power transfer from electrons to phonons. A 



array of superconducting islands is characterized by the 
junction Josephson energy Ej and the island charging 
energy E c . The phase and charge dynamics of the ID 
array depends on the ratio Ej / E c . In the experiment by 
(Pop et al. 2010[ ) Ej/E c in a SQUID array was tuned 
in situ by applying a uniform magnetic flux through all 
SQUIDs. The state of the array was detected by an extra 
shunt Josephson junction. The authors deduced the ef- 
fect of the quantum phase slips on the ground state of the 
array by measuring the switching current distribution of 
the entire Josephson circuit as a function of the external 
magnetic flux for different values of Ej/E c , 



G. Other realizations and proposals 

In this subsection we cover various ideas that have been 
brought up fx experimental demonstration. Although 
their metrological relevance is still to be proven, we wish 
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to present them for the sake of their complementarity, 
potential and for completeness. 



1. Ac-current sources 



2eV (t) 



The current pumps described in sections |III.A| |III.B| 
and |III.C| can be considered as single-electron injectors 
generating dc current. Coulomb blockade ensures a good 
control of the electron number on each island during the 
charge transfer. However, there is no time control of the 
electronic injection in metallic and hybrid single-electron 
pumps because of the stochastic nature of tunneling. 
Time-controlled injection can be realized in quantum 
dots, however, the energy of emitted electrons spreads 
randomly in a wide range exceeding hT, where T is the 
electron tunneling rate. 

A time-controlled sing le-electron source generating ac 
current was reported by |Feve et al. (2007). The source 



was made of a quantum dot, realized in a GaAlAs/GaAs 
heteroj unction of nominal density n s = 1.7 x 10 15 m~ 2 
and mobility fj, — 260 V m s and tunnel- coupled 
to a large conductor through a quantum point contact 
(Fig. 27). The discrete energy levels of the quantum dot 
were controlled by the voltage applied to the capacitively 
coupled metallic gate located 100 nm above the 2DEG. 
When the applied voltage is high enough, the electron 
occupying the highest energy level tunnels to the con- 
ductor. This process is called electron emission. If the 
applied voltage is changed abruptly at the subnanosec- 
ond time scale (a sudden jump), the emission process is 
coherent and the final state of the electron is a coherent 
wave packet propagating away in the conductor. The en- 
ergy width of the wave packet is given by the inverse tun- 
neling time, as required for an on-demand single-particle 
source. It is also independent of temperature. Its mean 
energy is adjusted by the amplitude of the voltage step, 
V^ xc . For all measurements, the electronic temperature 
was about 200 mK. A magnetic field B s» 1.3 T was ap- 
plied to the sample so as to work in the quantum Hall 
regime with no spin degeneracy. The QPC dc gate volt- 
age Vq controlled the transmission D of a single edge 
state as well as the dc dot potential. As reported by 



Gabelli et al. ( 2006 1 , this circuit constitutes the effective 



quantum-coherent RC circuit, where coherence strongly 
affects the charge relaxation dynamics. The effective 
quantum resistance R and capacitance C are defined as 
R = h/2e 2 and C = e 2 {dN/de), where dN/de is the lo- 
cal density of states of the mode propagating in the dot, 
taken at the Fermi energy ( TPretre et al] |1996[ ). Inter- 
estingly, the QPC transmission probability D affects the 
quantum capacitance but not the quantum resistance. In 
the experiment, the addition energy A + e 2 /C rj 2.5 K 
was extracted. The large top-gate capacitance makes the 
Coulomb energy e 2 /2C unusually small, and the total ad- 
dition energy reduces to the energy- level spacing A. 

The single-charge injection was achieved by the appli- 
cation of a high-amplitude excitation voltage V cxc ~ A/e 
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FIG. 27 Schematic of the single-charge injector and its oper- 
ation principle. Starting from step 1 where the Fermi energy 
level of the conductor lies in between two energy levels of the 
dot, the dot potential is increased by A moving one occupied 
dot level above the Fermi energy (step 2). One electron then 
escapes from the dot. After that the dot potential is brought 
back to the initial value (step 3), where one electron can en- 
ter the dot leaving a hole in the conductor. One-edge channel 
of the quantum RC circuit is transmitted into the dot, with 
transmission D tuned by the QPC gate voltage Vg- The dot 
potential is controlled by a radio-frequency excitation V BXC ap- 



plied on a m acroscopic gate located on top of the dot (Feve 
et aZ.[|2007l. 



to the top gate. When an electron level is brought sud- 
denly above the Fermi energy of the lead, the electron 
escapes the dot at a typical tunnel rate r _1 = DA/h, 
where A/h is the attempt frequency. Typically, the 
tunnel rates are in the nanosecond time scales, and 
this makes the single-shot charge detection a challeng- 
ing task. In the experiment, a statistical average over 
many individual events was used by repeating cycles 
of single-electron emission with period T followed by 
single-electron absorption (or hole emission) as shown in 
Fig. |27| This was done by applying a periodic square 
wave voltage of amplitude « A/e to the top gate. Thus, 
the signal-to-noise ratio was increased. 

Under certain conditions single-electron repeatable in- 
jection leads to a good quantization of the ac current. \I U \ 
as a function of V exc for two values of the dc dot potential 
at D s» 0.2 and D ps 0.9 is shown in Fig. [28] Apparently, 
a good current plateau is expected when the charge on 
the dot is well defined. The latter depends on the trans- 
mission and the working point set by the gate voltage. 
For example, transmission D ss 0.2 is low enough and the 
electronic states in the dot are well resolved, as shown in 



the inset of Fig. 28 (left). On the other hand, the trans- 



mission is large enough for the escape time to be shorter 
than T/2. When the Fermi energy lies exactly in the 
middle of the density-of-states valley (blue vertical line 
in the left inset)), a well-pronounced | = 2e/ current 
plateau is observed centered at 2eV cxc /A — 1. At this 
working point, the plateau is rather flat. It is claimed 
that the current uncertainty at the plateau is 5% due 
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to the systematic calibration error. In contrast, if, with 
the same transmission, the Fermi energy lies on the peak 
(light blue vertical line in the left inset), there is still a 
current plateau, but it is not as flat and very sensitive 
to parameter variations. When the energy level is res- 
onant with the Fermi energy, the dot mean occupation 
at equilibrium is 1/2 and the height of the plateau is 
1/2 x 2ef = ef (Fig. [28] A, left), which makes this work- 



ing point unsuitable for single-electron pumping. When 
transmission is increased, the charge fluctuations become 
stronger and the plateau width gets smaller and finally 
nearly vanishes at D ai 0.9 even for the optimal work- 
ing point, as seen from Fig. 28 A, right). The pump- 



ing curves measured at different transmissions cross at 
2eV exc /A — 1, which corresponds to the pumped cur- 
rent \I W \ = 2ef. Finally, Fig. 28 B (top) shows domains 
of good charge quantization, which are presented as 2D 
color plots, as a function of the excitation voltage V exc 
and the gate voltage Vq- The white diamonds are the 
domains of constant current |7 W | = 2ef suitable for AC 
current pumping. At high transmissions, the diamonds 
become smeared by the dot charge fluctuations as dis- 
cussed previously. At small transmissions, even though 
the dot charge quantization is good, current quantization 
is lost because of the long escape time tot 1, and the 
current vanishes. At 180 MHz, the optimal pumping is 
observed at D as 0.2. In Fig. [28] the experimental results 
are compared with the theoretical model (solid lines in 
Fig. 28 A and lower plot in Fig. 28 B, respectively) (ID 
modeling of the circuit (density of states, transmission, 



200G) 



and dot-gate coupling) described in ( Gabelli et al 
was used) , showing excellent agreement between the two 
The device described above is the electron analog of 
the single-photon gun. It is not a source of quantized dc 
current as the dot emitting the electron can be recharged 
only though the reverse process of the electron absorp- 
tion. However, periodic sequences of single-electron emis- 
sion and absorption generate a quantized alternating cur- 
rent |7J »2e/. 



2. Self-assembled quantum dots in charge pumping 

The idea of using self- assembled quantum dots for 
charge pumping is based on conversion of optical excita- 
tion into deterministic electric current, see |Nevou et aL] 
(2011). Unlike in the quantum dots formed in 2DEG 
that were described in |III.C[ in self-assembled QDs the 
electrons are confined not due to the Coulomb blockade 
effect, but because of the energy level quantization in 
the ultrasmall island. For the well-studied InAs/GaAs 
quantum dots, the characteristic energy scale, i.e., the 
energy difference between the ground state and the ex- 
cited state, is appro ximately 55 meV and can be tuned 
to even larger values ( Liu and Hopkinson 2003 1 . As a re- 




FIG. 28 AC current quantization. (A) |7 TO | as a function 
of 2eK xc /A for different dot potentials at D as 0.2 (left) and 
D as 0.9 (right). Points are measured values and lines are the- 
oretical predictions. Insets show schematically the dot density 
of states N(e). The color vertical lines indicate the dot poten- 
tial for the corresponding experimental data. (B) Measured 
(top) and modeled (bottom) | I w \ as a function of 2eV C xc/A 
and Vg shown as color plots ( Fcv e et q£) |2007[ ) . 



given operation condition. 



In the experiment of Nevou et al. (2011) a plane of 



self- assembled InAs quantum dots is coupled to an In- 
GaAs quantum well reservoir through an Alo.33Gao.67As 
barrier (see Fig. 29'a)-(c)). The structure is sandwiched 
between two n-doped GaAs regions. The device basi- 
cally works as a strongly asymmetric quantum dot in- 
frared photodetector (Nevou et aLl|2010|). In the absence 



suit, an increase in the electron confinement energy from 
the values ~ 1 meV typical for Coulomb blockade devices 
to 55 meV should improve the potential accuracy for a 



2010]) 

of any optical excitation, electrical conduction is inhib- 
ited by the AlGaAs barrier. When a laser pulse ionizes 
the quantum dots, electrons are excited out of the dot 
and then swept away by the applied bias voltage, giv- 
ing rise to a photocurrent. Electrons are then collected 
by the contact and thermalized, contributing to the cur- 
rent. Once all of the electrons are ionized, they will be 
refilled from the electron reservoir by tunneling through 
the AlGaAs barrier. The refilling time r can be varied 
over many orders of magnitude by changing the AlGaAs 
barrier width or the Al concentration. As a result, when 
the pulse energy is high enough to completely ionize the 
dot and its duration t p is much smaller than r, a fixed 
number of electrons will be photo-excited. If the process 
is repeated at a frequency / , the current will be given by 
equation 7 = nef , where n is determined by the number 
of dots and the number of electrons per dot. 

In order to suppress possible sources of error affecting 
the accuracy of the generated photocurrent, the proba- 
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FIG. 29 (a) Schematic layout of the self-assembled-quantum- 
dot electron pump; (b) transmission electron microscopy im- 
age of the quantum dots; (c) 3D sketch of the conduction band 
profile of the structure under zero bias; (d) Saturation current 
for two different pump wavelengths (A = 6.7/xm: curve A and 
A = 10fj,m: curve B. The difference provides a current plateau 
that should be 2ef (thick horizontal line). Inset: Variations 
of the measured current with respect to the average value 
UNevou et aL| 12010b. 



saturate the quantum dots during one pulse. The gen- 
erated current was measured continuously as a function 
of the bias voltage. The saturation current depended on 
the bias voltage, power and pulse wavelength, as shown in 
Fig. 29 [d) , which is not suitable for the quantized current 
source. By subtracting the curve measured for A = 10/im 
from the one measured for A = 6.7/xm, a current plateau 
is obtained in the voltage range - 150 mV to - 1.07 V. By 
rescaling the current in electron flux per dot (taking into 
account the mesa size, the measured dot density and the 
repetition frequency of the laser pulses) , one obtains the 
current plateau at the absolute level of 1.8 electrons per 
dot instead of the expected value of 2 electrons per dot, 
giving a pumping accuracy of 10%. This is explained 
by the uncertainty in the dot density of 13.6 % mea- 
sured by the transmission electron microscopy. Using the 
same subtraction method, the dependence of the current 
plateau height was measured as a function of the pulse 
repetition frequency in the range 10 Hz to 100 kHz. As 
expected, a linear dependence was observed. The proce- 
dure of subtracting the pumped currents measured at dif- 
ferent wavelengths was used to compensate for the strong 
inhomogeneous broadening of the quantum dot sizes and 
must be avoided in a real metrological device. For this, 
a more uniform set of dots must be fabricated. 



3. Mechanical single-electron shuttles 



bility that some dots are doubly ionized during one op- 
tical pulse and the probability that not all dots are filled 
during the refilling phase, one needs to satisfy two con- 
ditions: t p /r < e and exp(— 1/r/) < e, where e = 10~ 8 
is the desired accuracy. Both conditions are easily sat- 
isfied in the experiment. It follows from the two con- 
ditions above that / should not exceed kilohertz range. 
Although this seems a relatively low frequency, one must 
keep in mind that in the proposed device millions of dots 
are running in parallel, so even this rate could lead to 
the generation of currents in the nanoampere range. For 
example, in the performed pumping experiment, a 210 x 
210 /im 2 mesa was fabricated by photolithography, and 
with a dot density of the order of 10 9 cm~ 2 it contained 
about 44 million quantum dots. The final source of er- 
ror, independent of the device design, arises from the 
uncertainty in the number of quantum dots contained 
in one device, as well as the variability in the quantum 
dot transition energy. It is claimed that this is more of 
a technological problem than a fundamental limit and a 
variety of solutions can be found to solve this problem. 



In the pumping experiment by Nevou et al. (2011 1, re- 



filling time r varying from 3 to 5.5 fj,s was measured. For 
the excitation of the dots and generation of the current, 
a laser-pumped optical parametric system with a wave- 
length A continuously tunable between 3 and 16 fiia was 
used. It delivered 100-ps-long pulses with a repetition 
rate of 1 kHz. The laser source power was sufficient to 



Besides charge pumps with entirely electronic control, 
there is a group of devices in which a mechanical degree 
of freedom is involved. They are called mechanical charge 
shuttles, because they transfer either single charges (elec- 
trons or Cooper pairs) or portions of charges between the 
two electrodes due to the mechanical back and forth mo- 
tion of a small island between the electrodes. This results 
in current flow, either incoherent or coherent. The con- 
cept of the mechanical electron shuttle was introduced in 
1998 by (|Gorelik et oT| |1998| |Isacsson et al.||1998[). 



The proposed device has a small conducting island, 
which is mechanically attached to electrical leads with a 
help of an elastic insulator, see Fig. 30 'a) . The dc volt- 



age applied between the leads and elastic properties of 
the insulator together with charging and discharging of 
the island create instability and make the island oscil- 
late (see Fig. 30 'b)). For the proper operation of the 
shuttle, two assumptions were made: the amplitude of 
the mechanical oscillations is much larger than the elec- 
tron tunneling distance and the number of electrons on 
the island is limited. With these assumptions the island 
motion and charge fluctuations become strongly coupled. 
As the shuttle comes into contact (direct or tunneling) 
with an electrode, it will pick up charge q and accelerate 
in the electric field toward the opposite electrode. Upon 
contact with this second electrode, the shuttle transfers 
charge and acquires an equal but opposite charge —q. 
The shuttle then accelerates back towards the first elec- 
trode to begin the cycle anew, as shown in Fig. 30 'c). In 
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ity was reported by Tuominen et al. (1999). This was 
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FIG. 30 (a) Model of the shuttle device proposed by (Gore 



lik et al. 19981, where a small metallic grain is linked to two 



electrodes by elastically deformable organic molecular links; 
(b) Dynamic instabilities occur since in the presence of a suf- 
ficiently large bias voltage V the grain is accelerated by the 
electrostatic force towards the first electrode, then towards 
the other one. A cyclic change in direction is caused by the 
repeated loading of electrons near the negatively biased elec- 
trode and the subsequent unloading of the same charge at the 
positively biased electrode. As a result, the sign of the net 
grain charge alternates leading to an oscillatory grain motion 
and charge transport; (c) Charge variations on a cyclically 
moving metallic island. The force exerted on the island by 
the electrostatic field along the island displacement between 
source and drain creates instability making the island move. 
The dashed lines in the middle describe a simplified trajec- 
tory in the charge-position plane, when the island motion by 
Sx and discharge by 2Sq occur instantaneously. The solid 
trajectory describes the island motion at large oscillation am- 
plitudes. Periodic exchange of the charge 2q = 2CV + 1 
between the island and the leads results in charge pumping. 
From dGorelik et al.\ 119981) 



this manner, a net shuttle current of I = 2qf is estab- 
lished, where / is the shuttle frequency and q the charge 
transferred in one cycle. The dc voltage applied between 
the source and drain electrodes is sufficient to drive the 
island. There is a Coulomb force pushing the island pe- 
riodically to one lead or the other. Depending on the 
shuttle details, two regimes can be distinguished: classi- 



cal (Gorelik et al. 1998 Isacsson et al. 1998 Weiss and 



Zwerger 1999) and quantum mechanical (Armour and 



MacKinnon} |2002| |Cohen et al.\ [20091 |Fedorets et~aJ. 



2004 Johansson et al. 2008). The shuttle, when made 



superconducting, can transfer not only electrons but also 



Cooper pairs (Gorelik et al. 2001 Shekhter et al. 2003). 



The first experimental realization of the mechanical 
charge shuttle that operated due to a shuttle instabil- 



a rather bulky device even though it was scaled down 
considerably in size and operating voltage in compari- 
son to the earlier electrostatic bell versions. The mea- 
sured device contained a resonating beam of dimensions 
40 mm x 22 mm x 1.6 mm and a shuttle of mass 0.157 g 
and radius of about 2 mm. The effective mass of the 
bending beam was 30 g, its fundamental vibrational fre- 
quency 210 Hz and quality factor 37. The shuttle at- 
tached to the beam was placed in between two voltage 
biased massive electrodes having a gap adjustable by a 
micrometer. Besides measuring the current with a sensi- 
tive electrometer, the authors also measured the shuttle 
frequency using a diode laser and phototransistor sensor 
in a photogate arrangement. Both current and frequency 
were measured as a function of voltage. The observed 
current jumps as a function of the bias voltage as well as 
hysteresis in the transport characteristics were the main 
indications of the shuttling regime of the device. 

A nanoscale version of the instability-based electron 
shuttle was implemented by Kim et al. (2010), see 



Fig. |3l|a)). The device consisted of a 250-nm-tall Si pil- 
lar with diameter of 60 nm, covered on top with a 45 nm 
thick gold layer. The pillar was placed in the middle of a 
110 nm gap between two electrodes, source and drain, of 
the central line of a coplanar waveguide. For reference, 
the authors also fabricated and measured a similar device 
without a pillar in the gap. The samples were measured 
at room temperature in a probe station evacuated to a 
pressure of 10 -4 mbar. The pilar was actuated by ap- 
plying a small rf signal together with a dc bias voltage 
across the source and drain electrodes. The amplitude of 
the rf signal alone was not sufficient to make the pillar 
resonate, however, it was sufficient for driving in combi- 
nation with the dc bias. A clear frequency dependence 
was observed for the sample with a pillar in the gap, with 
the resonance frequency of 10.5 MHz and quality factor 
of about 2.5. It was estimated that the device shuttles 
on average 100 electrons per cycle. The resonance was 
also detected in the mixing experiment. 

Another realization of the nanoelectronic shuttling de- 



vice was reported by Moskalenko et al.\ ( |2009a|b[ ) . It had 
a configuration of a single-electron transistor, whose is- 
land was a 20-nm large gold nanoparticle placed in be- 
tween the Au source and drain electrodes by means of 
an atomic force microscope (AFM). First, planar 30- 
nm-thick gold electrodes were made by electron-beam 
lithography and separated by a gap of 10-20 nm. Then 
they were covered by flexible organic molecules, 1,8- 
octanedithiol, with a length of about 1.2 nm. The 
nanoparticle was linked to the electrodes through flexi- 
ble organic molecules, 1,8-octanedithiol, with length of 
about 1.2 nm. Then nominally 20-nm-diameter gold 
nanoparticles were adsorbed by immersion of the elec- 
trode assembly into aqueous gold solution. After that, an 
AFM image of the sample was taken and one nanopar- 
ticle in the vicinity of the gap was identified. It was 
pushed to the gap by the AFM tip. Current-voltage 
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FIG. 31 (a) Scanning electron micrograph of a nanopillar be- 
tween two electrodes ( |Kim et al.\ 2010| |; (b) Electron micro- 
graph of the quantum bell: The Si beam (clapper) is clamped 
on the upper side of the structure. It can be actuated by ac 
voltage applied to the gates on the left- and right-hand sides 
(Gl and G2) of the clapper (C). Electron transport is then 
observed from source (S) to drain (D) through the island on 
top of the clapper. The island is electrically isolated from the 
rest of the clapper which is grounded (Erbe et al. 2001); (c) A 



false-color SEM image taken at an angle to reveal the three- 
dimensional character of the nanomechanical SET. A gold 
(yellow) island is located at the center of a doubly clamped 
freely suspended silicon nitride (red) string. The gold island 
can shuttle electrons between the source and drain electrodes 



when excited by ultrasonic waves ( Koenig et al. 2008 1 



characteristics of the fabricated devices were measured 
at room temperature in a shielded dry box to protect 
the structure from moisture and electromagnetic noise. 
The authors observed characteristic current jumps in the 
current-voltage curves and attributed them to the shut- 
tling effect. They also compared the transport character- 
istics of their devices with and without the nanoparticle. 
They also compared characteristics of a working shuttle 
device with the same device but from which the nanopar- 
ticle was removed. After this procedure, current through 
the device dropped below the noise level. 

The effect of the mechanical vibrational modes on 
charge transport in a nanoelectronic device was observed 
in the Ceo single-electron transistor ((Park et al. |2000 )) . 
In this device, a single Cgo molecule was placed in the 
narrow gap between the two gold electrodes. It was found 
that the current flowing through the device increases 
sharply whenever the applied voltage was sufficient to 
excite vibrations of the molecule. Although mechanical 
vibrations were observed in this structure, they are not 
related to the shuttling of electrons. 

Externally driven resonant shuttles may be easier to 
implement in comparison to the instability-based shut- 
tles, because a much larger displacement amplitude can 
be achieved. However, the drawback of using a mechani- 
cal resonator is the discrete set of eigenfrequencies, which 
are determined by geometry and materials. Therefore, 
only a limited number of frequencies are available for 
electron transfer. In the experiment performed by |Erbe[ 



et al. ( 1998 ) a Si nanomechanical resonator was placed in 



between two contacts. The whole device, being a scaled- 
down version of the classical bell, was functioning as a 
mechanically flexible tunneling contact operating at radio 
frequencies. The sample was machined out of a single- 
crystal silicon-on-insulator substrate by a combined dry- 
and wet-etch process. It consisted of a metallized singly- 
clamped beam (the clapper), whose free end is aligned 
between the source and the drain. The beam was 1 /im 
long, 150 nm wide and 190 nm thick. The beam was 
driven by an oscillating voltage applied on two gates at 
frequencies up to 100 MHz. There was a ir phase shift be- 
tween the two gate voltages. The clapper was biased with 
a dc voltage and the source was grounded. The resulting 
clapper-drain dc current as well as rf response were mea- 
sured at 300 K and 4.2 K. The current-frequency depen- 
dence contained strong peaks, which were interpreted as 
being due to the mechanical resonances of the beam, in- 
dicating that shuttling was occurring. The peaks had low 
quality factors, ranging from 100 to 15, only. The num- 
ber of electrons N shuttled per cycle was estimated from 
the current peak height / using the relation N — I/ef. 
Below 20 MHz, 10 3 -10 4 electrons were shuttled in each 
cycle. On the 73 MHz peak the number was decreased 
to about 130 electrons per cycle. 



One issue in the device presented by Erbe et al. (1998) 



is that there was no island in it, which makes this system 
different from the one envisaged by Gorelik et al. (1998 1 . 



This issue was addressed in the later work by the same 
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group ( |Erbe et al. 2001). Using a similar fabrication 
procedure, the authors fabricated a singly-clamped beam 
(clapper) with a metal island on its end (see Fig.pHfTj)). 
First, transport characteristics of the device were mea- 
sured at 300 K. It was found that there was no detectable 
current through the device unless a driving ac voltage (± 
3 V) was applied to the driving gates. Under the exter- 
nal drive, the current exhibited several peaks, similar to 



those in the earlier device (Erbe et al. 1998), which was 



attributed to the beam motion. The background current 
was explained by the thermal motion of the beam. At 
a temperature of 4.2 K all the current peaks were sup- 
pressed except one at about 120 MHz with much smaller 
height (only 2.3 pA). This corresponds to shuttling of an 
average of 0.11 electrons per cycle. 



Koenig et al. (2008) implemented electromechanical 



single-electron transistors (MSETs) with a metallic is- 
land placed on a doubly clamped SiN beam (Fig. 31 c)). 



An array of 44 mechanical SETs was measured in one run. 
The sample was placed in a vacuum tube filled with He 
exchange gas at a pressure of 7.5 x 10~ 4 mbar. The mea- 
surements were performed at a temperature of 20 K. The 
beams were actuated using ultrasonic waves produced by 
a piezo unit. The actuation frequency ranged from 3.5 up 
to 6 MHz. To shield the sample from the electrical sig- 
nal required for the piezo actuation, the sample chip was 
placed inside a Faraday cage with solid 3-mm titanium 
walls. Measurements revealed several resonance features 
in the SET dc current, which were attributed to the me- 
chanical resonances. It was argued that the mechanical 
motion of the resonator was strongly nonlinear. This was 
imposed by the side electrodes constituting the impact- 
ing boundary conditions. The nonlinear nature of the 
system resulted in the shape of the resonance curves dif- 
ferent from Lorentzian. Although the expected step-like 
dependence of the SET current on the source-drain volt- 
age was not observed (because of the high measurement 
temperature as compared to the charging energy), the 
authors made an optimistic conclusion, that the device 
may be useful for quantum metrology. 



4. Electron pumping with graphene mechanical resonators 

An electron pump based on a graphene mechanical res- 
onator in the fundamental flexural mode was introduced 
by Low et al. (2012 1. The resonator is actuated elec- 



trostatically by a gate electrode. Time varying deforma- 
tion of graphene modifies its electronic energy spectrum 
and in-plain strain. Cyclic variation of these two prop- 
erties constitutes the scheme for quantum pumping. To 
have a nonzero pumping current, spatial asymmetry must 
be introduced. It is assumed that the contacts between 
the graphene layer and the left and right electrodes are 
not equivalent, which is modeled by different densities of 
states. This can be achieved in the experiment by using 
different materials for the two electrodes. It is empha- 
sized that Coulomb blockade effects will favor the trans- 



fer of an integer number of electrons per cycle, so that 
the relation between current and frequency will be quan- 
tized. This is just a proposal and the applicability of this 
approach for quantum metrology is still to be verified. 



5. Magnetic field driven single-electron pump 

Another proposal, not implemented though, is based 
on using for electron pumping a ferromagnetic three- 
tunnel-junction device ( jShimada and Ootuka[ 2001). Its 
islands and leads are made of ferromagnetic metals with 
different coercive forces. Such a device can be operated as 
a single-electron pump if controlled not by the gate volt- 
ages, but by ac magnetic fields. In addition to the charg- 
ing effects, it makes use of the magnetic-field-induced 
shift of the chemical potential and magnetization rever- 
sal in the ferromagnetic electrode. 

The proposed device has intrinsic limitations of the 
pumping speed, which are determined by the physical 
time constraints of the ferromagnet. The pump operation 
frequency must be much lower than the characteristic re- 
laxation times. The prospects of this type of an electron 
pump for quantum metrology are still to be understood. 



6. Device parallelization 



As discussed in Sect. HI. A it is possible to reach pre- 
cise electron pumping with a normal metal single-electron 
pump consisting of a sufficiently long array of islands. 
With six islands and seven junctions, the accuracy of the 
pumped current is at 10 -7 level. However, the maximum 
current is limited to a few pA. To get the current scaled 
up to 100 pA level, a requirement for practical metrologi- 
cal applications, see Sect. |IV.D.l[ approximately hundred 
pumps should be operated in parallel. The main reason 
why parallelization is impractical for normal-metal de- 



vices is the tuning of the offset charges ( Camarota et al. 



2012 Keller et al. 1996). Each island has an individ- 



ual offset charge that has to be compensated separately. 
Therefore, a metrological current source implemented as 
parallellized normal metal pumps would require of the 
order of 1,000 dc lines. 

Compared to normal metal pumps, quantum dot-based 
pumps allow for higher pumping speeds using fewer con- 
trol lines thanks to the tunability, see Sect jllLC Accu- 
racy of 1.2 ppm at output current of 150 pA has been 
already demonstrated with a single quantum dot ( |Giblin| 



et al. 2012). Therefore, parallelization of such pumps 



may not be even required if the accuracy can be im- 
proved without a loss in speed. Nevertheless, paralleliza- 
tion of semiconducting pumps has been considered in 
the literature. With two pumps, invariance with respect 
to gate variations has been shown to be below 20 ppm 
level ( Wright et al. 2009 1 with output current exceeding 
100 pA. In this case, all signals were individually tuned 
for each device requiring two dc and one rf signal per de- 
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FIG. 32 (a) Scanning el ectron micrograph of parallel turn- 



stiles (Maisi et al. 2009). The tunstiles are biased with a 
common bias 14 and driven with a common rf gate voltage 
Vrf. Gate offset charges are compensated by individual gate 
voltages Vg,i. (b) Output current I for ten parallel devices 
tuned to the same operating point producing current plateaus 
at I = WNef. The curves are taken at different 14 shown in 
the top left part of the panel. 



vice. However, it is possible to use common signals for rf 
drive and for the barrier voltages (Mirovsky et al. 20101. 



In this case, only one dc voltage per device is required for 
tuning the other barrier and possible offset charges. The 
obtainable accuracy, depending on device uniformity, is 
still an open question for this approach. 

For the hybrid NIS turnstiles, the maximum current 
per device is limited to few tens of picoamperes, as dis- 
cussed in Sect. |III.B.2| Hence, at least ten devices are to 
be run in parallel, which has been shown to be experi- 
mentally feasible ( Maisi et al. . |2009 ) . In Fig. 32 we show 
a scanning electron micrograph of a sample used in that 
work and the main experimental findings. The turnstiles 
in these experiments suffered from photon-assisted tun- 
neling due to insufficient electromagnetic protection, see 



Section II. E and hence the quantization accuracy was 
only on 1CP^ level. Improved accuracy is expected for a 



new generation of turnstile devices (Pekola et al. 20101. 
For parallel turnstiles, a common bias voltage can be used 
as it is determined by the superconducting gap A, which 
is a material constant and varies only very little across a 
film deposited. Also, the rf drive can be common if the 
devices have roughly equal Rt, Eq, and coupling from 
the rf line to the island. As the error processes that set 
the ultimate limit on a single turnstile accuracy are not 
yet determined, the exact requirements on device unifor- 
mity cannot be fully resolved. 



H. Single-electron readout and error correction schemes 

I. Techniques for electrometry 

The electrometer used to detect the presence or ab- 
sence of individual charge quanta is a central component 
in schemes for assessing pumping errors and error correc- 
tion. Figure 33 (a) introduces the essential components 
of an electron counting setup. In order to observe proper- 
charge quantization, the counting island is connected to 
other conductors only via low-transparency tunnel con- 
tacts. The electrometer is capacitively coupled to the 
counting island and biased in such a manner that the 
small voltage drop of the counting island due to change 
of its charge state by one electron induces a measur- 
able change in the electrical transport through the de- 
tector. The readout performance can be characterized in 
terms of response time (bandwidth), charge sensitivity, 
and back-action to the system under measurement. In 
the present context of electron counting, we define back- 
action to include all mechanisms by which the presence 
of the detector changes the charge transport in the mea- 
sured system. 

The two basic electrometer realizations providing suf- 
ficient charge sensitivity for electron counting applica- 
tions are the single-electron transistor (SET) (Fulton and 



Dolan 1987 Kuzmin et al. 1989 and the quantum point 



contact (QPC) (|Berggren et al. 



1986; Field et al. 1993 



Thornton et al. 1986p . From a sample fabrication point 



of view, it is convenient when the electrometer and the 
charge pump can be defined in the same process; hence, 
QPC is the natural charge detector for quantum dots in 
semiconductor 2DEGs, whereas metallic single-electron 
devices are typically probed with SETs. There exist also 
studies where a metallic superconducting SET has been 



used as the electrometer for a semiconductor QD ( 


Fuji- 


sawa et al. 


2004 Lu et al. 


2003 


Yuan et al. 


2011|), and 


the SET can be realized in the 2DEG as well (Morello 



et al. 2010). 



The charge sensitivity Sq is determined by the noise 
of the system as a whole (Korotkov 1994) and is con- 
veniently expressed in units of e/\/Hz for electrometry 
applications. For metallic SETs, output voltage fluc- 
tuations 8V ou t can be related to the charge coupled to 
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FIG. 33 (a) Circuit diagram of a charge counting device. 
Electric charge Q on the island on the left is monitored. The 
island is coupled to an electrometer island via capacitor C x , 
and also tunnel coupled to an external conductor. The single- 
electron box configuration illustrated here requires only one 
tunnel junction with capacitance Cj. In addition, there is 
capacitance Co to ground, which accounts also for gate elec- 
trodes and any parasitic capacitances. The probing current 
7det through the detector is sensitive to the charge on the cou- 
pling capacitor, which is a fraction C x /Cs of the total charge 
Q, where Ce = C x +Co- The detector is a single-electron tun- 
neling transistor based on Coulomb blockade, and hence the 
total capacitance of the detector island Cdet is of the order of 
1 fF or less, (b) Circuit diagram of a general noisy electrical 
amplifier that can also be adapted to describe the electrom- 
eters of single-electron experiments. This figure is based on 
( Devoret and Schoelkopf| |2000). For the configuration shown 
in panel (a), one has for input impedance Z- m (iS) = l/(jwCm), 
where C^ 1 = CJ 1 + C^ t . The input voltage is related to the 
island charge Q through Vi n = Q/Ce- Noise source 7jv repre- 
sents back-action and Vn the noise added by the electrometer 
at the output referred to the input. Gain of the amplifier is 
given by G. Output impedance Z out equals the differential 
resistance at the amplifier operation point. 



the electrometer according to Sq = C g SV ou t/(dV out /dV g ), 
where V g is the voltage of the SET gate electrode and C Q 
is its capacitance to the SET island ( |Kuzmin et al. 1989 1 . 
Here, C g can be determined reliably in the experiment 



from the period of Coulomb oscillations. Similar cali- 
bration cannot be performed for a QPC and hence the 
charge sensitivity is expressed in relation to the charge of 
the neighboring QD (Cassidy et al. 20071, correspond- 
la) 



ing to Q in Fig. |33 

as Sq = k 



Variations of Q and q are related 
where k = C x /C% is the fraction of the 
island charge that is coupled to the electrometer. For 
charge counting applications, the relevant parameter is 
SQ. The RMS charge noise for a given detection band 

equals Aq 



' J duj SQ out (uj), which reduces to Sq x \IB 

in the white noise limit, where B is the readout band- 
width. It is possible to pose the charge detection prob- 
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FIG. 34 Detector back-action mechanisms. The back-action 
can originate by direct electromagnetic (EM) coupling either 
by variations in pump biasing or by high frequency photon 
assisted tunneling (PAT). Another source of back-action is 
via heat conduction. The detector located in proximity of the 
device typically heats up. The heat can be then conducted to 
the device by either phononic or photonic coupling. 



lem in the language of quantum linear amplifiers as shown 
in Fig. |33tb) (|Averin| 



2003 


Clerk et al. 


2010 


Devoret 



and Schoelkopf] |2000 ). When such a detector is mod- 
eled as a linear voltage amplifier, Jjv an d Vn chracterize 
the input and output noise, respectively, and the quan- 
tum theory limit for the spect ral density o f fluctuations 
at signal frequency u) reads \J Sv{w)Si(uj) > hco/2. In- 
formation about the electronic back-action is contained 
in the correlator (SVi n (t) SVi D (t')) of the induced voltage 
fluctuations on the counting island. Denoting the total 
capacitance of the counting island by Ce, the fluctuations 
in the output charge signal equal SQ out (u>) — CeVjv(w), 
and the voltage fluctuations on the counting island are 
given by 5Vi n (u)) — 7jv(w)/(jwCs). One thus finds 
\j SQout (^Sv^ (w) > h/2 as the quantum limit. 

In theory, quantum-limited operation can be achieved 
with normal-state SETs operated in the co-tunneling 



regime ( Averin 


20011, superconducting SETs ( 


Zorin 


1996 


2001D, and QPCs ( 


Averin and SukhorukovJ 2005 


Clerk et al. 2003 Korotkov 1999 ) . In practical devices 



however, the noise spectrum up to 1-100 kHz depending 
on the setup is dominated by 1 / /-like charge noise that 
is intrinsic to the sample but whose microscopic physi- 



cal origin is still debated ( 


Buehler et al 


2005 


Starmark 


et al.\ 


1999 


Vandersypen et al. 


2004 


. Above 1 kHz, 



the charge noise level is usually set by the preamplifier 
noise, but studies exist where the intrinsic shot noise of 
the electrometer was comparable to the noise of the read- 



2006 



Kafanov and Dels- 



out electro nics ( Brenning et al. 
ing 2009). For the normal-state SET, sensitivities of 
the order of 10 -7 e/y/Hz are attainable in theory with 
present-day fabrication technology, where the intrinsic 
noise is due to stochastic character of the tunneling pro- 
cesses and includes both shot and thermal noise (Ko- 
rotkov 1994 Korotkov and Paalanen 1999). The best 
charge sensitivities reported to date for a single-electron 



:-!S 



transistor in ( Brenning et al. 2006 ) were almost iden- 



tical in normal and superconducting states, namely 1.0 
and 0.9 x 10~ 6 e/ yTIz, respectively , at a signal frequency 
of 1.5 MHz. In (Xue et al. 2009), also the back-action 



of a superconducting SET was measured and the prod- 
uct of noise and back-action was found to be 3.6 times 
the quantum limit. For QD charge detection with QPCs, 
charge sensitivity of 2 x 10~ 4 e/yTlz referred to the QD 
charge has been demonstrated ( Cassidy et al. 2007 ) . It 



appears to be easier to realize large charge coupling frac- 



tion k with metallic SETs than with QPCs (Yuan et al. 

20TTI). 



Let us now discuss the back-action mechanisms in 
more detail. Despite the above quantum theory result 
connecting back-action and noise, the electronic back- 
action of electron counter can be addressed in princi- 
ple independently of its charge noise, as the readout 
bandwidth (at most 100 MHz, see below) is much be- 
low the microwave frequencies that can promote charge 
transfer errors: Overcoming even a modest 100 fiV en- 
ergy barrier requires photon frequencies above 24 GHz if 
multi-photon processes are neglected. Nevertheless, volt- 
age fluctuations induced by the shot noise of the detec- 
tor usually have a non-negligible spectral density at mi- 
crowave frequencies. For SET detectors carrying quasi- 
particle current, the voltage fluctuations of the island 
have a spectral density Sv Aeii {w) = 2eidet/(w 2 C d ! et ) in the 
limit uj 27r/dct/e, of which fraction k 2 is transferred to 
the counting island via capacitive voltage division. This 
can easily be stronger than the equilibrium thermal noise 



from resistive components at the sample stage; cf. (Mar- 



tinis and Nahum 1993). A full quantum calculation has 



been presented in (Johansson et al. 2002). The high- 



frequency back-action can be attenuated by replacing the 
capacitive coupling by a lossy wire that acts as a low- 
pass filter for microwaves but does not affect the charge 



signal (Saira et al. 



2012a). For QPCs, Coulombic back- 



action can be divided into shot noise, which can be in 
principle eliminated by circuit design, and fundamental 



and Clerk 2010 ) 



charge noise (Aguado and Kouwenhoven 2000 Young 



The low-frequency part of detector back-action man- 
ifests itself as variation of the dc bias of the pump or 
turnstile device. The case of an SET electrometer cou- 



pled to a single-electron box was studied in ( Turek et al. 
2005 ) . In the limit of small coupling capacitance C x , the 
voltage swing on the counting island due to loading and 
unloading the detector island equals AV[ n = ne/Cdct- We 
note that this is just a fraction C x /Cd c t < 1 of the volt- 
age swing from loading or unloading the actual counting 
island with an electron. Hence, dc backaction of the de- 
tector does not necessarily place an additional constraint 
on the design of the counting circuit, as the voltage swing 
e/Cs needs to be kept small for the operation of the 
charge pump. 

In addition to the electronic back-action described 
above, one needs to consider the phononic heat conduc- 
tion from the detector to the charge pump. For reach- 



ing the ultimate accuracy, the charge pumps typically 
require temperatures of the order of 100 mK or lower, 
where small on-chip dissipation can raise the local tem- 
perature significantly due to vanishing heat conductivity 
in the low temperature limit ( Giazotto et al. , 2006 ; Kautz 



et al. 



patec 



1993|), see also Sec. II. F| The average power dissi- 

^dctKiot) 



by the detector equals P = (IdetKiet), and it needs 
to be transported away by the substrate phonons or elec- 
tronically via the leads. Requirement for a sufficiently 
large charge coupling coefficient k limits the distance by 
which the detector and charge pump can be separated. 
The temperature increase by dissipated power has been 



studied on a silicon substate in (Savin et al. 
they give the formula 



2006) and 



T = 



T 4 
1 o 



2fP 



1/4 



(36) 



where T is the substrate temperature at distance r from 
a point source of heating power P, To is the bath temper- 
ature, vv = 3600 Wm" 2 K" 4 is material parameter and 
/ = 0.72 is a fitting parameter for their experimental 
observations. For an illustrative example, we estimate 
that the dissipated power at the electrometer in the orig- 



inal RF-SET paper (Schoelkopf et al. 1998) was 120 fW 



based on the published numbers. According to Eq. (36), 



this will heat the substrate underneath nearby junctions 
(r = 200 nm) to 140 mK, which is high enough to dete- 
riorate the performance of many single-electron devices 



below the metrological requirements. In ( Sillanpaa, et al. 
2004), the readout was coupled to the Josephson in- 



ductance of a superconducting SET instead of conduc- 
tance, reducing the dissipation by two orders of mag- 
nitude. Usually, it is possible to assess the severity of 
detector back-action effects in the experiment by mea- 
suring the tunneling rates using d i fferent values of I& e t, 
see, e.g., (|Kemppinen et al] |2011| |Lotkhov et al] [2011 



Saira et al. 2012a I, so that any variation of the observed 
rates can be attributed to back-action. The picture is 
somewhat different in 2DEG systems due to significantly 
weakened electron-phonon coupling. Experimental study 



of phononic back-action in 2DEGs is presented in ( Har- 



busch et al. 2010 Schinner et al. 2009) 



Finally, let us consider the obtainable bandwidth in 
single-electron detection. Digital low-pass filtering is ap- 
plied to the digitized output of the last amplifier stage 
to reduce the noise level so that some form of thresh- 
old detection can be applied for identifying the individ- 
ual tunneling events. The optimal filter, also called the 
matched receiver [see, e.g., P^aylp98] )], in the case of 
a known transfer function can be represented in Fourier 
space as F(u) = A*(oj)/ [\A(uj)\ 2 + S n {lo) 2 / a 2 }, where 
A(uS) is the (complex) gain of the readout circuit, Sn(^>) 
is the spectral density of noise at the filter input, and a 
is r.m.s. signal power. The matched receiver will mini- 
mize the expected squared error ([Q ut(t) — Q{t)] 2 ), but 
in practice any low pass filter with an appropriate cut- 
off will work. The cutoff may have to be set lower than 
the intrinsic bandwidth in order to suppress charge noise 
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level to a sufficiently small fraction of e. The bandwidth 
of the readout, B, is commonly defined as the corner fre- 
quency of the gain from gate charge to output voltage 
(Visscher et al. 1996). The performance requirements 



for the charge readout depend on the particular charge 
counting scheme, but in general the bandwidth B places 
a limit on the fastest processes that can be detected and 
hence constrains the magnitude of the electric current 
that can be reliably monitored. 

Both the QPC and SET electrometers have an 
impedance of the order of Rk = h/e 2 25.8 M7. For the 
SET, R > Rk is required to realize Coulomb blockade 
accordi ng to the orthodox theory of single-electron tun- 
neling (Averin and Likharev 1991| Ingold and Nazarov 
1992). For QPC, the most charge senstive operation 



point is around a bias point where dV/dl — Rk, mid- 
way between the first c onductance plateau and pinch- 
off ( Cassidy et al. 2007 ) . The readout bandwidth is set 



by the RC constant of the electrometer's differential re- 
sistance and the capacitive loading on its outputs. As 
the barrier capacitance is of the order of 1 fF or less 
for the devices, the intrinsic bandwidth is in the GHz 
range. In practice, the capacitance of the biasing leads 
and the input capacitance of the preamplifier dominate. 
When the preamplifier is located at room temperature as 



in the pioneering experiments (Fulton and Dolan 1987 



Kuzmin et al. 19891, the wiring necessarily contributes 



a capacitance of the order of 0.1-1 nF and henceforth 



limits the readout bandwidth to the kHz range ( Petters- 



son et al. 1996 Visscher et al. 1996). Readout by a 



current amplifier from a voltage-biased SET avoids the 
RC cutoff on the gain, but the usable bandwidth is not 
substantially altered as current noise increases at high 
frequencies where the cabling capacitance shorts the cur- 



rent amplifier input (Starmark et al. 1999) 



In order to increase the effective readout bandwidth, 
the SET impedance has to be transformed down towards 
the cable impedance that is of the order of 50 il. Band- 
widths up to 700 kHz have been achieved by utilizing a 
HEMT amplifier with a low impedance output at the 
samp le stage (Pettersson et al. |1996| Visscher et al. 



1996 ). The dissipated power at the HEMT in these stud- 



ies was 1-10 fxW depending on the biasing, which can eas- 
ily result in overheating of the electrometer and/or the 
coupled single-electron device. The best readout configu- 
ration to date is the rf reflectometry technique, applicable 



to both SETs (Schoelkopf et all 1998) and QPCs (Qin 



and Williams 2006), where the electrometer is embedded 



in a radio frequency resonant circuit and the readout is 
achieved by measuring the damping of the resonator. A 
read-out bandwidth of 100 MHz was achieved in the orig- 
inal demonstration ( Schoelkopf et al. 1998 ). The authors 



also note that their charge sensitivity of 1 x 10~ 5 e/vHz 
yields Aq = O.le for the full detection bandwidth, i.e., 
electron counting at 100 MHz would have been possible 
in a scenario where the charge coupling fraction k was 
close to unity. 



2. Electron-counting schemes 

Realization of a current standard based on electron 
counting has been one of the key motivators for devel- 



opment of ultrasensitive electrometry ( Gustavsson et al. 



2008 Keller 2009 Schoelkopf et al. 1998). First, let 



us see why direct current measurement of uncorrelated 
tunneling events, like those produced by a voltage-biased 
tunnel junction, cannot be used for a high-precision cur- 
rent standard: Assume a noise-free charge detector that 
yields the charge state of the counting island with time 
resolution t = 1/B, and that Markovian (uncorrelated) 
tunneling events occur at rate F -C B. With probabil- 
ity Ft, a single tunneling event occurs during the time 
r and is correctly counted by the detector. With proba- 
bility (rr) 2 /2, two tunneling events occur within r and 
constitute a counting error. Hence, to achieve a relative 
error rate p, one needs T < 2pB. Even with a noiselees 
100 MHz rf-SET, one could not measure a direct cur- 
rent greater than 2e/s at metrological accuracy p = 10~ 8 
in this manner. Would it be practical to account sta- 
tistically for the missed events in a manner similar to 
( Naaman and Aumentado 2006 ) assuming truly Poisso- 



nian tunneling statistics and a well-characterized detec- 
tor? The answer is unfortunately negative: If N tunnel- 
ing events are observed, the number of missed events M 
is a Poissonian variable with a mean of NTt/2 and stan- 
dard deviation 8M = 0VTV/2. Requiring SM < pN 
gives N > r-r/(2p 2 ). For F = 1 MHz and r, p as above, 
one has to average over N > 5 x 10 13 events, which is im- 
practical. A more deta iled calculation based on Bayesian 
inference presented in ( Gustavsson et al , 2009 ) results in 



the same N dependence. 

Charge transport through a ID-array of tunnel junc- 
tions can take place in a form of solitons depending on 
device parameters ( Likharev 1988 Likharev et al. 19891. 



Propagation of the solitons promotes time correlation 
in the electron tunneling events, allowing the accuracy 
limitations of counting uncorrelated electrons presented 
above to be lifted. A proof-of-concept experimental re- 
alization has been presented in ( Bylander et al. 2005 1 . 



The array is terminated at the middle island of an SET, 
allowing for unity charge coupling, and a signal centered 
around frequency f c — I/e is expected. The authors 
claim a possible accuracy of 10 -6 based on the charge 
sensitivity of their electrometer only. However, the spec- 
tral peaks in the experimental data appear too wide for 
an accurate determination of the center frequency. Fac- 
tors not included in the accuracy estimate are instability 
of the bias current and SET background charge fluctua- 
tions. 

Single-electron electrometry can be used to count 
the much rarer pumping errors instead of the total 
pumping current. Such an approach has been used to 
study the accuracy of metallic multi-junction pumps that 
are used in the electron counting capacitance standard 
(ECCS) (|Keller et al] p99l [2007]. A circuit diagram 



of an ECCS experiment is shown in Fig. 35 Two cryo- 
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FIG. 35 Circuit diagram of a practical implementation of 
the electron counting capacitance standard. Switches NS1 



and NS2 are cryogen ic needle switches. Figure from Ref. (Ca- 
marota et al. 2012 1. 



none have been presented to date. Measurements of two 
series-connected semiconductor QD pumps with a QPC 
charge detector coupled to an island in the middle are 
presented in (Fricke et al. . 2011 ), although observation of 



quantized pumping errors was not demonstrated there. 



I. Device fabrication 

Fabrication of charge pumps, regardless of their opera- 
tional principle, requires advanced nanofabrication meth- 
ods. These include, for example, electron-beam lithogra- 
phy, various dry etching techniques, and MBE growth 
of semiconductor heterostructures. In general, pumping 
devices can have small feature sizes in multiple layers 
that must be aligned with each other accurately. We be- 
gin with the description of the fabrication procedure for 
the metallic single-electron and Cooper pair pumps and 
turnstiles described in Sees. |IH.A[ \lTLB\ and[TTLE] Sub- 
sequently, we present the fabrication methods of quan- 
tum dot pumps and turnstiles, the operation of which is 
discussed in Sec. MI. CI 



genie needle switches are required to operate the device 
in different modes: determining coupling capacitances 
and tuning the pump drive signal (NS1 and NS2 closed), 
operating the pump to charge C cryo (NS1 closed, NS2 
open), and comparing C cryo with an external traceable 
capacitor (NS1 open, NS2 closed). With NS1 open, the 
pump can be operated in a shuttle mode: a charge of e 
is repeatedly pumped back and forth across the pump at 
the optimal operation frequency (which is above detec- 
tor bandwidth), and pumping errors appear as discrete 
jumps in the electrometer output. A relative error rate of 
1.5 x 10~ 8 at a shuttling frequency of 5.05 MHz has been 
demonstrated for a 7-junction pump shown in Fig. |12 



which was used in the NIST ECCS setup (Keller et al. 
|1996| ). AtPTB, accuracy of the order of 10 was demon- 
strated for a 5-j unctio n R-pump operated at a shuttling 
frequency of 0.5 MHz (Camarota et al. 2012). In the ac- 



tual ECCS operation, the SET electrometer is used as a 
part of a feedback loop that maintains the voltage of the 
node at the end of the pump constant, but time-resolved 
single-electron detection is not needed then. 

(Wulf and Zorin 2008) propose various strategies 



for error accounting utilizing modestly accurate charge 
pumps and single-electron electrometry. They consider 
the ECCS, and in addition propose schemes involving 
several charge pumps connected in series. The advantage 
of the more complex schemes is that they do not require 
cryogenic needle switches, and in theory allow one to ac- 
count for the sign of pumping errors, thus improving the 
accuracy of the setup beyond the intrinsic accuracy of 
the pumps. Even though the accuracy requirements on 
the pumps are relaxed somewhat, the added complexity 
from operating several pumps and electrometers make 
the practical realization of these ideas challenging, and 



1. Metallic devices 

Metallic single-electron and Cooper pair pumps and 
turnstiles are typically made by the angle deposi- 
tion technique, which was first introduced by G. J. 
Dolan Dolan (19771] for the photolithography process 
and then later adapted by Dolan and Dunsmuir for the 



electron-beam lithography process Dolan and Dunsmuir 



( 1988 )]. We note that there is a myriad of different ways 
fabricating these devices. Below, we describe only a cer- 
tain fabrication process for these devices in great detail 
instead of giving a thorough study of all possible varia- 
tions. 

The process starts with the deposition of Au layer on 
an Si wafer covered by a native silicon oxide. The Au 
pattern is formed by a standard photolithography and 
liftoff process using photoresist S1813 and contains con- 
tact pads and on-chip wiring as well as alignment mark- 
ers for the deposition of the subsequent layers. Next, a 
tri-layer resist structure is built (from bottom to top): 
copolymer/Ge/poly-methyl-methacrylate (PMMA) with 
the thicknesses 200 nm, 20 nm and 50 nm, respectively 
(See Fig. [36^a). The polymer layers are spin-coated on 
the wafer and baked in a nitrogen oven, and the Ge layer 
is deposited in an electron gun evaporator. The wafer 
is then cleaved into smaller pieces which are exposed 
and processed separately. After the exposure of the top 
PMMA layer on one of the pieces in the electron beam 
writer, e. g., JEOL JBX-5FE, the piece is developed at 
room temperature in isopropyl alcohol (IPA) mixed with 
methyl isobutyl ketone (MIBK) at a ratio of 3:1. Thus, 
a desired pattern is formed in the PMMA layer (Fig. 
36jb)). The pattern is transferred into the Ge layer by 
reactive ion etching in CF4 (Fig. [36[c)). The sample 
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is then placed in an electron cyclotron resonance (ECR) 
etcher, in which an undercut is formed by oxygen plasma. 
The undercut depth is controlled by the tilt of the sam- 
ple stage in the ECR machine. At the same time, the 
top PMMA layer is etched away. At this stage, each chip 
has a Ge mask supported by the copolymer layer (Fig. 
36|(d)) . Some parts of the mask are suspended form- 
ing the Dolan bridges. Although we described above a 
method with three layers, in many cases a bi-layer mask 
composed of copolymer and PMMA resists is sufficient. 

The chips with masks are placed in an electron gun 
evaporator equipped with a tilting stage. Two consec- 
utive depositions of metal through the same mask are 
carried out at different angles to create a partial overlap 
between the metal layers (Fig. |36|e)). If the surface of 
the bottom layer (typically Al) is oxidized by introducing 
oxygen into the evaporation chamber, after the deposi- 
tion of the top electrode, the sandwich structure compos- 
ing of the overlapping metal layers with a thin oxide in 
between form small tunnel junctions (Fig. |36[ f)). 

The normal-metal or superconducting charge pumps 
are made entirely of Al, which can be turned normal 



bottom electrode was Al and t 
or AuPd. 



Geerligs et al. 
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FIG. 36 Fabrication of metallic devices, (a) Build-up of a 
tri-layer resist structure and exposure in the electron-beam 
writer; (b) Development of the top PMMA layer; (c) Transfer 
of the pattern formed in the resist into the Ge layer by reactive 
ion etching; (d) Creation of undercut in the bottom resist 
and removal of the top resist by oxygen plasma; (e) angle 
deposition of metals with an oxidation in between; (f) the 
resulting structure after the liftoff process. 



2. Quantum dots 

The gate structure of the charge pumps based on quan- 
tum dots is also fabricated using electron beam lithog- 
raphy. The main differences in the fabrication com- 
pared with metallic devices are the following: Ohmic con- 
tacts have to be made between metallic bonding pads on 
the surface of the chip and the 2DEG located typically 
^100 nm below the surface. Furthermore, the 2DEG 
has to be either depleted with negative gate voltage from 
the unwanted positions in the case of GaAs devices [see 
Fig. [17] (d)] or accumulated with positive gate voltage 
in the case of metal-oxide-semiconductor (MOS) silicon 
devices [see Fig. 17 (c)]. For GaAs, also etching tech- 
niques have been employed to dispose some parts of the 
2DEG leading to a smaller number of required gates [see 
Fig. [19] (a)]. In GaAs devices, typically a single depo- 
sition of metal through a mono-layer PMMA resist is 
sufficient to create the gate structure. For MOS sili- 
con dots, several aligned layers of gate material are often 
used. However, only a single layer is typically deposited 
with each mask in contrast to metallic devices employing 
angle evaporation. 

Let us describe in detail a fabrication process for MOS 
silicon quantum dots. We begin with a high-resistivity 
(p > 10 kSlcm at 300 K) near-intrinsic silicon wafer. 
Phosphorous atoms are deposited on the silicon surface 



using standard photolithography and they diffuse to a 
depth of roughly 1.5 /mi during the growth of a 200 nm 
field silicon oxide on top. All the following process steps 
involving etching or deposition have to be aligned with 
the previous ones with the help of alignment markers, a 
routine we do not discuss separately below. Then a win- 
dow with size 30 x 30 pm 2 is opened to the field oxide 
and replaced by 8-nm-thick high-quality Si02 gate oxide 
that is grown in ultradry oxidation furnace at 800 °C in 
O2 and dichlorocthylene. This thin-oxide window over- 
laps by a few micrometers with the ends of the metallic 
phosphorous-rich n + regions. The field oxide is etched se- 
lectively above the other ends of the n + regions formed in 
the previous process. The ohmic contacts and the bond- 
ing pads are made by depositing metal on these etched 
regions forming a connection to the n + silicon. Subse- 
quent annealing is employed to avoid the formation of 
Schottky barriers. 

At this stage, we have bonding pads connected to the 
metallic n + regions that extend some hundred microme- 
ters away from the pads to the thin-oxide window with 
the line width of four micrometers. Electron beam lithog- 
raphy with a 200-nm PMMA resist and metal evapora- 
tion with an electron gun evaporator is employed to de- 
posit the first layer of aluminum gates inside the window 
and their bonding pads outside the window. After the 
lift-off, the gates are passivated by an Al^Oj, layer formed 
by oxidizing the aluminum gates either by oxygen plasma 
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or thermally on a hot plate (150 °C, 5 min). The oxide 
layer electrically insulates completely the following over- 
lapping layers of aluminum gates that are deposited in 
the same way with alignment accuracy of ~20 nm. 

At least one gate has to overlap with areas where n + 
regions extend to the thin-oxide window. By applying 
positive voltage on these reservoir gates, the electrons 
from the n + are attracted to the Si/SiC>2 interface below 
the reservoir gates forming the source and drain reser- 
voirs of the device. For example, the device shown in 



in quantized voltage plateaus 



Fig. 22 (b) is composed of one or two layers of gates: 
one top gate that induces the source and drain reservoir, 
two barrier gates below the top gate defining the quan- 
tum dot, and a plunger gate in the same layer with the 
barrier gates. Finally, a forming gas (95%N2, 5%H2) an- 
neal is carried out for the sample at 400 °C for 15 min 
to reduce the Si/SiC>2 interface trap density to a level 
~5xl0 10 cm _2 eV~ 1 near the conduction band edge. Sil- 
icon quantum dots can also be fabricated with all-silicon 
process, in which the aluminum gates are replaced by 
conducting polysilicon gates shown in Fig. 22 (a). 



IV. QUANTUM STANDARDS OF ELECTRIC 
QUANTITIES AND THE QUANTUM METROLOGY 
TRIANGLE 



The ampere is one 
of The International 



of the seven 
System of 



base units 
Units (SI) 



(Bureau International des Poids et Mesures 



2006) 



and is defined as follows: "The ampere is that constant 
current which, if maintained in two straight parallel 
conductors of infinite length, of negligible circular 
cross-section, and placed 1 m apart in vacuum, would 
produce between these conductors a force equal to 
2 x 10 7 newton per metre of length." The present 
definition is problematic for several reasons: (i) The 
experiments required for its realization are beyond the 
resources of most of the National Metrology Institutes, 
(ii) The lowest demonst rated uncertainties are not better 
than about 3 x 10 



(Clothier et al. 



1989 



Funck and 



Sienknecht 19911. (iii) The definition involves the unit 
of newton, kg x m/s 2 , and thus the prototype of the 
kilogram, which is shown to drift in time (Quinn 1991 1. 



In practice, electric metrologists are working outside the 
SI and employing quantum standards of voltage and 
resistance, based on the Josephson and quantum Hall 
effects, respectively. 



A. The conventional system of electric units 

According to the ac Josephson effect, V = (h/2e) x 
(d(j)/dt), the voltage V applied over the Josephson junc- 
tion induces oscillations of the phase difference <j> between 
the macroscopic wavefunctions of the superconductin g 
electrodes on each side of the junction (Josephson 1962 ). 
Phase-locking <j> by a high-frequency (/j) signal results 



Vr 



h njfj 



(37) 



which are often called Shapiro steps (Shapiro 19631. 



Here, Kj is called Josephson constant and nj is the in- 
teger number of cycles of 2ir that <fi evolves during one 
period of the high-frequency signal. 

The Josephson voltage standards ( JVS) have been used 
in electric metrology since the 1970s, see, e.g., |Jeanneret| 



and Benz| ( |2009[ ); |Kohlmann et al.\ ( |2003| ) for reviews 
The first standards consisted of a single junction and 
generated voltages only up to about 10 mV. Arrays of 
more than 10,000 junctions with the maximum output of 
10 V were developed in the 1980s. They were based on 
hysteretic junctions where Shapiro steps with different 
nj can exist at the same bias current. Since 1990s, the 
research has focused on arrays of nonhysteretic junctions 
where nj can be chosen by the applied current bias. Ar- 
rays divided in sections of 2 m junctions (m = 0,1,2...) 
are called programmable since one can digitally select any 
multiple of fj/Kj up to the number of junctions as the 
output voltage ( |Hamilton et al. 1995 Kohlmann et al. 
2007). They are practical for dc voltage metrology, but 



are especially developed for generating digitized ac volt- 
age waveforms up to about 1 kHz, which is an active 
research topic ( Behr et al. 2005 1 . Voltage waveforms 



at higher frequencies can be generated by pulse-driven 
Josephson junction arrays where the desired ac wave- 
form is synthesized by the delta-sigma modulation of fast 
voltage pulses, each having the time integral of one flux 
quantum h/(2e) (Benz and Hamilton 1996). 



The quantum Hall resistance standard (QHR) consists 
of a two-dimensional electron gas, which, when placed 
in a high perpendicular magnetic field, exhibits plateaus 
in the Hall voltage Vh — RhI over the sample in the 
direction perpendicular to both the field and the bias 
current /. Here, 



Rh = 



1 h 

i K e 2 



Rk 

IK 



(38) 



is the quantized resistance which is proportional to the 
von Klitzing constant Rk and inversely proportional to 
the integer ik (von Klitzing et al. 1980). The plateau 



index ix can be chosen by tuning the magnetic field. 
Usually, the best results are obtained at ix = 2. 

Quantum Hall standards based on Si MOSFET or 
GaAs/AlGaAs heterostructures were taken to routine use 
quickly during the 1980s, see, e.g. , IJeckelmann and Jean- 



neret| (|2001[); |Poirier and Schopfer| ( |2009[ ) ; |Weis and von 
Klitzingj ( |2011[ ) and Issue 4 of C. R. Physique, vol. 369 
(2011) for reviews. Different resistances can be calibrated 
against the QHR by using the cryogenic current com- 
parator (CCC). It is essentially a transformer with an 
exact transform ratio due to the Mcissner effect of the 



superconducting loop around t he windings (Gallop and 
Piquemal 2006 Harvey 1972). Another way to divide 
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or multiply Rh are parallel and series quantum Hall ar- 
rays, respectively, which are permitted by the technique 
of mult iple connections that suppresses the contact resis- 
tances ( |Delahaye]|1993[ ). One rapidly developing research 
topic are ac quantum Hall techniques, which can be used 
in impedance standards to expand the traceability to ca- 



pacitance and inductance (Schurr et al. 2011). An im 



portant recent discovery is that graphene can be used 
to realize an accurate and very robust QHR standard 
( Janssen et al. 2011[ Novoselov et al.\ |2007| |Tzalenchuk 
et aZ.||2010| |Zhang et al.\\200ty . 



The most precise measurement of Kj within the SI 
was performed by a device called liqui d-mercury elec- 
trometer with the uncertainty 2.7 x 1CP 7 (Clothier et al. 



[19891 ■ The SI value of r k 

can by obtained by comparing 

the impedance of the QHR and that of the Thompson 



Lampard calculable capacitor (Bachmair 2009 Thomp 



son and Lampard 19561. The lowest re ported uncer 



tainty of such comparison is 2.4 x 10 8 (Jeffery et al. 



1997). However, both JVS and QHR are much more 
reproducible than their uncertainties in the SI, see 
Sec. |IV.B| Therefore the consistency of electric measure- 
ments could be improved by defining conventional values 
for Rk and Kj. Based on the best available data by 
June 1988, the member states of the Metre Convention 
made an agreement of the values that came into effect in 
1990: 



Rk- 



00 

on 



483597.9 GHz/V 
25812.807 n 



(39) 



Since then, electric measurements have in practice been 
performed using this conventional system which is some- 
times emphasized by denoting the units by Vgo, O go , A 90 
etc., and where JVS and QHR are called representations 
of the units. 



B. Universality and exactness of electric quantum 
standards 

A theory can never be proven by theory, but, as argued 
by Gallop ( 2005 ) , theories based on very general princi- 



ples such as thermodynamics and gauge invariance are 
more convincing than microscopic theories such as the 



original derivation of the Josephson effect (Josephson 



1962). There are rather strong theoretical arguments for 
the exactness of JVS: Bloch has shown that if a Josephson 
junction is placed in a superconducting ring, the exact- 



ness of Kj can be derived from gauge invariance ( Bloch 



1968 1970). Furthermore, Fulton showed that a depen- 
dence of Kj on materials would violate Faraday's law 
(Fulton 1973). For quantum Hall devices, early theo- 



retical works argued that the exactness of Rk is a con- 



sequence of gauge invariance (Laughlin 



19811 Thouless 



et al. 1982). However, it is very complicated to model 
real quantum Hall bars, including dissipation, interac- 
tions etc. and thus the universality and exactness of 
QHR has sometimes been described as a " continuing sur- 
prise" flKellerl 120081 iMohr and Taylorl 120051) . Extensive 



theoretical work, e.g., on topological Chern numbers, has 
strengthened the confidence on the exactness of Rk, see 
Avron et al] ( |2003| ); |Bieri and Fr5hlich| ( |2011[ ); |Doucot 



( |201ip for introductory reviews. Recent quantum electro- 
dynamics (QED) based theoretical work predicts that the 
vacuum polarizatio n can lead to a mag netic field depen- 
dence of both R K (IPeninl 120091 [2010a| and Kj (IPenin 



2010bj), but only at the level of 10- zu . The case of single- 
electron transport has been studied much less and there 
are no such strong theoretical arguments for the lack of 



any corrections for the transported charge (Gallop 2005 



Keller 



2008 



Stock and Witt 2006 1. 1 



On the experimental side, comparisons between Si and 
GaAs quantum Hall bars show no deviations at the exper - 
imental uncertainty of ~ 3x 10 -10 (Hartland et al. 1991 1. 



Recently, an agreement at the uncertainty of 8.6 x 10 



was found between graphene and GaAs devices (Janssen 



et al. 2011). This is an extremely important demon- 
stration of the universality of Rk because the physics of 
the charge carriers is notably different in graphene and 
semiconductors (Goerbig 2011). Comparisons between 



JVSs have been summarized recently in | Wood and Solve] 
( 2009 ). The lowest uncertainties obtained in comparisons 
between two JVSs are in the range 10 -11 . Even much 
smaller uncertainties have been obtained in universality 
tests of the frequency-to- voltage conversion by applying 
the same frequency to two different junctions or junction 
arrays and detecting the voltage difference by a SQUID- 
based null detector. Several accurate experiments have 
indicated that the conversion is invariant of, e.g., the su- 
perconducting material and the junction geometry. The 
lowest demonstrated uncertainty is astonishing: 3x 10 _1 ' 



(Clarke 


1968 


Jain et al. 


1987 


Kautz and Lloyd 


1987 


Tsai et al. 


1983 


I- 



The reproducibility and universality of the quantum 
standards are an indication that Eqs. 



(|37j) and (|38j are 

exact, but a proof can only be obtained by comparisons 
to other standards. Any one of the electric quantities 
V, I or R can be tested by comparison to the other 
two via Ohm's law. If all three are quantum standards, 
such experiment is called quantum metrology triangle 



(QMT) (Likharev and Zorin 19851. It is a major goal 
in metrology, but the insufficient performance of single- 
electron devices has to date prevented the reaching of 
low uncertainties, see Sec. |IVD| However, the exactness 
of Eqs. ( 37 1 and ( 38 1 can also be studied in the framework 
of the adjustment of fundamental constants. The most 
thorough treatment has been performed by the Com- 
mittee on Data for Science and Technology (CODATA). 
Up date papers are nowadays publis h ed every four years, 
see IMohr and Taylorl pOOfjl [2005|; IMohr et al.\ (T20081 



1 A condensed-matter correction of ~ 10 10 e for the charge of 
the electron was suggested by theory based on QED l|Nordtvedt 
|1970|), but it was refuted in |Hartle et aL| (l971 1; Langcnb erg and] 
Schrieffer|([T97Ti. 
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FIG. 37 Simplified sketch of the most accurate routes to in- 
formation on ej,K- Direct measurement of Rk together with 
an independent measurement of the fine structure constant 
(a) yields a value on ek- Values for the sum of ej and ck can 
be obtained from the combination of the so called watt bal- 
ance experiment (KjRk) and a measurement of the Avogadro 
constant (Na), or from the combination a and measurements 
of low-field gyromagnetic ratios (710). Less accurate informa- 
tion is provided by measurements of high-field gyromagnetic 
ratios, Kj, and the Faraday constant F = eNA, and by the 
QMT. 



2012 ) 2 . Karshenboim (2009) provides a useful overview. 
We review here the most accurate (< 10~ 7 ) routes to in- 
formation on the electric quantum standards. They are 



also illustrated in Fig. 37 Most of the equations in this 
chapter assume that Eqs. (37 1 and (38) are exact, but 



when referring to possible deviations, we describe them 
by symbols er K ,s- 



R K = (l + cjf)4 



K,, = (1 
(1 
(1 



(40) 



es)e. 



In this context, the current generated by the single- 
electron current source is Ig =< kg > Qgf where 

< kg > is the average number of electrons transported 
per cycle. Thus any determination of eg will require both 
a measurement of the macroscopic current and counting 
the number of transported charges. 

There are a number of fundamental constants that 
are known by much smaller uncertainties than those re- 
lated to electric metrology. Some constants, e.g., the 
permeability, permittivity, and the speed of light in 
vacuum, and the molar mass constant, fj,o, cq, c, and 
M u = 1 g/mol, respectively, are fixed by the present 
SI. Examples of constants known with an uncertainty 

< 10~ 10 are the Rydberg constant and several rel- 
ative atomic masses, e.g., that of an electron, A r (e). In 
the past few years, there has been tremendous progress in 
the determination of the fine structure constant a. First, 
the electron magnetic moment anomaly, a e , was mea- 
sured with high accuracy. A separate calculation based 
on QED gives the function a(a e ). Together these results 



2 The adjustments are named after the deadline for the included 
data, e.g., CODATA-10 is based on experimental and theoretical 
results that were available by December 31, 2010. The values 
of the constants and much more information is available at the 
Web site physics.nist.gov/constants/. 



yield a value on a with an uncertainty 0.37 x 10 9 (Han 



neke et all 2008 ). Soon after, a measurement of the recoil 



velocity of the rubidium atom, when it absorbs a photon, 
yielded a value for a with an uncertainty of 0.66 x 10~ 9 
(Bouchendira et al. 2011 ). These two results are in good 



agreement. Together they give a validity check for QED 
since the first result is completely dependent and the lat- 
ter practically independent of it. 

The fine structure constant is related to Rk by the 
exact constants //q arL d c: 



2h 



2R 



K 



(41) 



This relationship means that when Rk is measured with 
a calculable capacitor, it also yields an estimate for a. 
Thus, a measurement of Rk could also test QED, but 
in practice, the atomic recoil measurement is more accu- 
rate by about factor 30. A metrologically more impor- 
tant interpretation of this relation is that a comparison 
between a and the weighted mean of the measure ments 
of Rk yields an estimate of €k = (26 ± 18) x 10~ 9 (Mohr 



et al. 2012 1. There is thus no proof of a nonzero ck > but 



several groups are developing calculable capacitors in or- 
der to determine Rk with the uncertainty below 10~ 8 



(Poirier and Schopfer 2009 Poirier et al. 2011) 



The existing data that yields information on ej is 
more discrepant. It turns out that ej is related to mea- 



surements of gyromagnetic ratios (see Mohr and Taylor 



(2000) for a deta iled descripti o n) and to th e realization 



of the kilogram (Mohr et al. 2008 2012). The gyro- 



magnetic ratio 7 determines the spin-flip frequency / of 
a free particle when it is placed in a magnetic field B: 
7 = 2tt f/B. The gyromagnetic ratios of a helium nu- 
cleus and a proton are accessible in nuclear and atomic 
magnetic resonance experiments. These ratios can be 
related to the gyromagnetic ratio of an electron that is 
linked to a and h. There are two methods to produce 
the magnetic field: In the low-field method, it is gener- 
ated by electric current in a coil and determined from the 
current and the geometry. In the high-field method, the 
field is generated by a permanent magnet and measured 
from a current induced in a coil. When the electric cur- 
rent is dete rmined in t erms JVS and QHR, the product 
KjRk = \/2/iQc/(ha) will appear either in the numera- 
tor or the denominator of 7, depending on which method 
is used. It turns out that in the low-field method, h can- 
cels out from the equations of the gyromagnetic ratios, 
and the experiment yields a value for a. The high-field 
results also depend on a, but since it is known much more 
precisely than h, they essentially yield a value for h. The 
high-field results on h are in good agreement with other 
experiments, albeit their uncertainty is not better than 
about 10 -6 . However, the low-field results are discrepant 
from the CODATA value of a. By substituting Eqs. (|40l 



to the observational equations of the low-field data, one 
obtains the estimate 6k + ej = (—25 4 ± 93) x 10~ 9 
(ICadoret et al.\ 120111 IMohr et al.\ 120121) . Since the mea- 



surements of Rk yield much smaller value for ck, the 
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gyromagnetic data seem to imply a significant negative 
ej. However, as explained below a positive ej can be 
found from measurements aiming at the redefinition of 
the kilogram. 

There are essentially two candidate methods for the 
future realization of the kilogram: the watt balance and 
silicon sphere methods. The first, suggested in |Kibble| 
(1975), relates electric power to the time derivative of 



£j,k can be found from the COD AT A papers IMohr et al. 



potential energy: 



mgv 



V 



R KjR K 



= h/A. 



(42) 



When the mass m, its velocity v and the gravitational 
acceleration g are traceable to the SI, the watt bal- 
ance yields a value on h. Watt balance results have al- 
ready been published by NPL, NIST, METAS, and NRC, 
and several devices are under devel opment. See Steele 
et al.\ (|2012|) for the latest status and |Eichenberger et al. 



( 2009[ ); |Li eFaT\ ( [2012] ); |Stock| ( |201l| for reviews. The 
silicon sphere approach is so demanding that it is em- 
ployed only by the international Avogadro consortium 
(IAC). The res ults were published in 2011, see Andreas 
et al. ( 2011a|b ) and the whole no. 2 of Metrologia, vol. 48 
(2011). This project determines the Avogadro constant 
Na by fabricating spheres of enriched 28 Si whose mass 
is compared to the prototype of the kilogram and whose 
volume is measured by laser interferometry. The lattice 
parameter and the relative atomic mass of 28 Si are mea- 
sured in different experiments, and the ratio of the rela- 
tive and absolute mass densities yields Na- 

Results on h and Na can be compared precisely with 
the help of the molar Planck constant 



N A h 



,A r (e)M u c 

2i?oo 



Its uncertainty is only 0.7 x 10~ 9 ( |Mohr et aL||2012| ) and 
depends mainly on those of A r (e) and a. This equa- 
tion can be derived from the definition of the Rydberg 
constant by writing the inaccurate absolute mass of the 
electron in terms of its relative mass and Na which links 
microscopic and macroscopic masses. The IAC 2011 re- 
sult resolved the discrepancy of 1.2 x 10 -6 between watt 
balances and the Avogadro constant determined from a 
sphere of natural Si that had puzzled metrologists since 
1998 flMohr and Taylor[ |2000D . E specially after the new 
NRC results (Steele et al. 2012) there is no longer any 



clear discrepancy between the two methods, but the two 
most accurate watt balances deviate by factor 260 x 10 -9 
which is 3.5 times the uncertainty of their difference 
dSteele et al.\ |2012| |Steiner et al.\ |2007|). Also the mea- 



surements of the isot ope ratio of the sili con sphere deviate 
more than expected (Yang et al. 2012). Nevertheless, by 



combining the Planck constants obtained from watt bal- 
ance and silicon sphere experiments, h w and ft-Avo, one 
can obtain an estimate of ej + ex/2 ~ (h,Avo/h w — 1)/2 = 
(77±18) x 10 -9 . Here, we neglected correlations between 
experiments. More detailed analysis on the existence of 



(2008 


2012 


), see also 


Keller 


(2008) 



least-squares analysis of fundamental constants several 
times allowing cither nonzero ex or ej, and including 
only part of the data. When they excluded the most 
accurate but discrepant data, the remaining less precise 
but consistent data yielded the conservative estimates 
e K = (28±18) x 10~ 9 and ej = (150±490) x 10" 9 . Thus 
the exactness of the quantum Hall effect is confirmed 
much better than that of the Josephson effect. 



C. The future SI 

Modernizing the SI towards a system based on funda- 
mental constants or other true invariants of nature has 
long been a major goal, tracing back to a proposal by 
Maxwell in the 19th century, see, e.g., Flowers (2004) 
and references therein. Atomic clocks and laser interfer- 



ometry permitted such a revision of the second and of 
the meter. The development of quantum electric stan- 
dards, watt balance experiments, the Avogadro project, 
and measurements of the Boltzmann constant have made 
the reform of the ampere, kilogram, mole, and kelvin re- 



alistic in the near future. Specially, suggestions in Mills 



et al. 


( 2005 ) launched an active debate among metrolo- 


gists 


Becker et al. 


2007 


Mills et al. 


2006 


Milton et al. 



tered before there are at least three independent experi- 
ments (both from watt balance and Avogadro constant) 
with uncertainties < 50 x 10 -9 that are consistent within 
the 95% confidence intervals, and at least one of them 
has the uncertainty < 20 x 10~ 9 (Glaser et al. 2010). 



There have also been requests to await better results from 



(43) 

single-electron and QMT experiments (iBordel 20051 |Mil- 



ton et al. 2007 1 , and to solve the discrepancy of low-field 



gyromagnetic experiments (Cadoret et al. 2011) 



There is already a draft chapter for the SI brochure 
that would adopt the new definitions: BIPM (2010), see 
also the whole issue 1953 in Phil. Trans. Royal Soc. A 
369 (2011), especially [Mills et al.\ fl2011| ). In this draft, 
the whole system of units is scaled by a single sentence 
that fixes seven constants. The most substantial changes 
are that the base units ampere, kilogram, mole, and 
kelvin are defined by fixed values of e, h, Na, and fcs, 
respectively. The new definition for the ampere would 
read: "The ampere, A, is the unit of electric current; its 
magnitude is set by fixing the numerical value of the ele- 
mentary charge to be equal to exactly 1.60217X x 10~ 19 
when it is expressed in the unit s A, which is equal to C." 
The new definitions will not imply any particular meth- 
ods for the realizations of the units. They will be guided 
by raises en pratique, e.g., the a mpere could b e realized 
with the help of JVS and QHR (ICCEMl [20121. 



The new SI would significantly lower the uncertain- 
tiesof many fundamental constants, see e.g., |Mills et al.\ 
(2011) for evaluations. One should note, however, that 



choosing the optimal set of fixed constants is always a 



4G 



tradeoff. For example, since a is a dimensionless number 
and thus independent of the choice of units, one can see 



from Eq. ( 41 ) that fixing e and h would make /j,q (and 
eo) a quantity that is determined by a measurement of a. 
Presently, and eo are fixed by the definition of ampere. 
However, their uncertainty woul d be very low, the same 
as that of a, which is 3.2 xltT 10 (jMohr et al. \\2012h. One 



alternative suggestion is to fix h and the Planck charge 
q p = v^ep/ic, w hich would keep fio and eo exact (Stock 
and Witt 2006 1 . It is also worth noting that out of h, Na, 
and the molar mass of carbon 12, M( 12 C) = A r ( 12 C)M u , 
only two can be fixed. The suggested SI would release 
the equality M( 12 C) = 0.012 kg/mol, which has raised 
criticism. Especially, there have been claims that the def- 
inition of the kilogram based on h would not be under- 
standable for the wider audience, and a definition based 
on the mass of a number of elementary particles would 
be better in this respect (|Becker et al. 2007[ 



2011 ILeonard 2010; Milton et al. 2007D. Milton et al. 



Hill et al. 



(120101) studies two alternatives, fixing either Na and h 
or Na and the atomic mass constant m u = M u /Na, and 
shows that this choice has little effect on the uncertain- 
ties of fundamental constants, mainly because the ratio 
h/m u is known well from the atomic recoil experiments. 



D. Quantum metrology triangle 



Phase-locked Bloch ( |Averi n et al. 
( jAverin and Likharev 



1985) and SET 



1986) oscillations in supercon- 



ducting and normal-state tunnel junctions, respectively, 
were proposed as a source of quantized electric current 
in the mid 1980s, soon after the discovery of the QHR. 
Already Likharev and Zorin (19851 suggested that the 



quantum current standard could provide a consistency 
check for the existing two electric quantum standards in 
an experiment they named as "quantum metrology tri- 
angle" . However, the quantized current turned out to 
be a much greater challenge than JVS and QHR. Still, 
after a quarter of a century, quantum current standards 
are yet to take their place in metrology. On the other 
hand, the progress on the knowledge on Kj and Rk has 
also been rather slow: in CODATA-86 the u ncertainties 
were 300 x 10 -9 and 45 x 10 -9 , respectively (Cohen and 



Taylor 1987). These uncertainties are essentially on the 
same level as in CODATA-10 if the discrepancy of the 
data is taken into account. 

The QMT experiment and its impact has been dis- 



cussed, e.g., in 


Felt in and Piquema 


( 


2009); Gallop 


(2005 


; Keller ( 


2008); Keller et al. 


( 


2008 


1; |Piquemal 


(2004 


2012J) ; Piquemal and Geneves 


200C 


): 


Scherer and 


Camarota (|2012 


)■ 


Zimmerman and Keller 


(2003). There 



are essentially two strategies to close the QMT: by ap- 
plying the Ohm's law V = RI or by so-called electron 
counting capacitance standard (ECCS) which utilizes the 
definition of capacitance C = Q/V. They are sometimes 
called direct and indirect QMT, respectively. Regardless 
of the variant, QMT compares the macroscopic current 



generated by a single-electron source to the other stan- 
dards. 



1. Triangle by Ohm's law 

Applying the Ohm's law is the most obvious way to 
compare the three quantum electric standards. It can be 
realized either as a voltage balance Vj — RhIs °r as a 
current balance Vj / Rh — Is- In both cases, substituting 
Eqs. (|40l to Vj = R H Is yields 



njiK fj 
2k s X Is 



1 + ej + ^ + e s . 



(44) 



The major difficulty in the realization of the QMT can 
be outlined as follows. Let us consider the ideal case 
where the noise of the experiment is dominated by the 
Johnson noise of the resistor. The relative standard de- 
viation of the measurement result is 



Sis 
Is 



tRI 2 s ' 



(45) 



By substituting realistic estimates t = 24 h and T = 
100 mK for the averaging time and the temperature of 
the resistor, respectively, and by assuming that R — 
Rk/2 and 7g = 100 pA, one obtains the uncertainty 
SIs/Is ~ 7 x 10~ 7 . In practical experiments, l/f noise 
and noise of the null detection circuit make the mea- 
surement even more demanding, but this simple model 
demonstrates that the magnitude of the current should 
be at least 100 pA. 

Another problem is that the product RhIs yields a 
very small voltage, e.g., 12.9 Ml x 100 pA = 1.29 ^V. 3 
Even the voltage of a JVS with only one junction is 
typically of the order of 70 GHz/ Kj « 140 fiV. Such 
low voltages are also vulnerable to thermoelectric effects. 
One way to overcome this problem is to multiply the 
small current of the SET by a CCC with a very high 
winding rat io , ^10,000, as was suggested i n |Hartland 



et al. 



(1991); Piquemal and Geneves (2000); Sese et al. 



( |1999[ ), see Fig. |38| (a). It allows room-temperature de- 
tection, and that JVS, SET, and QHR can be operated 
in different refrigerators. This type of efforts have been 
descri b ed in Feltin and Piquemal (2009); Feltin et al. 
(2011); Piquemal (2004). Another approach is to use 
a high-value cryogenic resistor that is calibrated against 
QHR with the help of a CCC flElmquist et al. I [20031 pan" 



ninen et al. 2010[ ). All parts of the Ohm's law are in the 



same cryostat which can reduce thermoelectric effects, 
see Fig. 38 'b). Only the difference current Vj/R clyo — Is 
needs to De multiplied. Despite of persistent efforts, these 
approaches have not yet produced any significant results. 



3 A quantum voltage standard based on integr ating a semiconduc t 
ing pump with QHR has been pioneered in Hohls et al. (2011 1. 
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(d) A 




FIG. 38 (a-c) Variants of Ohm's law triangles where the 
quantized current (Is) is compared to resistance (R) cali- 
brated against QHR and to JVS. (a) Quantized current is 
magnified by a CCC, which allows room temperature null de- 
tection of the voltage difference (AT/), (b) Triangle with a 
high- value cryogenic resistor. The current balance Al can be 
determined, e.g., with the help of a CCC. (c) QMT experi- 
ment where the null detection is performed by a room tem- 
perature transimpedance amplifier, (d) ECCS experiment. In 
the first phase (A), the electron pump charges the cryocapac- 
itor C « 2 pF. An SET electrometer (E) is used to generate a 
feedback voltage ( V) that maintains the potential of the island 
at zero. Hence all the charge is accumulated to the cryocapac- 
itor and not to the stray capacitance. The feedback voltage 
constitutes the third part of the indirect QMT Q = CV. In 
the second phase (B), the cryocapacitor is calibrated against 
the re ference C re f which i s traceable to a calculable capacitor. 
From |Keller et al] |l999| . 



The best result from the Ohm's law triangle so far were 
obtained recently in Giblin et al. ( |2012| ) where a CCC 
with high winding ratio was used to calibrate a precision 
1 Gil room temperature resistor, see Fig. 



3S 



This 

experiment benefited from the relatively large current of 
150 pA that was generated by a semiconducting quantum 
dot pump. The uncertainty of the QMT experiment was 
1.2 x 1(T 6 . 



2. Electron counting capacitance standard 



The ECCS experiment was first suggested in Williams 
et al. (1992). A single-electron current source is used to 
charge a cryogenic capacitor C cryo by a known number 
Ns of electrons. The generated voltage is compared to 
JVS. The result 



C c: 



N s Qs 
V 



(46) 



thus yields a quantum capacit ance standard. The E CCS 
exp erim ent was pioneered in Keller et al. (1999), see 
Fig. 38 d), where an uncertainty of 0.3 x 10 _D was ob- 
tained for the ECCS capacitance. In this approach, 
ECCS was compared to a calculable capacitor. Then 



the observational equation corresponding to Eq. ( 44 ) 
HocnjfjC 



AaN s 



! + £./ + es, 



(47) 



does not include ex- However, calculable capacitors have 
been compared to QHR with very low uncertainty, and 



ac QH R techniques (Camarota et al. 2012 Schurr et al. 
2011) allow C cryo to be compared directly against Rk 
One should thus obtain an uncertainty of ~ 10~ 8 before 
there is any significant difference between the implica- 
tions of the two QMT versions. 

A major weakness of the ECCS is that it calibrates 
C cryo at ~ 0.01 Hz, but commercial capacitance bridges 
that are used to compare C cry o to the calcula ble capacitor 
(and also a c QHR) operate at 



Zimmerman 



1000 Hz 

et al. (2006) presents a model for the dielectric dispersion 
of insulating films at the surface of the electrodes of the 
capacitor. They fit this model to measurements of the 
frequency dependence and its temperature dependence 
in the range 100-3000 Hz and 4-300 K. The frequency 
dependence decreases at low temperatures. They evalu- 
ate that it yields an uncertainty component of 2 x 10~ 7 
for the QMT. In |Keller et cd] ( |2007[ ) this estimate was 
used to finish the uncertainty budget of the NIST ECCS 
experiment that closes the QMT at the uncertainty of 
0.9 x 10~ 6 . 

Very recently, PTB reached the uncertainty of 1.7 x 
10~ 6 in the ECCS experiment (Camarota et al, 2012). 



They used a five-junction R-pump instead of the seven- 
junction pump without resistors employed by NIST, 
but otherwise the experiments are quite similar. PTB 
presents their result as "preliminary" and plan both a 
more detailed uncertainty budget and several improve- 
ments to the experiment. Besides NIST and PTB, the 



ECC S has been pursued at METAS (Riifenacht et al. 
20l0b. 



3. Metrological implications of single-electron transport and 
QMT 

It is clear that when a quantum current standard will 
be available, it must be compared against other stan- 
dards in the QMT. One should, however, be conservative 
when interpreting the results. The single-electron sources 
are proven much more sensitive to failures than JVS or 
QHR, and unless the pumping errors are detected sepa- 
rately, QMT should only be considered as a way to char- 
acterize the current source. This is the case for Giblin 



et al. ( 2012 ) . Only the combination of error counting and 
QMT allows to separate (kg) from ej.K.s- So far, it has 
been achieved with a reasonable uncertainty (~ 10 -6 ) 
only in the EC CS experiments of NIS T ( |Keller et al] 
2007[ ) and PTB ( [Camarota et aL||2012 ). 

As shown in Sec. IV. B an uncertainty of < 2 x 10 -8 is 



required to give information on ck, and an uncertainty 
of the order of 10~ 7 would strengthen the knowledge on 
ej. Thus, the NIST and PTB results, e s = (-0.10 ± 
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0.92) x 1CT 6 and eg = (-0.3 ± 1.7) x 10~ 6 , respectively, 
only measure ej ( |Keller[[2008| . 



V. PERSPECTIVES AND OTHER APPLICATIONS 



Milton et al. (2010) analyze a scenario where ej is an 
adjusted parameter and eg = ex = 0. They study the 
effect of QMT on the uncertainties of fundamental con- 
stants. It turns out that when QMT is inaccurate, the 
uncertainties of h, e and m u are mainly determined by 
the Avogadro experiment, but their uncertainties will be 
dominated by those of the watt balance and the direct 
measurement of Rk when QMT is improved. 

One problem of the QMT is that it only gives a value 
for the sum of the errors of the quantum standards, and in 
principle, they could cancel each other. It is thus useful 
to have independent test for each standard, and those 
for JVS and QHR are discussed in Sec. IV. B A test for 
the current standard only, i.e., an SI value for Qs, can 
be obtained by combining results from three experiments: 
QMT, measurement of Rk by a calculable capacitor, and 
watt balance ( |Keller et aT| |2008l ). Applying Eqs. (42) and 
(46), and substituting R by l/(wC), one obtains 



Qs = 



N S 



mgvC 



(48) 



Also the Ohm's law triangle can be used to yield a similar 
result, but in a less direct way. One should note, that 
JVS and QHR are used here only as transfer standards. 
Keller et al.\ ( |2008[ ) derives a result based on the NIST 
ECCS: Q s = 1.6021763 x 10~ 19 ± 1.5 x 10~ 25 C. This 
could be compared to the CODATA value on e, which, 
however, depends strongly on h and on the exactness of 
Kj and Rk- Instead, Qs can be compared to another 
v alue of e that is also indep endent of JVS and QHR. e = 
v / a 3 A r (e)M u /(iJ, R oc> N A ) ( Feltin and Piquemal| |2009| , 
whose uncertainty ~ 1.5 x 10 -8 is dominated by that of 
N A . Using the NIST ECCS result and the IAC value on 
N A , one obtains e s = (-0.2 ± 0.9) x 10~ 6 . 



We note that according to Eq. (|47|, QMT also yields 
a value for a independently of QHR, which was one of 
the early motivations for ECCS (Williams et al. [1992). 
This fact, however, has little importance before the un- 
certainty is competitive with the atomic recoil experi- 
ments (< 10 -9 ). Then QMT would strengthen the veri- 
fication of QED. 

Although single-electron transport would be the con- 
ceptually most straightforward realization of the ampere 
in the future SI, it is not likely that it would replace JVS 
and QHR as the typical realization in the near future. 
The exception is naturally the growing field of metrology 
for small electric currents, where single electronics is ex- 
pected to yield major improvements of accuracy. On the 
other hand, when the accuracy of single-electron trans- 
port accuracy improves, it can yield vital information on 
other standards and on fundamental constants. 



The quantum dot pump dGiblin et al] |2012[ |K"aestner| 



et al. 2008b) discussed in Sec. 



III.C 



has definitely proven 
its potential to be the basis of the future quantum stan- 
dard of ampere. The verified uncertainty of the 150 pA 
output current on the level of 1 ppm and the theoretically 
predicted 0.01 ppm uncertainty of the present device are 
truly remarkable figures of merit. On the other hand, 
a few important questions remain to be answered be- 
fore one would realize the ampere with the quantum dot 
pump: The superior device performance depends crit- 
ically on applying a strong > 10 T magnetic field on 
it. This dependence is not fully understood, and the ex- 
act magnetic field characteristics seem sample dependent. 
The reproducibility of the highly accurate pumping re- 
sults with samples from different fabrication runs remains 
to be shown. Importantly, error counting experiments on 
the dot samples have not been carried out, which also 
prevents to study possible errors of other quantum stan- 
dards in the QMT. Future experiments will likely show 
whether all the relevant error processes have been ac- 
counted for in predicting the obtainable accuracy to be 
on the level of 10~ 8 . However, even if not in the case 
of a bare device, the quantum dot pump may perhaps 
be applicable to the realization of ampere, if the error 



correction techniques that were described in Sec. III.H 
become feasible experimentally. 

Another important development and potential future 
realization of ampere is the SINIS turnstile introduced 



in Sec. III.B Although presently inferior to the quantum 
dot pump in the level of current output, and consequently 
with less definite assessment of proven accuracy (present 
verified uncertainty below 10 -4 ), this device does not 
suffer from known obstacles on the way of achieving the 
required accuracy. Currently, the main error mecha- 
nisms have been assessed theoretically and experimen- 
tally, including photon assisted tunneling, Andreev cur- 
rent, co-tunneling, residual and generated quasiparticles, 
and possible residual density of states in a superconduc- 
tor. Positive conclusions can be drawn from individual 
experiments with respect to suppressing them in an op- 
timized device. Sample fabrication and reproducibility is 
currently on a high level, and it has been demonstrated 
that the requested magnitude of current can be achieved 
by running many turnstiles in parallel. For the SINIS 
turnstile, as for the quantum dot pump, the ultimate test 
would be an error counting experiment and the quantum 
metrological triangle. At the time of writing this review, 
such experiments have not been performed for either of 
these promising new single-electron sources. 

Presently several other new proposals are being pushed 
towards critical tests to study their applicability in cur- 
rent metrology: these include superconducting phase-slip 
wires, Josephson junction arrays, and mechanical shut- 
tles, just to mention a few less conventional ideas. Al- 
though it is not in the horizon at present, it is pos- 
sible that eventually one of these devices beats the 
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present Coulomb-blockade-based realizations both in cur- 
rent yield and in their robustness against transfer errors. 

Developing ever more accurate current sources has con- 
stantly been a driving force for understanding the under- 
lying physical phenomena. On the other hand, the stud- 
ies for the precise control of single electrons and Cooper 
pairs have created special expertise that is also applicable 
in variety of other research topics. 

In addition to the charge degree of freedom, the elec- 
trons hold information in their spin states which have 



were very recently tested by Saira et al. (Saira et al. 



been envisi oned (|Hollenberg et al. 2006 Kane 1998 ) to 



be utilized (|Morello et al 
processing 



2010 ) for quantum information 
Although the electron transport is typically 



incoherent in the electron pumps, the spin-encoded infor- 
mation can potentially remain coherent, and hence this 
information can possibly be transported from the mem- 
ory cell of the computer to the qubit-qubit interaction 
cell and back. The transport cycle has to be carried out 
with high accuracy for fault tolerant computing to be 
possible, which creates a close connection to the metro- 
logical electron pumps. 

Geometric phases ( |Shapere and Wilczek (edsjj 1989) 
in quantum mechanics have been studied extensively due 
to both fundamental scientific curiosity and their appli- 
cations in geometric quantum computing (|Zanardi and| 



the Berry, phase has already been measured in the super- 



conducting sluice pump (Mottonen et al. 2006 Mottonen 



et al. , 2008 ) thanks to the development of the sluice for 
metrology. Some theoretical work on the more complex 
phases referred to as holonomies has been put forward 



in the framework of Cooper pair pumps ( Pirkkalainen 
et al. 2010 Solinas et al. 2010) but it remains to be 



seen if these ideas will be implemented experimentally. 
The main obstacle in practice is perhaps the high level of 
precision required for the control signals of the pumps, a 
problem that can possibly be solved with the help of the 
work on the metrological current source. 

Detecting single-electrons and Cooper pairs by single- 
electron transistors and quantum point contacts has 
been largely motivated by the need for tests of the 
charge-transport errors in metrology. During the past 
decade, these techniques have also been successfully im- 
plemented, e.g., in experiments on full counting statistics 
and noise of charge transport. The experiments on the 
full counting statistics of current fluctuations in a semi- 
conductor quantum dot by real-time detection of single 
electron tunneling with a quantum point contact have 
been successfully performed for instance by Gustavsson 
et al. (Gustavsson et al. 2007 2006). In these experi- 



ments, moments of current up to the fifth one and be- 
yond could be reliably measured. Very recently, single- 
charge counting experiments have been applied to study 
energy fluctuation relations (Averin and Pekola 2011 



Crooks 1999 Evans et al. 1993 Jarzynski 1997) in sta- 



tistical mechanics. The experiments in stea dy-state non- 
equilibrium were performed by Kiing et al. ( Kung et al. 
2012), and the celebrated Jarzynski and Crooks relations 



2012b). Single-charge counting experiments allow one 



to test fundamental statistical mechanics and thermody- 
namics of classical and quantum systems. 

The variety of spin-offs out of the development of 
single-charge current sources for metrology is certainly 
expanding. This way the benefits of this research will be 
obvious not only for the community interested in the sys- 
tem of units and in traceable measurements, but also for 
other researchers working in basic and applied sciences 
looking for new tools for measurements that need precise 
control. 
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